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Die Integralgleichungen der Kippung gerader Trager 
mit diinnwandigen, offenen und doppeltsymmetrischen Profilen 


Von L. Kraus 


1. Einleitung. Die Erforschung der Mechanik des Kipp-Problems bei geraden Tragern mit 
diinnwandigen, offenen und doppeltsymmetrischen Profilen, jenes Sonderfalles des Biegedrill- 
knickens, kann man heute als abgeschlossen betrachten. Seit dem Erscheinen der Dissertation von 
L. Prandtl’, in welcher zum ersten Male des Auskippen von Tragern behandelt wurde, entstand 
eine Fille von Literatur”. Eine vollstandige Darstellung des Knickens, Biegedrillknickens und 
Kippens mit zahlreichen Literaturangaben ist wohl in dem neuen Lehrbuch von Kollbrunner und 
Meister* zu finden. Abgesehen von einfachen Sonderfallen sind die mathematischen Schwierig- 
keiten bei der Lésung von Kipp-Problemen naturgema8 sehr groB. Auch sind hierbei das Ritz- 
und das Differenzen-Verfahren im allgemeinen recht miihevoll. Deshalb liegen zahlenmaBige Ergeb- 
nisse bei den oben erwahnten Veréffentlichungen nicht fir alle in der Praxis auftretenden Fille vor. 

Im folgenden wird daher ein Weg gezeigt, auf dem man mit ertraglichem Rechenaufwand, 
wenn auch nur naherungsweise, so doch mit beliebig steigerungsfahiger Scharfe, kritische Kippwerte 
errechnen kann. Fiir den in zwei Gabeln gelagerten, geraden Trager und fiir den geraden Krag- 
trager mit diinnwandigem, offenem, doppeltsymmetrischem und iiber die Trigerlange konstant 
bleibendem Profil werden bei allgemeinster, stetiger Momentenverteilung die Integralgleichungen 
mitgeteilt. Die Integralgleichungen fiir den Sonderfall des Kragtragers, bei dem die Kipplast im 
Schwerpunkt des Endquerschnittes angreift, stellte E. 4. Deuker® auf. In der vorliegenden Arbeit 
gilt zuerst einschrankend, daB die 4uBeren Lasten ihre Angriffspunkte in der Schwerlinie des Tra- 
gers haben. Dann wird diese Einschrankung fallengelassen. Das heiBt, die Wirkungslinie der auf 
den Trager wirkenden, stetig verlaufenden Streckenlast q besitzt einen Abstand e von der Trager- 
schwerlinie. Aus den zuerst unter der Einschrankung gewonnenen Integralgleichungen erhalt man 
die Integralgleichung fiir e 4 0 durch Umformung des Kernes. Wie immer bei Kippaufgaben, so 
ist auch in den folgenden Untersuchungen die Wirkungsebene der auBeren Krafte gleich der aus 
Schwer- und Minimumstragheitsachse gebildeten Ebene. 


2. Die Grundgleichung. Die selbstadjungierte, einfache Differentialgleichung vierter Ordnung, 
welche den Kippvorgang eines Tragers mit den oben aufgezahlten mechanischen und geometrischen 
Eigenschaften und der oben gemachten Einschrankung beschreibt, Jautet in den dimensionslosen 
Koordinaten @ und s folgendermafen® 


1 a g 
— az Oasas + O45 +A M%s)d = 0. (1) 
Dabei bedeuten 
2 BP 

ae VEG, 
den dimensionslosen ,,W6lbkennwert“ des Tragerprofiles, 

ea atat the 

eras SED. 


den dimensionslosen Eigenwert der Differentialgleichung, | die Tragerlange [cm], D = G Jr die 
Torsionssteifigkeit nach de St. Venant [kpcm?], C, den auf den Schwerpunkt bezogenen Wélbwider- 


1 L. Prandtl, Kipperscheinungen, Diss. Miinchen 1899. 

2 §. Timoshenko, Theory of Elastic Stability S. 239—286. New York 1936. 

3 E. Chwalla, Forschh. a. d. Geb. Stahlbau Nr. 2. Berlin 1939. 

4 C, Kollbrunner u. M. Meister, Knicken. S.97—147. Berlin/Gottingen/Heidelberg 1955. 
5 E. A. Deuker, Ing.-Arch. 21 (1953) S. 399. 

6 Siehe FuRnote 3 S. 32, Gl. (D 10). 
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stand des Tragerprofiles [em®], B= E Jin die Biegesteifigkeit der Minimumstragheitsachse 


[kpem?], « die Langenkoordinate [cm] zur Beschreibung des Abstandes vom linken Tragerende 
(die x-Achse ist mit der Schwerachse des Tragers identisch), s = x/l die dimensionslose, auf die — 
Lange des Tragers bezogene Koordinate, 9 den Drillwinkel, M* die ausgezeichnete Ordinate im 

Momentendiagramm [kpcm], z. B. q/?/8 beim Trager auf 2 Stiitzen mit Gleichstreckenlast q, _ 
M(s) den qualitativen, dimensionslosen Verlauf des Momentendiagrammes, z. B. 4 (x/l) (L—x/l) | 
—4s(1—s) beim Trager auf zwei Stiitzen mit Gleichstreckenlast q, und Striche bedeuten Ab- — 


lJeitungen nach s. 


3. Die Integralgleichung fiir den Rechteckquerschnitt. Fiir Rechteckprofile nimmt der Wélb- | 


widerstand C, den Wert 0 an. Der Wélbkennwert a wird unendlich groB, und das Glied mit @,,,, 
in (1) verschwindet. Gleichung (1) reduziert sich somit auf eine Differentialgleichung zweiter 
Ordnung: 


8,, +A MXs) 8 =0. (2) 
Wir tberfithren nun (2) fiir verschiedene Randbedingungen in Integralgleichungen. Zu diesem 


Zwecke denken wir uns in (2) den mit # behafteten Term auf die rechte Seite geschafft und mit f(s) 
bezeichnet. Die verbleibende linke Seite nennen wir L[?]. Es gilt also+ 


f(s) =—AM%s) 0, 
Lid) =o (3) 
1{O)] = f(s) 


Gleichung (3) kann man somit als eine inhomogene Differentialgleichung zweiter Ordnung mit _ 


konstanten Koeffizienten auffassen. Die Lésungen u, und u, der homogenen Differentialgleichung (3) 
Jauten beim Fortlassen der Integrationskonstanten: 


ao | () 


Aus den n berechneten Liésungen,u; einer homogenen Differentialgleichung n-ter Ordnung L[?] =0 
kann man bekanntlich durch Variation der Konstanten nach Lagrange die Lésungen der inhomo- 
genen Differentialgleichung L[?] = f(s) wie folgt gewinnen. 


Man macht einen Ansatz: 


nt 
VON De 
i=1 
In allen Ableitungen von # auBer der n-ten, setzt man die mit v; behafteten Glieder Null und 
arbeitet sie in die inhomogene Differentialgleichung ein. Damit erhalt man ein lineares Gleichungs- 
system in v;, woraus man die v; berechnen kann und dann durch Integration die v;. Die aus dieser 
Integration herriihrenden n Integrationskonstanten sind mit den Integrationskonstanten der voll- 
standigen Differentialgleichung identisch. 


In unserem Falle errechnen wir vj und v; zu 
v, =f(s), 
v =—sf(s), 


daraus durch Integration: 


v,(s) = J f(t) dt + C,, 


v(s) = —ftf(t) dt + C,, 
f 
und schlieBlich die vollstandige Lésung der inhomogenen Differentialgleichung (3) 
2 s s 
Hs) = as v, =f sf dt + Crs—fltf() de +1. (7) 
= 0 0 
Bedenken wir, daB 


aes f(t) =—/A M(t) 0) ; (7a) 
' F. Schwank, Randwertprobleme, S. 232—324, Leipzig 1951. 
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dann folgt aus (7) und (7a) die Integralgleichung fiir (2): 
Ws) = —ifs M(t) P(t) dt + Aft M*(t) O(t) dt ++ Cys + C,. (8) 
0 0 


C, und C, berechnen sich aus den Randbedingungen. Diese lauten fiir einen Trager auf zwei Stiitzen 
mit Gabellagerung! 


? 
0 ux vs ae, 0) 


und fiir einen Kragtrager mit einer Einzellast P im Abstand ep iiber dem Schwerpunkt des End- 
querschnittes, wobei die Abszisse s vom freien Ende (Lastende) aus gezihlt ist, 


0=O0" fur tell 


o —Oetan os —— le | 


= 10 
# +Vlad=0 fiir sao, | eo) 


In (10) ist der Eigenwert A definiert durch \A=P P/V B D und « gegeben durch « = (ep/1)/ B/D. 
Eine eventuell zusatzlich auf den Kragtrager wirkende Streckenlast q miiBte bei der Bestim- 
mung der Kippstabilitat in Abhangigkeit von P ausgedriickt werden (z. B. g = konst. = y* P/l). 
a) Der Trager auf zwei Stiitzen. Wir betrachten zuerst den Trager auf zwei Stiitzen. Wir 
erhalten aus (8) und (9) 
C= 0 


und 


C= a M2(t) H(t) dt —A f t M%(t) B(t) de. 
0 0 


Setzen wir nun die errechneten Werte fiir C, und C, in (8) ein und spalten die Integrale auf in 


1 s 1 
fapts. 
0 0 s 
dann erhalten wir 
8 piel 
Hs) =2 ft (1—s) M%(t) H(t) dt +2 f s (1 —2t) M%X(t) H(t) de. (8a) 
0 s 
Multiplizieren wir nun (8a) mit M(s) und definieren eine neue Funktion 
g(s) = M(s) Hs) ; (8b) 
dann formt sich (8a) um in 
xis Mt 7 
p(s) =A | ft (1 —s) M(s) M(t) y(t) dt + f s(1—t) M(t) M(s) y(t) dt} . (8c) 
0 s 


Die Integralgleichung fiir (2) mit den Randbedingungen (9) lautet somit 


ie, 
p(s) =A f K(s, #) ge) at 
mit (11) 
K(s, t) =t M(t) (1— s) M(s) fir oStSs, 
K(s, t) = s M(s) (1—t) M(t) fir sStsl. 
Fir den Sonderfall, daB das Moment iiber den Trager konstant ist, gilt 
Ms) ==21., 
Somit geht (2) in eine Differentialgleichung iiber, die in der Literatur als Knickgleichung fiir den 
an beiden Enden gelenkig gelagerten Stab bekannt ist 
o,, +19 =0. 


Setzen wir in (11) 


1 Siehe S. 61 und S. 58 der in FuRnote 3 von Seite 1 genannten Arbeit. 
De 
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dann ist sofort zu erkennen, daB sich der Kern in den von Kléppel und Lie aufgestellten Kern des _ 


beiderseitig gelagerten Knickstabes verwandelt!. Damit haben wir eine wertvolle Kontrolle fir | 


die Richtigkeit von (11) gewonnen. 
b) Der Kragtraiger. Beim Kragtrager verlauft die Aufstellung der Integralgleichung analog. 
® ergibt sich zu 


2 s 
O'(s) = Duz,y,=1lff@d+1¢,. (12) 
i=I 6 


Gleichungen (7), (10) und (12) liefern 


A gente funy Wf d 


Gees | (¢—1) f(y) de. 
lee toe 
0 


Mit den errechneten Werten fiir C, und C, gehen wir in (7) unter Beachtung von (7a) und (8b) 
und erhalten : 


mit 


K(s,1) > —(1— Rat) (8) MG) M@) fir OStSs, (13) 


K(s, 1) = —— (1 = ins) (aM) MG) et 
1—Viu 


Fir M(s) = s, M(t) = t und « = 0 geht der Kern aus (13) in den unabhangig von uns aufgestellten | 
Kern (3.25) der Deukerschen Arbeit iiber. Durch diese Betrachtung ist die Richtigkeit von (13) | 


wohl auBer Zweifel gestellt. 


4. Die Integralgleichung fiir den I-Querschnitt. Bei I-Profilen ist der Wélbwiderstand von 
Null verschieden. Somit gilt fiir das reziproke Quadrat des Wolbkennwertes der Bereich 


zs EE <9) 
a 
und das Glied mit 1/a? verschwindet nicht in (1). Wiederum berechnen wir die Integralgleichung 
fiir die vollstandige Gleichung (1) nach der Methode von Lagrange. Es gilt 
i 
Lid} re qi vesss Uses 
f(s) = +2 MX) 8, 
L{O] = f(s) » 
4 
Maxie STH Ute 
i=l 
Die vier Lésungen u,; der Differentialgleichung L[?] = 0 lauten 
yet Oe, =O OS. 
US Oe ls 


Analog zu den in Abschnitt 3 dargestellten Operationen erhalten wir die Lésungen fiir das in v; 
lineare Gleichungssystem 


, 1 : 
v, = darn Seer 


a 


a K. Kléppel u. K.-H. Lie, Stahlbau 16 (1943) S. 17; dort Gl. (14) und (17). 


| 
| 


| 


(14) | 
| 
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v1, =—fi(s), 
v= +s f(s) 
und daraus durch Integration 
1 s 
w= +5 [ro edt + C,, 
t) 
1 , Y 
eo [r@ ett dr 4 G,, 
0 
es) ee Xe 
t= + Sef det G. 
a 
Somit ergibt sich fiir 9, 0, 0 und #’ nach Einfiihrung der Hyperbelfunktionen 
1 : mn 
H(s) = 2 Lae [ «Sin [a (s — t)]| —a(s —t)\ f(t) dt 
a 
+ Cet *! 4 CG, e-22 4. Gs + CL, (15a) 


eS) = > uv; = {Coy [a (s — t)] — 1} f(s) dt + C,aet+**— C,ae—** + Cy-1, (15b) 
0 


Ve) >t, 0, = a f¢Sin [a (s — t)]} f(t) dt + C, a® e+ ** + C,a®%e—*s, (15c) 
i 6 

OB Yu 1 = a [ oj [a (s — t)]} f(t) dt + C, a® et **— C, a e—*:. (15d) 
i @ 


Gleichung (15a) stellt im wesentlichen schon die vollstandige Integralgleichung des Kipp-Problems 
dar. Die vier Integrationskonstanten errechnen sich aus folgenden Randbedingungen ! 


1. bei einem in zwei Gabeln gelagerten Trager 
a) mit Wélbfreiheit an beiden Tragerenden 
P=2 
Vi 
b) mit Wolbfreiheit am linken und Wélbverhinderung am rechten Tragerende 
&=0 fir s=0 und s=1, 
0) erase =— ly (16b) 
i=) eatin si—= 0), 
c) mit Wélbverhinderung an beiden Enden 
Oo = | 
wv =of 
2. beim Kragtrager mit einer Last Pim Abstand ep itber dem Schwerpunkt des Endquerschnittes, 
wo & = (ep/L) B/D und 2 = (PP)2/B D 
a) mit Wolbfreiheit an beiden Enden 


1 Grit i 
OW + ae nae era ome 


fir $0 undss.== |, (16a) 


firs-=0. und -$: ==); (16c) 


i 3 | tticese——al es 


1 Siehe S. 34 Gl. (D 13), S. 33 Gl. (D 12) und S. 35 Gl. (D 22) der zitierten Arbeit von Chwaila. 
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b) mit Wélbfreiheit am freien Ende (s = 0) und Wolbverhinderung am Einspannende (s = 1): 


t 1 (hts aa ; 
—# +38 ces iceeans | 
oe” = 0 (16e) 
A) es 
er fiir s=1. 


a) Der Trager auf zwei Stitzen. Wir untersuchen im Folgenden den praktisch wichtigsten 
Fall beim Trager auf zwei Gabel-Lagern. Er ist durch die Randbedingungen (16a) gekennzeichnet 
und liefert die tiefsten Eigenwerte. 


Mit den sich aus (16a) ergebenden Integrationskonstanten 


G=—G=—scemay | Gin led —onsO a, 
0 


C, =f —a fe a, 
0 
C0 


und mit den Gleichungen (15a) und (8b) erhalten wir die Integralgleichung 
1 


g(s) =A f K(s, t) p(t) dt 


mit 


(17) 


Reg aaa a) Ms MG) fix 00 eee 


__, __ Sin [a(1 —1)] Sin [a s] s sa ee 
s(1— 1) eee | mc MG)" tar sores 
Fiir den Grenziibergang a—> oc geht (17) in (11) ttber. Dies erkennt man leicht, wenn man im 
Kern von (17) die hyperbolischen Kreisfunktionen durch e-Funktionen ausdriickt und Nenner und 
Zahler durch a e* teilt. Die Unterfalle (16b) und (16c) kénnen analog zum Unterfall (16a) berechnet | 


werden. 


b) Der Kragtrager. Fiir den Kragtrager mit den Randbedingungen (16d) ergeben sich die 
Integrationskonstanten zu 


1 


C, =— (= — ea | {Sin [a (1 — a} f@) de, 
1 1 : 
C, = — C, ae rere Se l i d 
Yaa (= Vaz [ HO a 


Daraus folgt mit (15a) und (8b) die Integralgleichung 


v(s) = K(s, t) p(t) dt, 


mit 


ee i. i = Sin [a(1 —s)] Sin [ot] 
K(s, t) eee eee a SE |M M 
ere vt) (ls) ae (s) M@ ng 
fir 6 0/Ssiaas. 
ices ee __) _. Sin [a — 9] Sin [a 5] 
(5) =| — a (== Viies) (ee) + Shea | M(t) M(s) 
fica ae 
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Beim Kragtrager mit den Randbedingungen (16e) erhalten wir die Integrationskonstanten 


nach Kinfithrung der HilfsgréBen P, Q und N zu 


Goat mes — al Sp ee Vi« 
c= — 1 ¢, = af? Sit fal — 0604 fa} 


fiw. at . 
mit 
P= [Go| fa (1—9]— Ysa, 
= J {Sin [a (1 —t)] —a (1 —0)} f(t) de, 
N =a {a(1— JZ) Gof [a] + Vix Gin [a)} 
Daraus folgt mit (15a) und (8b) nach langerer Rechnung die Integralgleichung 
g(s) =A f K(s, t) p(t) dt 
mit 
a (s,-t)i = Ke — Go} [a (1 — s)]} Sin [at] + Gin [as] — Gin [a] -+ a (1 — 5) Cof (a) ) 
+ Viw ¢ {60 [a (1 — s)] — = Gin [a (1 — s)] — s} Sin [at] 


+ #{Gin [«]—a(1—s) 60} [a] — Gin [asf] MONO __ 
<f Na (U2 a) Col [a] eV le Sit [a] (19) 
bie () 7k Se 


K(s, t) = Ku — Goj [a (1 — t)]} Gin [as] + Gin [at] — Gin [a] + a(1 — 2) Coj (ly 
si View ¢ {ooh [a (1 —1t)] — “ Cin [o(1— 1) i Gin [as] 
+ s{Gin [a] — a (1—t) Cj. [a] — Gin [at]} 5 Mi) M(s) 


Vara @ Vie Cu leleov ie Bin ial 
fir s<t<1.j 


Mit M(s) = s, M(t) =t und « = 0 stimmt (18) mit (3.23) der Deukerschen Arbeit iiberein. Fiir 
a-—> oo gehen (18) und (19) in (13) tiber. Da die Randbedingungen (16d) den geringsten Zwang 
ausiiben, liefert (18) die tiefsten Eigenwerte. Dieser Fall ist theoretischer Natur; denn in der 
Praxis wird man stets versuchen, durch konstruktive MaBnahmen die Flanschen im Einspann- 
querschnitt an der Verwélbung zu hindern, um dadurch die Kippsicherheit zu erhéhen [siehe (19) ]. 


5. Die erweiterte Grundgleichung. Nun wollen wir die Kerne der von uns aufgestellten Integral- 
gleichungen fiir die Falle abandern, bei’ denen die Streckenlast q nicht in der Schwerlinie angreift. 
In (1) muB infolge des oben Gesagten der Term mit # umgeformt werden; es gilt dann + 

— 1 

Se, es act MX) ae e ats) oO, (la) 

a Va 
Dabei bedeutet y** einen dimensionslosen Faktor, durch den eine ausgezeichnete Ordinate q* 
im Streckenlast-Diagramm mit dem fiir den Eigenwert A maBgebenden Biegemoment M* verbun- 
den ist, g* = (M*/I?) y**, und q(s) bedeutet den qualitativen, dimensionslosen Verlauf der 
Streckenlast (z. B. g(s) = 1 bei einer Gleichstreckenlast). Es besteht der Zusammenhang 

M* 
q=4" 4s) = y** als) - 

Ferner ist e* eine ausgezeichnete Ordinate [em] im Diagramm des Abstandes der Streckenlast 
von der Schwerlinie, e(s) der qualitative, dimensionslose Verlauf des Abstandes der Streckenlast 


1 Siehe S. 32 Gl. (D 10) der Arbeit von Chwalla. 
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von der Schwerlinie. Es gilt e = e* e(s). SchlieBlich bedeutet 6 = (e*/I) VB/D y** eine dimensions- 
lose HilfsgréBe. 


Es ist also in (1) anstelle von M*(s) der Term [:12(s) -f- (iva) 6 e(s) a(s)| getreten. Die Sym- 


metrisierung der entsprechenden Integralgleichungen erfolgt nun durch Multiplikation von O(s) mit 


/M%s) -}- (iva) B e(s) q(s) und Einfithrung von ¢(s) = 0s) Ms) + (7) f e(s) Rey Wir sehen 


daraus, daB sich die entsprechenden Kerne nur insofern verandern, als anstelle von M der Term 


ave + fe q/VA steht, wobei hier M, e und q von s baw. t abhangen. Ist e(s) = 0, dann gehen 
die abgednderten Kerne wieder in die urspriinglichen iiber. 

Alle Integralgleichungen (11), (13), (17), (18) und (19) wurden durch fortgesetzte Differen- 
tiationen kontrolliert, die sowohl die Randbedingungen erfiillen als auch die zugehérigen Dif- 
ferentialgleichungen liefern. 


6. Die naherungsweise Berechnung der kritischen Kipplasten. In vielen Fallen lassen sich die 
Eigenwerte der Differentialgleichungen (1) und (2) nicht exakt berechnen, weil es fiir die Dif- 
ferentialgleichungen keine einfachen Lésungen gibt. Der stets als Lésung fiir (1) und (2) geltende 
Potenzreihenansatz 

1S rete tees” 

Del) 
steigert die Arbeit bei der Auswertung der Kippdeterminanten im allgemeinen sehr, so da seine 
Anwendbarkeit relativ beschrankt ist. Bei Anwendung des Ritzschen Verfahrens zur Naherungs- 
berechnung der Kippwerte ist fiir jeden Momentenverlauf die Auswertung bestimmter Integrale 
von meist nicht ganz einfacher Bauart erforderlich. Wegen der groBen Rechenarbeit beim Auf- 
stellen der Determinantenkoeffizienten erscheint uns auch dieses Verfahren mit einem grofen 
Rechenaufwand verbunden und daher fiir numerische Berechnungen bei komplizierteren Rand- 
bedingungen weniger empfehlenswert. Dagegen liegt bei Verwendung der Integralgleichungs- 
methoden! die Formel fiir die Koeffizienten der spater noch zu erlauternden Matrix mit dem Vor- 
handensein der Integralgleichungskernes bereits vor. In den Fallen leicht integrierbarer Kerne ist 
es angebracht, die Erhard-Schmidtsche Ungleichung 


a ed Saiip: 
A> If [K(s, t)]? dt as} (20) 
0 0 J 

zu gebrauchen. Man gewinnt mit ihr, wie bekannt, eine untere Schranke fiir den kleinsten Eigen- 
wert. 


Das Operieren mit (20) ist dann zweckmaBig, wenn die Kerne (11) oder (13) gelten, sich das 


Biegemoment in Form von Potenzen von s darstellen la48t und der Eigenwert A nicht im Kern 
auftritt [also « = 0 in (13), (18) und (19)]. Bei den Integralgleichungen (17), (18) und (19) wird 
man (20) wegen Integrationsschwierigkeiten nur in Sonderfallen anwenden kénnen. Hier empfiehlt 
es sich, die Integralgleichung nach Fredholm als ein lineares homogenes Gleichungssystem der 
Stiitzstellen s; aufzufassen: 


IMs 


Hs) =A X Ks) pl) (/ = 01,2, p) (21) 


Mit 2 = p/A kann man (21) in Form eines speziellen Matrizeneigenwertproblems 

(K—AC) xr —0 (22) 
schreiben. Dabei sind die Komponenten des Vektors ¢ gleich den Naherungswerten der gesuchten 
Eigenfunktion g an den Stiitzstellen s;, und & ist die Einheitsmatrix und & die Kernmatrix. Das 
Verschwinden der Nennerdeterminante des homogenen Gleichungssystems (22) stellt die Kipp- 
bedingung dar. Die j-Werte lassen sich mit Hilfe des ,,Matrizen-Iterationsverfahrens“ oder des 
Hessenberg-Verfahrens? aus (22) ermitteln. Bei dem unten errechneten, ersten Beispiel zeigte sich, 
daB man alle Kigenwerte zum Beispiel mit Hilfe des Matrizen-Iterationsverfahrens finden kann, 
wenn man mit geeigneten Ansatzen fiir die Eigenfunktionen @ in die Rechnung eingeht. Ferner 
stellt hierbei der Rayleigh-Quotient keine obere Schranke dar, sondern nihert den Eigenwert von 
unten. Dies hat seinen Grund darin, daB (22) durch die Umwandlung eines Integrals in eine Summe 


1 Siehe FuBnote 1 von Seite 2. 
* R. Zurmiihl, Praktische Mathematik, Berlin/Gottingen/Heidelberg 1953. 
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entsteht. Hingegen fiihrt das Ritz-Verfahren (Energiemethode) auf das allgemeine Matrizen- 
Eigenwertproblem 

QU=-429) b=. (23) 
Dabei stellt B giinstigstenfalls eine Diagonalmatrix dar. Die Vorteile der Integralgleichungs- 
methoden gegeniiber dem Ritz-Verfahren fassen wir in den vier wesentlichen Tatsachen zusammen: 

1. Die Iteration des speziellen Matrizen-Eigenwertproblems ist mit weniger Rechenaufwand 
verbunden. 

2. Wie oben ausgefiihrt, ist mit Aufstellung des Kernes bereits die Koeffizientenformel fiir § 
gegeben. 

3. Man umgeht das beim Ritz-Verfahren notwendige Aufsuchen eines passenden, den tat- 
sachlichen Kurvenverlauf gut beschreibenden Naherungsansatzes, was bei komplizierteren Rand- 
bedingungen auf Schwierigkeiten stéBt. 

4. Fiir Integralgleichungen stehen Verfahren bereit, mit denen man die J-Werte nach oben 
und unten einschranken kann!. 

Tritt jedoch der Eigenwert im Kern der Integralgleichung auf, wie etwa in (13), (18) oder (19), 
versagen die oben geschilderten, nur auf lineare Eigenwertprobleme anwendbaren Naherungs- 
verfahren. Wir schlagen in diesem Falle das folgende Verfahren zur Ermittlung der Figenwerte 


vor: Fiir diskrete 2 bzw. A berechnet man die zugehérigen Werte AN der Nennerdeterminante von 
(22) und tragt sie in einem (4N,/)-Diagramm auf. Diejenigen /-Werte, fiir welche die zugehérigen 
AN-Werte die GréBe Null annehmen, gewinnt man naherungsweise durch Ablesen und scharfer 
durch lineare oder quadratische Interpolation. Die so gefundenen J-Werte stellen die naherungs- 
weise ermittelten Eigenwerte der Integralgleichung dar. 


7. Erstes Beispiel. AbschlieBend mége nun die Berechnung der Kigenwerte einiger einfacherer 
Falle die oben gewonnenen Erkenntnisse veranschaulichen. 


a) Untere Schranken fir die Kippwerte nach der Integralgleichungsmethode. Als 
erstes Beispiel berechnen wir nach (11) und (20) naherungsweise den tiefsten Eigenwert fiir einen 
Balken auf zwei Stiitzen mit Rechteckquer- M, 
schnitt, der durch zwei voneinander ver- 
schiedene Endmomente belast ist (Abb. 1). Gon oe 


Das Biegemoment ist 


le 


Der Kern hat die Bauart Abb. 1. Systemskizze fiir Beispiel 1. 


; >| |M. IWEN . 25 E 
IK (5, 2) = Is (! 2 i) s (! =| * Ine tS ( -¥) “| Hue W) <p SS Ge 


M,\ ,|[M ite ots hates 
K() =| (1 2) #—(1 mt) *| mes +(1 a)" fits re, 
Ferner ist 
= _ Mb 
= 5D 


Zieht man in der Integralgleichung (M,/M,)? aus dem Kern, dann ist der Kigenwert 


M,\2> _ /M,\2(M, 1) 
(mn) *= (wi) “aD 


und der Kern schreibt sich wie folgt: 
K(s,t) = +ms+ns*)¢—n?) fir 0StSs,) 
K(s, t) = (1 + mt + n?#*) (s —ns*) in OSes ting 


(24) 
mit Pant ‘ 
nt Ned z wa, bavel © 70. a) - 


M,/M, M,/M, a 


ine ahnli ic ; sei hier auf ‘as, Ei ehwingungen 
1 Als eine ahnliche Anwendung der Integralgleichungsmethode sei hier auf K. Karas, Eigense 4 
von Saiten mit elastisch befestigten Enden, Federhofer-Girkmann-Festschrift S. 37, Wien 1950; und Osterr. 


Ing.-Arch, 9 (1955) S. 352, verwiesen. 
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Das Doppelintegral wird nach langerer Rechnung 


et 1l{s 1 
ff [K(s,0)P dtds = an [K(s, t)]? dt + f LK(s, t)]? dt; ds 
00 0 46 3 
eral 
~ 6300 
— 900 m2? n + 140 n — 855 m n? + 315 m? ni? 
— 855 n* — 70 n® + 252 n4 + 560 m n°) . 


Daraus berechnet sich das kritische Moment M, zu 


(1050 + 1680 m -++ 700 m? — 900 mn (24a) 


_ /BD j 6300 
ig pe ee 


a (25) 


Yaw ema 


M, 
Der Inhalt der eckigen Klammer ist wegen seiner Lange in (25) nur angedeutet und mit dem der 
eckigen Klammer in (24a) identisch. 


Als Sonderfalle betrachten wir 


1. M, = M,, somit m = —1 und n = 0. Es ergibt sich M, = 3,0801 VBDIl. Der Fehler von 
3,0801 gegeniiber dem exakten Wert z betragt 1,96%. 
2. M, = 0. Unter der Wurzel in (25) verbleiben nur die Terme, die in m und n biquadratisch sind 


1/6300 
/315 + 252 — 560 


Wir bekommen M, = 5,4772 /BD/I. 
Der Fehler von 5,4772 gegeniiber dem 
exakten Wert 5,5618 betragt 1,52%. 

In Abb. 2 ist M,]///BD aus (25) tiber 
M,/M, aufgetragen. Fiir die Verhalt- 
nisse M,/M,=+1, +1/2, 0, —1/2 
und —1 wurden aus Vergleichsgriinden 
die Kippwerte gemaB (21) mit 4- und 
9-reihigen Matrizen berechnet. Die da- 
mit gleichzeitig erhaltenen Naherungen 
fiir den Verlauf des Drillwinkels sind in 
Abb. 3 dargestellt. 


b) Obere Schranken fiir die 
Kippwerte nach der Energie- 
methode. Zur Gewinnung — oberer 
Schranken fir die naherungsweise 
berechneten Kippwerte wollen wir nun 


SI SATIZE 


7 = : die Energiemethode heranziehen. Mit 
: den Abkiirzungen 
b pons M,/M, —_ 1 
My 1 — M,/M, i)” 

: | %@  =M = M,(1— M,/M,) (s +b),4 (26) 
“4 -$ 0 +4 +7 1 

Abb, 2. Beispiel 1, Kippwertverlauf. Y io BD [M, l (1 oss M,/M,))? 2 

rata SpE ena i rainy schreibt sich die Differentialgleichung (2) 


O15 PAs bP BG =O. (27) 
Durch Multiplizieren von (27) mit 9/2 und nachfolgendes partielles Integrieren in den Grenzen 
von 0 bis 1 unter Beachtung von #0) = 8a) = 0 erhalten wir bekanntlich den Energieausdruck JJ zu 


a | [(9,)2 —A (s + 6)? #2] ds = 0. 


0 


| Me 
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Nach Einfithrung des die Randbedingungen befriedigenden Naherungsansatzes 
Ks) = >) asin (mz s) 


driickt sich das Verschwinden der ersten Variation von // folgendermaBen aus 


1 
ee wimor es » 9 00 
wi 2 sage >| 2. le Ce 
m m 0 
Daraus folgt beim Weglassen der Zwischenrechnungen das in a,, homogene Gleichungssystem 


ee a k* m?* | ay, + Di Bin Cia 0 > 


(n#m) 


at O0,,2= Or 


mit! 
6 a4 
vk 

ae A 

2m? x? — 3 2 \ (28) 

A,, Ss aE a 6 7? b (1 +- b) 5 = 
Bae = | A8 ae fiir m+ 1 gerade 9 
Brn = —48 aa (1+26) fiir m+n ungerade . 


Die Forderung, da die Nennerdeterminante von (28) verschwinden muB, liefert die obere 
Schranke fiir die Kippwerte. Die Rechnung zeigt, da8 bei M,/M, = 0 ein Naherungsansatz fiir 
# mit m = | und 2 Eigenwerte ergibt, die um weniger als 1° von den exakten Werten abweichen. 
Bei M,/M, = —1 zerfallt die Matrix (28) wegen (1-++2b) = 0 in zwei Teilmatrizen. -Diejenige 
mit ungeradem m liefert bei einem zur Achse s = 1/2 symmetrischen Drillwinkelverlauf den tiefsten 
Eigenwert (ein 3-gliedriger Ansatz ergibt eine Abweichung von 0,02°%). Die Teilmatrix mit geradem 
m gilt fiir den zur Achse s = 1/2 antimetrischen Drillwinkelverlauf. Ihr kleinster Eigenwert ist 
gréBer als der aus der ersten Teilmatrix gewonnene. Bei M,/M, = — 1/2 betragt der Fehler des 
aus einer 4-reihigen Matrix berechneten Naherungswertes gegeniiber dem genauen Wert 0,06%. 

c) Die strengen Kippwerte. Um ein Bild von der Giite der oben auf mehreren Wegen 
erhaltenen Naherungsergebnisse zu bekommen, mégen nun abschlieBend fiir die oben behandelten 
Verhaltnisse M,/M, die Kippwerte durch exakte Lésung der Differentialgleichung (27) berechnet 
werden. 

Nach Substitution von z = (s + b) wandelt sich (27) um in 


O, tAPO =0. (29) 

Fiir das Verhaltnis M,/M, = 1, somit M = konstant, ergibt die Lésung von (2) fiir J mit Be- 

achtung der Randbedingungen (9) sofort den trigonometrischen Sinus mit dem kleinsten Eigen- 

wert z, so dafi in diesem Falle nicht auf die Lésung durch Bessel-Funktionen zuriickgegriffen zu 
werden braucht?. 

Die Randbedingungen fiir (29) lauten gemaB (9) 
G0 \firr7, = .0% | 
Ofer = ee hoe | 


Fiir die oben notierten Verhaltnisse folgt mit Beachtung von (26) und (30) 


(30) 


1 
a) been cent k und 4, = 2 fir =z: ) 
VE 
b) Bo Pea O und 275: |. far i ee | | 
l 9 M. 1 c (30a) 
< es 1 
¢) eet UN 2, es x? 
1 Lee Ml, 
d) a b= unde, — 5 oe i ; 


1 Verwechslungen mit den in (24) eingefiihrten Abkiirzungen m und n sind hier wohl nicht zu befurchten, 
2 Siehe S. 24 Gl. (C 38) der Arbeit von Chwalla, wobei dort f = 0 und D = 0 gesetat werden miissen. 
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Die Lésung der Differentialgleichung (29) erhalt man mit Hilfe der Lommel-Transformation 


2) =e Ic Jus & 2) = OR aioe, ie 2). (31) 


wo C, und C, die aus den Randbedingungen (30) zu bestimmenden Integrationskonstanten dar- 
Hellen Das Finarbaten der Randbedingungen (30) liefert zwei homogene und lineare Gleichungen 
in C, und C,. Durch Nullsetzen ihrer Nennerdeterminante 


fee rape A re a 
| (Cy Yo Jia E jb J—14 ie v | 
ito Jal (1 + vp jib smP (1 +o] 
erhalten wir bekanntlich die Kippbedingung. Fiir M,/M, > 0 kann man die Wurzeln yo und 
y1 + 6 ohne weiteres aus (32) herausziehen. 
Die Kippbedingung lautet dann? 
Jia eo +) Jos (ee »p| — a dus E (1 + 6) Jay 


Dagegen erfordert der Bereich M,/M, < 0 eine gesonderte Betrachtung: Wir iibersehen den 


(32) 


Fe) ow 


Verlauf von #(z) sofort, wenn wir die ersten Glieder der mit /z multiplizierten Bessel-Funktionen 
von der Ordnung ++ 1/4 und — 1/4 in (31) explizit niederschreiben: 


ALB A318 
2) = G |g pleas da a 
O=G hares) * BPP TOA™ | | a 


CL? ieee | 
F218 (3/4) 277? P(2) TU * i 


In (31a) beschreibt die mit C, multiplizierte unendliche Reihe eine zur (z = 0)-Achse antimetrische 
und die mit C, Grduplemere tnendlicle Reihe eine zur (z = 0)-Achse symmetrische Funktion. 
Somit gilt 


y/ tulle | = al a “Syn (Gr 2, 


\—#J—1) i (— 2 =+y-aJiys ee 2} 


Deswegen miissen in (32) und in (32a) fiir den Bereich — 1 < M,/M, < 0 die ors in Klammern 
angedeuteten Minuszeichen stehen. 


Fir M,/M, = 0 folgt aus (30a) und (31a) 
OF = 0 > 


und der kleinste Kippwert ergibt sich nach Prandtl zum doppelten Wert der ersten Nullstelle der 
Bessel-Funktion von der Ordnung +-1/4 


Aes selon 
D 


Fiir M,/M, = — 1 liefern (30) und (30a) 
b? = (1 + b)*, 


und (32a) unter Beachtung des dort eingeklammerten Minuszeichens 


Jiys Z e Ay 2 | =0. (32b) 


Durch Nullsetzen von J 14 Ie i in (32b) erhalten wir den kleinsten Kippwert, und damit folgt 
aus dem (32) zugeordneten, in C, und C, homogenen Gleichungssystem C, = 0. Der kleinste 


? R. Rothe u. I. Szabo, Héhere Mathematik Teil VI, S. 143—169, Stuttgart 1953. 

* Das Auswerten von (32a) geschah mit Hilfe des Buches Tables of Besse] Functions of Fractional 
Order Vol. I, New York 1948. Dieses Tabellenwerk stellte das Institut fiir Praktische Mathematik (Prof. Dr. 
A. Walther), T. H. Darmstadt, dem Verfasser entgegenkommend zur Verfiigung. 
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Kippwert ergibt sich somit zum 4-fachen Wert der ersten Nullstelle der Bessel-Funktion von der 
Ordnung —1/4 | 


M, 1 
—————— 8,0252 . 

BD 
Die Zusammenstellung der oben auf verschiedenen Wegen erhaltenen tiefsten Kippwerte ist in 


der heigefiigten Tabelle 1 zu finden. Die eingeklammerten Zahlenwerte stellen die prozentuale 
Abweichung yom exakten Wert dar. 


Tabelle 1. Zusammenstellung der berechneten Kippwerte tal 


eS 
) MT 1 

M,/M, —l meies 0 a si ai 
a 


exakte Werte | 


gema (32a) ) 0252 =) 73241 5,5618 4,1207 | I 

Nah.-Werte ‘ | 7,8627 7,1100 5,4369 4.0490 | 3,0902 
mit p=5 gemaB (21) | (2,03) (2,92) | (2,25) (1,74) (1,64) 
Nah.-Werte ie Qo 7,2581 5,5288 4,1026 3,1287 


mit p=10 gemaB (21) (0,92) (0,90) | (0,59) | (0,44) | (0,41) 


untere Schranken | 17,4917 | 17,2246 5,4772 4,0437 3,0801 
gemaB (25) (6,65) | (1536)"-/)_ (1,52) (1,87) (1,96) 
obere Schranken 8,0265 | 7,3282 | 5,5804 4.1211 1 
gemaB (28) | (0,02) | (0,06) (0,33) (0,01) (0) 


d) Die Berechnung der Verformungen. Der genaue Drillwinkelverlauf fiir M,/M, ~ 1, 
—1 und 0 lautet mit der stets bei Eigenwertproblemen auftretenden willkirlichen Konstanten 


K=(CG 


c yA , 

ae I Sirs ie (1 + 6)? a 

Hz) = Ky23Siya 2 a & Sed a a . (31b) 
Jo 144 EB ly vy 


Fir M,/M, = | ist, wie nach (29) erlautert wurde, eine Halbwelle des trigonometrischen Sinus 
die erste Eigenfunktion. Fiir M,/M, = 0 wird die erste Eigenfunktion durch die ,,erste Halb- 
welle** der zur Achse z = 0 antimetrischen, in (3la) mit C, multiplizierten Funktion, hingegen fiir 
M,/M, = — 1 durch die ,,erste Halbwelle der zur Achse z = 0 symmetrischen, mit C, multipli- 
zierten Funktion dargestellt. 

Fir die untersuchten Verhaltnisse M,/M, sind die ersten Eigenfunktionen mit den zugehérigen 
Naherungen des Drillwinkelverlaufes in Abb. 3 aufgetragen. Der Unterschied zwischen den exakt 
berechneten Eigenfunktionen und den aus 9-reihigen Matrizen berechneten Naherungen ist bei 
der zur Verfiigung stehenden Zeichengenauigkeit nicht mehr erkenntlich. Die strichpunktierte 
und die ausgezogene Linie fallen daher in Abb. 3 zusammen. Die Durchbiegungen w in Richtung 
der Minimumtragheitsachse ergeben sich nach der Kipptheorie durch Integration von 


Wx x a (33) 


unter Beachtung der Randbedingungen w(0) = w(l) = 0 fiir unser Beispiel exakt zu 


i les (an) op + Wp, (33a) 
wo @p und , die Abkirzungen fiir die folgenden Funktionen bedeuten: 
Op = 2s sas os und @, = 's (l= 5%)’. 
Die Durchbiegungen wu in Richtung der Maximumtragheitsachse lassen sich aus der Formel? 
M(x) (2) Os 


u == i 
id FE Imin 


1 Siehe S. 9 GI. (A 23) der Arbeit von Chwalla. 
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mit u(0) = u(l) = 0 bestimmen. Dies geschieht am besten naiherungsweise nach den baustati- 
schen Methoden entweder graphisch durch Zeichnen eines geeignet verzerrten Seilpolygons fiir 
einen Trager auf 2 Stiitzen mit einer Streckenlast von der Intensitat w(x) = — M(x) O(x)/E Iinin 
oder analytisch mit Hilfe der Y&-Gewichte. Die Diagramme fiir w und wu des behandelten Bei- 
spieles sind der Vollstandigkeit und Anschauung wegen in Abb. 4 aufgezeichnet. 


Hs) 


Abb. 3. Beispiel 1, Drillwinkelverlauf, naherungsweise aus 4- und Abb. 4. Beismel 1, Verlauf von w exakt und von u naherungsweise 
9-reihigen Matrizen und exakt. mit ¥8-Gewichten. 
—_— — ® aus 4-reihiger Matrix gemaB (21) und (24); ni Berne 
9 aus 9-reihiger Matrix gemiB (21) und (24); peal OT ea gemaB (33a); 
—++—-++— # exakt gemaf (31). u gemafs (34). 


8. Zweites Beispiel. Als zweites Beispiel berechnen wir den kritischen Kippwert einer Last P, 
die am Ende eines Kragtragers mit Rechteckquerschnitt sitzt und zuerst im Schwerpunkt des 
Endquerschnittes angreift?. 


a) Untere Schranken fiir die Kippwerte nach der Integralgleichungsmethode. 


Es ist M(s) =s, 
und es gilt (13) mit x = 0. Das Doppelintegral des Kernquadrates in (20) hat den Wert 
Lae safes . 
| | [K(s, t)]2 dt ds = | | [(s—s®) ¢]2 dt + | [(t—#) s}? dtbds = ae 
d 0 6 bo 3 


Somit ist P = 3,9843 VBDIP gegeniiber 4,0126 VBDIP der exakten Lésung (0,71°% Fehler). 


1 Siehe FuBnote 1 von Seite 1. 
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Greift die Last P ober- oder unterhalb des Schwerpunktes des Endquerschnittes an, dann gilt 
(13) mit M(s) =sunda #0. Der Eigenwert J tritt nun auch im Kern der Integralgleichung auf. 
Fir verschiedene Werte « wurden die zugehdrigen kleinsten Higenwerte aus 4-reihigen Matrizen 
gemaB (22) nach dem oben vorgeschlagenen Verfahren berechnet und in Tabelle 2 zusammen- 
gestellt. Fiir « —> oo lautet der Kern von (13) 

K(s,t) =t(l—s)st fir 0 
K(s,t)=s(l1—it)ts fir s 


t 


t 


WA IIA 
IA IIA 


Hierbei, sowie in dem Falle « = 0, laBt sich das Matrizen-Iterationsverfahren anwenden. Wahrend 
man fiir ~« = — oo auf diese Weise auch den tiefsten EKigenwert gewinnen kann, liefert dasselbe 
Verfahren fiir « = -++ oo nur die héheren Eigenwerte. Das liegt daran, daB die Integralgleichung 
und die zugehérige Differentialgleichung fiir den ersten Eigenwert entartet. Wie unten ausgefithrt 
wird, ist im Falle ~ = + oo der tiefste Eigenwert A gleich Null, und das unbestimmte Produkt 


Aa = 0: co nimmt den Wert 1 an. Fir die Falle x = — oo, — 0,5182, 0 und + 0,4254 sind die 
aus 4-reihigen Matrizen ermittelten Naherungen fiir den Drillwinkelverlauf in Abb.6 dargestellt. 


b) Die strengen Kippwerte und die zugehérigen Verformungen. Um die Scharfe 
der oben erwahnten, naherungsweise berechneten Eigenwerte zu erkennen, wird nun die strenge 
Kippgleichung hergeleitet und ausgewertet. Mit dem Eigenwert 


eed IBY 
h —) = BD 


und der Koordinate z = s ist die Kipp-Differentialgleichung durch (29) und ihre Lésung durch 
(31) bzw. (31a) gegeben. Die zugehérigen Randbedingungen sind in (10) verzeichnet. Diejenige 
an der Stelle z = 0 liefert unter Beachtung von (31a) 


(3/4) 


fa 2 « 15/4) 


SRG (35) 


und diejenige an der Stelle z = 1 unter Beachtung von (31) 


Be 5 
C, Jia (7) oy Ge J 44 2 a0 (36) 


Die Forderung, daB die Nennerdeterminante des durch (35) und (36) gebildeten, homogenen und 
in C, und C, linearen Gleichungssystems fiir nichttriviale, von Null verschiedene Lésungen ver- 
schwinden muB, fiihrt auf die strenge Kippgleichung 
es 
1 32 LG/4) Va tN 2, 
[2 


1 _ gan POI) /Va NT (37) 
E T/A) V 2 Al 
all 
Daraus folgt mit /'(3/4) = 1,225417, [°(5/4) = 0,906403 und der Abkiirzung yal2 = 
1 & Sin) ( 
— = 2,092099 ya =. 37a 
a ibs J_144(%) ( 


Die aus (37) berechneten, zu gewahlten yA-Werten gehdrigen Werte « findet man in Tabelle 2 
(erster, zweiter und dritter Eigenwert). Fiir einen Wert )/A, welcher mit dem ihm entsprechenden 
Wert « durch (37) verbunden ist, ergibt sich aus (36) das Verhaltnis der Integrationskonstanten. 
Der Drillwinkelverlauf lautet dann gemaB (31) mit der nicht bestimmbaren Konstanten C, 


fA g 
Hz) = CG yz jus) i (7) Foil? ) (38) 


Fiir ~ = 0 vereinfacht sich (37) zu der schon von Prandtl berechneten Gleichung 


i B —0, | (37b) 
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ep, /B ; 
Tabelle 2. Zusammenstellung der in Abhdngigkeit von « = Vz berechneten ersten, zweiten und 


dritten Kippwerte VA = /ED 


| 5 1, (nah.w.) | Fehler exakt A, (exakt 
re | Vales an oe | in % ie Va (e an S Foe an) 
; 7 | 
fore) | 0 foe) 5,9618 foe) | 11,8123 
19,9984 0,05 10,6588 5,6 3,2878 11,9 
9.9963 0,1 3,0628 il 1,6001 12,0 
4,9933 0,2 0,2000 0 1,8400 5,8 0,3283 13,0 
1,9832 0,5 0,5000 0 1,0626 6,0 0,1934 14,0 
0,9658 1,0 1,0000 0 0,3845 7,0 0,1202 15,0 
0,4254 | 2,0 1,9991 0,05 0,2293 8,0 0,0500 | 16,0 
0,2004 | 3,0 2.9909 0,30 0,1335 9,0 0 16,5159 
0,0031 4,0 | 0,0334 | 10,0 —0,0762 | 17,0 
0 4,0126 3,9755 0,93 0 10,2461 —0,3463 17,6 
—0,5182 5,0 4,9101 80 —0,1898 11,0 —0,4690 | 17,7 
—0,9236 5,2 —0,7520 LS —0,6761 | 17,8 
—1],3524 33 | —1,1734 11,6 SSS 3) 1a) 
—2,3043 5,4 —2,3402 EET —2,5406 18,0 
—6,3234 Deo ROS oow! 2,20 —22,4440 11,8 —3,9922 18,03 
nae 5,5618 | 5,4369 9,25 — oo 11,8123 —oo | 18,0848 


mit den Nullstellen ya= 4,0126, 10,2461, 16,5159,... Gleichung (36) liefert C, = 0, und der 
Drillwinkelverlauf wird 


rs 2, 
OG) = One lila 2) (38a) 


Fiir « = + co werden in (37) die beiden Bessel-Funktionen durch unendliche Reihen (Potenz- 
reihen in A) ausgedriickt [siehe (3la)]. Durch Umformung ergibt sich dann aus (37) 


c FAECES Fey 

2 24 7(2) (9/4) | ; 1 

secre TE) ae ie ak Gee 
SET) POY Oe | = 


Gleichung (37c) ist erstens erfiillt, wenn nach dem Grenziibergang ~ —> + co die Klammer des 
Zahlers verschwindet. Das ist fiir solche ya der Fall, die der Gleichung 


Jiyl(Z) =o 
geniigen, wobei deren erste Nullstelle ya = 0 ausgeschlossen ist 
yA =5,5618, 11,8123, 18,0848,... 
Gleichung (36) liefert C, = 0, und der Drillwinkelverlauf wird 
O(2) = C, V2 Jys ie 2) : (38b) 
Gleichung (37c) ist jedoch zweitens auch fir ya = 0 und yao = 1 erfillt. Das ist sofort einzu- 
sehen, wenn vor dem Grenziibergang (37c) mit « multipliziert, 2 in den eckigen Klammern Null 


gesetzt und nach dem Produkt Aa aufgelést wird. Um bei der Errechnung des Verhaltnisses 
C,/C, der Bestimmung von Grenzwerten unbestimmter Ausdriicke zu entgehen, werden die neuen 


Konstanten C, und C, so gewahlt, daB in den eckigen Klammern von (31a) nur ganzzahlige Potenzen 
von / auftreten 


‘ 1/8 7 = 
Gi = OE CG un C, = 4 ae — 

V2 T(3/4) AMS 1713/4) 
Mit 2 = 0 aus (37c), sowie C, und C, lautet der Drillwinkelverlauf nach (31a) 
O=C24+G. 


Q° 
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Die Randbedingungen (10) liefern G =— GC und ra = 1, somit 
Hg (12) - (38c) 


Der tiefste Eigenwert ya = 0 fiir ~ = + oo besagt, daB die Kippstabilitat gleich Null ist, daB also 
der Trager beim Aufbringen der geringsten Last sofort wegkippt. Diesem Falle kommt als Grenz- 
wert und als Ubergang des Kipp-Problems zu einem Spannungsproblem gewisse Bedeutung zu. 
Der wegen ya a = 0+ co = 1 unbestimmte, aber endliche Ausdruck 


ya oe D(e=0) — GC 


bedeutet eine mechanische GréBe, die einem Torsionsmoment proportional ist, welches am freien 
Tragerende wirkt und den Stab verdreht. Dieses Torsionsmoment bedingt nach der de St. Venant- 
schen Torsionstheorie einen linearen Drillwinkelverlauf, dem (38c) genau entspricht. 


Va 2.1? 
16,0848 A V 


16,5759 


15 


Verlaut des 3. Eigenwertes 
BOG Sie Eee Aa me —— 
70,2467 
gl 
5 Verlauf des 2. Figenwertes 
99618 Bee eee = 
== \ is a BS 4569 
=A \ 
\ 
M 40726 
5975S 
iN 
\ 
Verlaut des 7. Figenwerres _&_/B 
rae. 
eee 
6 -Y -2 0 % 4 6 


Abb. 5. Beispiel 2, Kippwertverlauf. 
Kippwertverlauf gema (37); 


By iam » (13) und (21); 
—— ¢——"s — aa ay (Gk) 

Die Berechnung der Higenwerte fiir groBe positive a zeigt, da dem Werte « = —oo der erste 
Eigenwert ya = 5,5618 und « = + oo der erste Eigenwert yA = (0) zugeordnet ist (Abb. 5). In 
diesem Bild ist auch die von Prandtl berechnete Naherungsformel fiir die Kippkraft bei kleinen 
Exzentrizitaten eingetragen. Sie lautet mit unseren Bezeichnungen 


VA = 4,0126 — 4,1334 0 (39) 

und liefert im Bereich —0,2 < « < + 0,6 hinreichend gute Naherungswerte. . 
Abb. 6 beinhaltet fiir die oben angegebenen Werte « den exakt und naherungsweise be- 
rechneten Drillwinkelverlauf. In Abb. 7 ist die w-Linie aufgetragen. Sie gilt fiir alle VA, auBer 


2 


18 Kraus: Die Integralgleichungen der Kippung gerader Trager Ingenieur-Archiv 


fiir den ersten Kigenwert ya —0 bei « = + © und ergibt sich aus (33) mit M(x) = Px und den 
Randbedingungen w = w, = 0 fiir x = 1 zu 
heey 2 


PP 
B! 6 E J max 
Abb. 7 enthalt weiterhin die mit ¥-Gewichten naherungsweise nach (34) ermittelten u-Linien 
(dieser Berechnung lag nach 0. Mohr ein Ersatzkragtrager mit der Einspannstelle am freien Ende 


des Kipptragerszugrunde). Fiir das Wertepaar « = ~- © und yr —0 tritt keine u-Verschiebung auf. 


Ww 


V/s) 


YA=4,0 


oc = 40,4254 


\A=4,0126 
o=0 


VA'=5,5678 und 0¢ =— 20 
VA'=5,0undoc =-0,5182 
VA= 40126 und oc=0 
VA'=20 und cc = +0,4254 
VA=0 und 0 = +00 


VA=5,5678 


@=-90 


0 OZ G4 G6 08 40 0 02 O4 06 G8 40 
Abb. 6. Beispiel 2, Drillwinkelverlauf, naherungsweise aus 4-reihigen Abb. 7. ‘Beispiel 2, Verlauf von w exakt gemaB (40) und von 
Matrizen und exakt. u naherungsweise mit ¥3-Gewichten gemaB (34). 
--- # aus 4-reihiger Matrix gema4B (13) und (21); ® exakt gemaf (38). 


9. Zusammenfassung. Es werden die Integralgleichungen des Kipp-Problems des geraden 
Tragers mit diinnwandigem, doppeltsymmetrischem und offenem Profil aufgestellt. Man verwendet 
hierfiir die Methode der Variation der Konstanten von Lagrange. Einfache Falle, bei denen die 
exakten Lésungen teils vorliegen, teils berechnet wurden, zeigen die Giite der Naherungen, welche 
die Erhard-Schmidische Ungleichung und 4- bzw. 9-reihige Matrizen liefern. Eine scharfere und 
allgemeinere Untersuchung, sowie die Berechnung der noch nicht vorliegenden, praktisch inter- 
essierenden Falle erfolgt demnachst. Hierbei werden auch die Integralgleichungen des Kipp- 
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tragers mit einfach-symmetrischem Profil mitgeteilt und insbesondere der Kragtrager behandelt. 
Weiterhin ist die Bestimmung der Biegedrillknick- und Kippstabilitat von Darchlauftragern 
mit doppeltsymmetrischen Profilen Gegenstand einer weiteren Arbeit 1. 

Es ist dem Verfasser eine angenehme Pflicht, seinem verehrten Lehrer, Herrn Prof. D 


K. Karas ferner Herrn Dr. W. Bérsch-Supan vom Institut fiir Praktische Math 
stutzung zu danken. 


r.-Ing. 


ematik fiir ihre Unter- 


Mitteilung aus dem Institut fiir Angewandte Mechanik und Technische Schwingungslehre 
der Technischen Hochschule Darmstadt (Prof. Dr.-Ing. K. Karas) 


(Eingegangen am 23. Februar 1957) 
Anschrift des Verfassers: L. Kraus, Wiesbaden, Idsteinerstr. 11 


i lita a lir ei a Bull. Nr. 10 
1 Die Kippstabilitat von Durchlauftragern wurde fiir einfache Sonderfalle von O. Pettersson, 
from The Div. of Building Statics and Struct. Eng. at the Royal Inst. of Techn., Stockholm, behandelt. a 
neueste Arbeit auf dem Gebiet des Kippens ist die Abhandlung von H. Nylander, Tekniska Hogskolans Hand- 
Iingar Nr. 102, zu benennen. Auf beide Veroffentlichungen wurde der Verfasser in freundlicher Weise von 
Herrn Prof. Dr. E. Chwalla, Graz, aufmerksam gemacht. 
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A Study on the Period of the Free Lateral Vibration of the Beam Bridge 
by the Theory of the Orthotropic Rectangular Plate | 


By M. Naruoka and H. Yonezawa 


1. Introduction. Important data necessary to discuss the dynamic rigidity of a bridge can 
be obtained by researches on the free lateral vibration of the bridge. The free lateral vibration of 
the beam bridge has been researched generally by the theory of the free lateral vibration of the 
beam and the co-operation between the beam and slab has been neglected in the calculation of the 
period of the free lateral vibration. Therefore, there are sometimes cases when the calculated period 
differs considerably from the measured value. As a more rational method of calculating the period 
of the free lateral vibration of beam bridges, the theory of the continuous isotropic plate supported 
by flexible beams or of the grillage beam can be applied to the analysis of the free lateral vibration 
of the beam bridge, but if the number of bridge beams increases or if the torsional rigidity of the 
beam and slab must be taken into consideration, the calculation by these theories is very difficult. 
On the other hand, only the theory of the orthotropic rectangular plate makes it possible to apply 
the theory of thevibration of the plate to the beam bridge without regard to the number of | 
beams. 


In this paper, the theory of the orthotropic rectangular plate is theoretically and experimentally 
applied to the analysis of the free lateral vibration of the right beam bridge, and the result of the’ 
analysis is compared with that by the conventional analysis. Thus, it was made obvious that the | 
theory of the orthotropic rectangular plate is very effective for the analysis of the free lateral 
vibration of the beam bridge. | 


2. Frequency equation of the orthotropic rectangular plate with two opposite edges simply 
supported and the other edges free. The differential equation for the free lateral vibration of the 
orthotropic plate is represented generally by the following well known equation by M. T. Huber: 

ov. aw ow ew 
Depa Te papa ee ee ee 


with (1) 
2H = 7, Da vy Da ee 


where, @ is the mass of the material of the plate per unit area. 
D,, and Dy, are the flexural rigidities in x- and y-directions. 
Equation (1) can be written in the following form: 


simply supported 


ow [Dy BW | Dy ow 0 &@Ww 
ri 2 Dee eee 
Ox’ DNOx7 ION 4 D,, dy* De set | 
in which | 

1g! 
9 Ss 4 
| Dz Dy 

simply supported When a plate performs a normal mode vibration, the 


deflection surface can be represented as follows: 
W =w (A cos pt + B sin pt) 


Fig. 1. Orthotropic Plate. 


where w is the normal function and p is the circular frequency. Substituting the above equation 
in eq. (2), we obtain 


4 d 2 
cay 2 m/s Ow , Dy dw gp? 


axt oo— 05 


D,, 0x2 dy? ' Dy dy! D, (3)) 


er eau solution of equation (3) varies according to the condition x > 1, x =1, x < 1} 
and 70; | 


We shall consider the particular case in which the boundary conditions are that the two edge| 
x = 0 and x =a of the plate are simply supported and the other edges are free at y = + b/2 as 
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~shown in Fig. 1. Then we have the following end conditions: 


i=. 
tees on aap for the edges x = 0 andx =a, 
x x\ 9x2 | Vy dy? 
02w oh 
ie D, (Fs i ra) =o for the ed b/ 
A onde = ’ for the edges y = + 0/2. 
i= Ow ! Vp (2% P EE ee = 
dy? Dy, y Dy Oy Ox? 
In these boundary conditions, we shall assume as Vv, = Vy = v to make it simple. 


For the sake of convenience, the frequency equation for only the two cases x = 1 and x = 0 
will be induced. The result of calculation is as follows: 


if 1/2 
a) Case reel OW | ae y —l1>0, 
INDY Nine 
2\ / a \4]1/2 
i ee ae ea v2 = |(2F) (na) ae 


Ey (2 | a ) | ale (2 5) cosh = sin a = 
[Dz é Ds. a0: 
iv (2 \/ Dy as ») ye — (2 =—=— 35) sinh = Peon = Tees =) (4) 


for the vibration symmetrical to x-axis, 


yx ( \z af br ae (2 —yY ie “| sinh 5 hoy Fa med 
—y (ya |/ Ps = 7 as cosh" sin =? = 0 (5) 


for the vibration reversely symmetrical to x-axis; 


2 1/2 
b) case | and (eZ 4 | SW), 
D; /\mx 


lf / de : £ ie 
ae ae a a Ayr =1+(E)(2)| 


yy! (2 \z —y bv (2 Y V5") cosh i sinh a = 


ENT O46) 
—7(y2/3e—») he —(2 —»|/ 3) sinh "cosh", = (6) 
for the vibration symmetrical to x-axis, 
+ DAV b } 
2) ¥ (7 VE) by (2 VF) sinh “> cosh “> — 
—y[y2 22-9) y? (2. VF 2) cosh sinh,” = = (7) 
\ y¥ 
for the vibration reversely symmetrical to x-axis 
2 4 
c) case rea) and y= (2F)(2)_1>0, Tey 
Dd; Dy na ga 
1) (2 Vo- »] be (2 v | cosh “5 sin 5 +- 
rb (3) 


VE Hey Bn 
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for the vibration symmetrical to x-axis, 


2) (2 V3 —y Le (2 —y V3) sinh fs cos Bs — 


Dy D: pli ce ele 
= (2 V3 + ») be (2 y V2) cosh =a sin ee 0 
for the vibration reversely symmetrical to x-axis; 
= os op*\/ a \4 _%y 
d) case ~*=0 and yt=1 (Be) ea) =0., . jz 


1) be (,2 V5 =r ) = (2 V3 y )(2 y V3) sin 3 cos = SE 


De Ds DAW Merb rb 
+ [2 (AV BE») + (AY Be + 9) (2—v VB] sin 7S cosh 2 — o 


for the vibration symmetrical to x-axis, 


TD D; D,\| . rb rb 
2 foley) (Veo Bette 
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(9) 


(10) 


Ds D, 1/Dz\|_. 4 3rb rb 
— bel Vax) + (V5 +9) rf Fe) nnn oooh P= 0 


for the vibration reversely symmetrical to x-axis, 


where 6 = —|/ =. 
a Dy, 


For the particular case when the plate is supported by flexible beams at both edges y = + 6/2, 
the frequency equation for the vibration symmetrical to x-axis becomes as follows, under the | 


assumption vy = 0: 


2 4|1/2 
ess lle eral ae —1>0, 


2 4)1/2 
r= oy!s yt = (2)-2)) wy, 


a) Case 


x 


; b “b ; ‘ eee Ts 
@ (y? + y'?) cosh 5 cos - Se (, sinh i cos = + y' cosh = sin _ a 
‘ek i ribo Tb Fas ele ead Po TDN ORs. 
2yy (7 cosh 7 sin =- — sinh 3 Cos 7 ns 


b) case 


4j1/2 
vor ot [leo 


41/2 
T=0Y, r=ay': aye =1+(()(4)] ; 


Mm I 


® (y? — y"?) cosh a cosh a + p22 (, sinh x cosh se —y’' cosh S sink ) aa 


2 
+2yy' (r cosh i. sin i: =! sinh fe cos 7) ==) ()R 
c) case == 00 Sand = (2 cae 1 S0eineaigwys 
2D — 4 (tanh = - tan 2 = 0% 
d) case x=0 and Miele! oP) 2 Un a yl 2 > 


: b A b b 3 
sin? 5 mes (sinh > cosh > -++ sin fs cos = = 0 


/2@ (cosh 


(6') 
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ee 
where 1 eel yee ees 
Xn a BD, > D — De —f X 
and E Jp = flexural rigidity of edge beam. 


3. A theoretical study on the period of the free lateral vibration of the beam bridge. Let us 
consider for instance the beam bridge as shown in Fig. 2. As the co-ordinate axis, the direction of 
the main beam and the direction of the floor beam or of the short span of the slab are assumed as 
x-axis and y-axis respectively. 

a) On the deflection surface. We shall consider the beam bridge of the T-section. When 
the fundamental frequency of the beam bridge is calculated by the beam theory, the vibration 


IZ 


Fig. 3. Fundamental Mode of Vibration of Simple Beam. 


b 2 | 


= + 


2 2 


Wig. 4. Vibration Mode Having one Node in y-Direction. 


ae 


Fig. 2. Beam Bridge Considered as Orthotropic Plate. Fig. 5. Vibration Mode Having Two Nodes in y-Direction. 


without nodes as shown in Fig. 3 is most important and the frequency and the corresponding period 
of this mode of vibration are obtained by the well known equation 
Che E ih re 
Pies 22 f= a ee — + = a pes (12) 
where I, is the moment of inertia of T-section of which the effective width of the compression 
flange is the distance from center to center of the slab, and g, the mass of the above T-section per 
unit length. 

Accordingly, when the span and the section of the beam are equal, the frequency becomes 
equal, even if the number of main beams or the width of the beam bridge or the relative stiffness 
of the beam to that of the slab are different. This is obviously irrational when we consider that 
the main beam, floor beam, stringer and slab resist the load as co-operating members and operate 
together, The theory of the orthotropic plate gives rational results to these problems as can be 
unterstood from the author’s study on the experimental stress analysis on the steel highway 
beam bridge?. 

In the analysis by the theory of the orthotropic plate, the deflection surface in x-direction is the 
same as with the case of the beam, but the deflection surface in y-direction should be discussed. 


1 W. Cornelius, Stahlbau, 21 (1952), S. 21, 43 und 60. 
2 M. Naruoka, Publication of I[ABSE, 14 (1954), p. 183; M. Naruoka and H. Yonezawa, Proceedings of the 
5th Congress of I[ABSE (Lisbon), Preliminary Publication 1956, p. 393. 
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~ 


With the beam bridge, of which the span is large compared with the width of the bridge and which 
vibrates analogously to the beam in the whole, the vibration which has no nodes in y-direction is 
important and equation (12) can be applied to this case. However, with the beam bridge having 
a width nearly equal to or larger than the span, it is supposed that the vibration having one or two 
nodes in y-direction as shown in Fig. 4 


oleae and 5 occurs. Which form of vibration 
7 occurs, the symmetrical or the reversely 
symmetrical to x-axis, must be solved 
@ dan by experimental researches, but we shall 
i sie t compare both forms by theoretical con- 
F ! a/b=z lh sideration. 
Seon tar ——- The spare apee Se of rs 
fe = (p @/x?) (9/D,)"! 
2 1 which satisfy equation (4) and (5) for 
; beam theor the various values of D,/D, = 4 ~ 400, 
Dalby under the asumption of a/b = 2,% = 1 
0 700 200 Ba aes Wo and y = 0 and the values of y calcula- 
Fig. 6. Circular Frequency for Vibration Mode (Fig. 5). ted by the beam theory are shown in 


Fig. 6 and 7. It is a remarkable charac- 
teristic that the frequency of the vibra- 
tion having nodes in y-direction is con- 
siderably large in comparison with that 
obtained by the beam theory even if the 
vibrating form in x-direction is the same. 
Furthermore, the value of mw changes 
beam theory considerably with the values of relative 
stiffness D,/D,. 
Da/Oy b) On the effect of the value 
0 100 200 500 400 of x. The value of x in M. T. Huber’s 
Fig. 7. Circular Frequency for Vibration Mode (Fig. 4). differential equation can only be deter- 
mined experimentally and differs accord- 

_ pat, [9 ing to the structure of the beam bridge. 
Here, we shall discuss the effect of the 
value of x on the frequency. 

The values of uw which satisfy equa- 
tion (4) and (8) or (5) and (9) under the 
assumption x = 1 or x = 0 are shown 
in Fig. 8 and 9 for the values of a/b = 2, 
y =0 and D,/D, = 4 ~ 400. The fre- 
quency under the assumption of x = 1 
is large in comparison with the case of 

| bifly %=0, which is quite natural, considering 
0 700 200 300 400 that the torsional rigidity is neglected in 
Fig. 8, Circular Frequency for Vibration Mode (Fig. 5). the case of x = 0. Itis supposed fromthe 
results of our field test that the calcu- 
lation under the assumption of x = | is 
applicable to the reinforced or pre- 
stressed concrete T-beam bridge and 
that of x = 0 is applicable to the com- 
posite grillage beam bridge. But, the 


beam theory 


| c 
0 700 200 > 300 700 frequency shows no remarkable diffe- 
Fig. 9. Circular Frequency for Vibration Mode (Fig. 4). rence whether the assumption 1S. eh 
or = 0. 


c) On the effect of the values of the ratio of span to width a/b and the Pois- 
son’s ratio y. To investigate the effect of the values of a/b and v on the frequency, the frequency 
will be calculated for the cases of a/b = 1, 2,4, 8 and »y = 0, vy = 0,15 under the assumption of 
x = 1 and D,/D, = 4 ~ 400. The results are shown in Fig. 10, a ~d. It is natural that the 
frequency (vy = 0) is larger than the frequency (vy = 0.15). 


f 
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pa’, [o 
ff a Vor 


beam theory, 


DafDy 


0 100 200 300 400 
Fig. 10 a). Circular Frequency for Vibration Mode (Fig. 5). 


Fig. 11. Model Beam Bridge (G, ~ G;). 
b 


beam theory 


—— 
0 700 200 300 400 


Fig. 10 b). Circular Frequency for Vibration Mode (Fig. 5). 


A 
Fig. 12. Model Beam Bridge (G, ~ G;). 


Fig. 13. Dimension of Model Grillage Beam Bridge 


Fig. 10 d). Cireular Frequency for Vibration Mode (Fig. 5). G, (unit: cm). 
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Next, let us consider the relation between the fundamental and higher order frequencies. The 
ratio 1:4:9 is found in the frequency of the fundamental, 2nd and 3rd vibrations of the simple 
beam, but, this ratio becomes ramarkably different in the case of the orthotropic plate according 
to the values of a/b and D,/D,. This ratio for the case of the orthotropic rectangular plate is given 


in Table 1 for the various values of a/b 
Table 1. Ratio of Frequency of Ist, 2nd and 3rd Order Vibration dD./D.. 
of Rectangular Orthotropic Plate with Each Two Opposite Edges and D,} y 


Simply Supported and Free It is presumed from this table that 


the interval of the ratio approaches to 

x/Dy ; i Pn each other when the ratio b/a and the 
afb | relative stiffness are large, and, on the 
other hand the ratio of the frequency 


1 petletss:2:6 1:1.1:1.4 ig NE} becomes larger than that of the simple 
- | Tonia ace eee beam 1:4:9 when the ratio b/a becomes 
1:7.1:73 1:4.0:17 1:34:13 small. This suggests that the vibrations 


having nodes in y-direction generate 
easily in the beam bridge of which the width is large compared with the span. 


4, Experimental study on the model beam bridge. The results of the measurement of the period 
of the free transverse vibration of the model beam bridge shall be compared with the results of the 
analysis by the theory of the orthotropic plate or by the beam theory in the following. 


Table 2. Dimension of Model Beam Bridge 


a (cm) 40 40 40 
b (cm) 20 20 20 
2 (cm) 4 4 2.9 
e (cm) il 1 1 
h (cm) 0.5 0.5 0.5 
number of beams 5 3 8 
Table 3. Dimension of Model a) Model beam bridge. The model beam bridges are parallel 
Beam Bridge (unit: cm) beam bridges G,, G,,G3,G, and G,, grillage beam bridge G, and 
G, | Gs ladder type beam bridge G,. The models G, ~ G, are made of cast 


iron, Young’s mudulas of which is 1.03 x 10° kg/cm?, and G, und G, 


= si are made of steel. The dimensions of these models are shown in 
4 9 Table 2 and 3 and in Fig. 11, 12, 13 and 14. The periods of these 


z models are measured by cathode-ray oscilloscope with single sweep 
5) circuit and cristal pickup. 

b) Results of measurement and comparison with theo- 
retical results. The measured or calculated values of the period 
of the free lateral vibration of these models are shown in Table 4 and 5. The measured periods on 
G, ~ G,; coincide with the period calculated assuming the beam bridge as an orthotropic plate 
vibrating with two nodes, G; and G, coincide with the period as an orthotropic plate vibrating 
reversely symetrically with one node and G, coincides with the period of a simple beam. It is 
well expected that the period of G, will coincide with that of the simple beam, considering that 
the ratio of span to width is 5. 

We studied the nodal line on the vibrating plate by the fine and dry sand scattering on the 
plate (b/a = 1/2), and as the result of this experiment, it could be made clear that there are two 
nodal lines perpendicular to y-axis. 


omanove 


5. An experimental study on existing beam bridges. It was proved by the author’s experimen- 
tal study on many models that the theoretical results above mentioned are proper but in this 
chapter the experimental studies on the existing beam bridges will be described. 

a) The Kugenuma Bridge (reinforced concrete T-beam bridge). This bridge is 
shown in Fig. 15. Assuming FE, = 2.1105 kg/cm, x = 1, » = 0 or y = 0.15, the average height 
of beam = 62 cm and D,/D, = 56 and neglecting floor beam, the period of the free lateral vibration 
symmetrical to x-axis is calculated by eq. (4). 
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Table 4. Calculated and Measured Periods of Tested Mode Bridges (unit: sec) 


i ee 
ey Theeretical Value =i 
\ | Oribewepier Plate a 
ed | ia pep er on ee | Pye 
el ee ee 
\ | as | 
(a) | 0.00123 0.00160 0.00101 0.00142 | 0.00121 
G, (b) 32 | 222 | 115 | 175 | mary 
(c) $45.94 369 135 261 | 157 
(d) 147 | 1124 155 412 | 169 
(a) | 0.00115 0.00160 0.00100 0.00140 | 9.00120. 
G, \(b) 125 | 224 123 186 | 146 
(c) 130 | 424 139 263 162 
(d) PZ, 9) 1087 | 162 | 432 | 176 
(a) | 0.00093 0.00164 0.00106 0.00149 | 0.00127 
G, (db) 96 | 225, 123 181 | 149 
(c) 1G | 364 145 260 170 
(d) 135 | 1034 | 162 | 432 177 


Y ¥ 


Fig. 14. Dimension of Model Ladder Type Beam Bridge G, (unit: em). Fig. 15. Dimension of the Kugenuma Bridge (unit: em). 


Table 5. Calculated and Measured Periods of Tested Model Bridges (unit: sec) 


\ Theoretical Value 


Orthotropic Plate 


Measured | Dy Dy H? — D,, Dy H—=0 
Value | Beam a —= ————— = thf atte oS a 
i 5 
a | gaan ne < S 

Gy 0.0046 0.0050 0.0016 0.0031 

G; 0.0036 0.0113 | 0.0024 

G, 0.00063 0.0011 0.00017 | 0.00045 | 0.00018 

G, 0.0040 0.0040 | 0.0013 
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The periods measured or calculated are given in Table 6 with the values by the theory of 
the simple beam. The value calculated by the theory of the orthotropic plate can explain well] 


the measured value. 


Table 6. Calculated and Measured Period of the Kugenuma Bridge (unit: sec) 


| | Orthotropie Plate, H? = Dy Dy 


Measured Beam 10) y = 0.15 
Value 
0.040 ~ 0.065 | 0.103 ~ 0.106 — 0.055 | 0.059 


b) The Sagoshi Bridge (composite grillage beam bridge). This bridge consists of 
four main beams and one cross beam and the skelton is shown in Fig. 16. Assuming E, = 2.1 x 10° 
kg/cm?, E,/E, = 10 and y = 0, the period is calculated by eq. (4), (4’), (8), (8’) and (12). The measu- 
red calculated values are shown in Table 7. 


+ | 
Me a 
S Ri 
S| J=487x 10%cm* 1S) 
2 & 
o 1° 
ME =§57x10°cm* Hi 
i 75m ~75™ —-| 


Fig. 16. Arrangement of Main Beams and Cross Beam of the Sagoshi Bridge. 


c) The Shohei Bridge (prestressed concrete beam bridge). This bridge is composed 
of fifteen main beams and three cross beams, and the span and the width are 12 m and 11] m re- 
spectively. Assuming D,/D, = 9, E, = 4.0 x 10° kg/cm?, » = 0 and x = 1, the period is calculated 
by eq, (4), (5) and (12). The calculated and measured periods are given in Table8. The value 
calculated by the theory of the simple beam does not agree with the measured value, but, on the 
other hand, the value calculated by the theory of the orthotropic plate for the symmetrical vibrating 
orm with two nodes in x-direction can explain well the experimental result. 


J=7735x10%em* 


1 @ 14=7198m 


444 x10%em* 


Te 


A=905 x 10%cm* 


L 4(@50=2h0m Z| 


Fig. 17. Arrangement of Main Beams and Cross Beams of the Shigita Bridge. 


d) Shigita Bridge (composite grillage beam bridge). This bridge is a composite 
grillage beam bridge, the span and width of which are 20 m and 22 m respectively, and the numbers 
of main beams and floor beams are 12 and 3 as shown in Fig. 17. The period is calculated by equa- 
tions (4), (4’), (8), (8) and (12) and the calculated and measured periods are given in Table 7. 

It is obvious from these tables that the periods calculated by the theory of the orthotropic 
plate coincide with the measured values, but the values calculated by the beam theory do not. 
With the beam bridge of which the width is almost equal to the span such as Kugenuma, Shohei 
and Sigita Bridges, the calculated period for the symmetrical vibration having two nodes in x-di- 
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rection coincide well with the measured period. With the beam bridge of which the span is rela- 
tively large compared to the width such as the Sagoshi bridge, the calculated frequency calcu- 
lated by beam theory coincides with the measured frequency. 


Table 7. Calculated and Measured Period of the Sagoshi Bridge and Shigita 
Bridge (unit: sec) 
a ae a EE ee 
| 


Calculated Value 


= Remarks 
by orthotropic plate 


0.127 (A, C, x = 0) 


DED, An 

The Sagoshi 0.112 (A, C, x = 1) 

Bridge 0.136 | 0.147 

g 0.134 (B, C) x = 0) i 

0.116 (B, C, x = 1) hae Risa 
0.123 (4, C, x = 0) rare 
0.105 (A, C, x = 1) alyD = 6. 

The Shigita 0.153) Ce 0) 

: 0.127 , 0.189 | D,/Dy=87.5 
Bridge 0.132 (B, C, x = 1) 5 =|Dy 


0.150 (A, D, x = 0) | 
0.112 (4, D, x =1) 


D,[Dy = 6.4 


i 


Remarks: A means that the formula {D, = E J/a, Dy =n E J/l (x: number of cross beam)} 


is used in the calculation of the flexural rigidity of the orthotropic plate. — B means that the for- 
mula {D, = E J/a, Dy = Eh'/12 (h: thickness of slab, »y = 0)} is used in the calculation of the 


flexural rigidity of the orthotropic plate. — C means that the orthotropic plate is supported by the 
edge beams at two edges. — D means that the orthotropic plate is free at two edges. 


Table 8. Calculated and Measured Period of the Shohei Bridge (unit: sec) 


4 + aa Pelt) ree S D,, Dy : 
Measured Beam | | 
Value | 
__ eee ee 
0.063 0.185 0.071 | 0.123 


6. Conclusion. The frequency of the right, simple beam bridges is discussed by the labo- 
ratory study on the model beam bridges and the field tests on existing beam bridges in this 
paper, but furthermore, the application to the skew beam bridge or to the continuous beam 
bridge and the problem of forced vibration should be studied. 

Even if more exact studies are necessary, it is made sure by the study in this paper that the 
theory of the orthotropic plate is proper to the analysis of the free lateral vibration of the beam 
bridge and sufficient enough for practical purpose, especially for the beam bridge of which the width 
is almost equal to the span. 


(Eingegangen am 4. Marz 1957.) 
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Professor Masao Naruoka, Kyoto (Japan), Kyoto University, Department of Civil Engineering. 
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Praktische Methoden zur Lésung von Problemen der stationéren 
und instationéren Grundwasserstromungen 


Von G. Heinrich und K. Desoyer 


1. Einleitung. In zwei vorangegangenen Arbeiten! wurde die Theorie der Grundwasserstrémun- 
gen einerseits auf die Behandlung von instationaren Fliissigkeitsbewegungen, andererseits auf den 
Fall anisotroper inhomogener Festkérper erweitert. Im folgenden soll die Anwendung der ent- 
wickelten Gleichungen auf praktische Aufgaben der Grundwasserstrémungen gezeigt werden. 

Um zu einer definierten Lésung eines vorliegenden Problems der Grundwasserbewegung au 
gelangen, ist eine genaue Kenntnis der Rand- und Anfangsbedingungen sowie eine quantitative 
Erfassung der Durchlassigkeitseigenschaften des Festkérpers erforderlich. Auf die praktische Er- 
mittlung dieser Daten soll hier nicht eingegangen werden”. Die in Frage kommenden Rand- und 
Anfangsbedingungen wurden bereits in der ersten der vorangegangenen Arbeiten besprochen. 
Fiir die eindeutige Lésung einer praktischen Aufgabe bendtigt man ferner die Form und Lage der 
undurchlassigen Schichten sowie der Grenzflachen des Grundwassertragers gegen Luft und gegen 
ruhende Fliissigkeit. 

Wenn ein Problem ohne freien Grundwasserspiegel vorliegt, dann handelt es sich bei homogenen, 
isotropen Grundwassertragern im stationaren Fall um ein Potentialproblem mit festen vorgegebenen 
Randern und vorgegebenen Randbedingungen; im instationaren Fall kénnen die Randbedingungen 
vorgegebene oder aus der Problemstellung ermittelbare Funktionen der Zeit sein. Probleme dieser 
Art bieten meist keine besonderen Schwierigkeiten und sollen daher hier nicht weiter in Betracht 
gezogen werden. Sie kénnen aber auch nach der im folgenden genauer besprochenen Relaxations- 
methode behandelt werden. 

Mathematische Schwierigkeiten besonderer Art treten bei Vorhandensein von freien Grund- 
wasserspiegeln auf. Der freie Grundwasserspiegel ist als Randflache des Strémungsbereiches mit 
vorgegebenem Druck anzusehen, jedoch ist der Grundwasserspiegel bei stationaren Strémungen, 
bei denen er iiberdies eine Stromlinienflache sein mu’, nicht vorgegeben, sondern muf bei der 
Lésung des Problems mit bestimmt werden. Zur strengen Lésung von Aufgaben dieser Art wird 
vornehmlich die konforme Abbildung in Verbindung mit der Hodographenmethode verwendet’. 
Als Naherungslésung ist im hydraulischen Schrifttum der Dupuitsche Ansatz gebrauchlich, der 
jedoch, wie im folgenden genauer erértert wird, sowohl in grundsatzlicher wie in praktischer Hin- 
sicht sehr bedenklich erscheint. 

Bei der theoretischen Behandlung von instationéren Grundwasserstr6mungen mit freien 
Spiegelflachen wurde bisher fast ausschlieBlich die Dupuitsche Naherung verwendet’, die in diesem 
Fall ebenso bedenklich erscheinen muf8. 


2. Die Dupuitsche Niaherung als strenge Lésung fiir ein anisotropes Festkérpermodell. Im 
folgenden soll zunachst gezeigt werden, daB nur fiir ein homogenes, anisotropes Festkérpermodell, 
bei dem die Durchlassigkeit in vertikaler Richtung unendlich groB ist (,,einachsig widerstandsfreier 
Festkérper“), die Grundwasserstrémung dem Dupuitschen Ansatz streng gehorcht. 

Nach Dupuit wird bekanntlich naherungsweise angenommen, daf die Standrohrspiegelhdhe h 
fiir alle Punkte des Grundwassertragers, die auf einer Vertikalen liegen, gleich groB ist®. Die aus 
dieser Naherung hervorgehende Gleichung fiir die Spiegelbewegung soll nun fiir allgemeinere Vor- 
aussetzungen als bisher (beliebige Form der undurchlassigen Schicht, anisotroper Festkérper) 


 G. Heinrich und K. Desoyer, Ing.-Arch. 23 (1955) S. 73; Ing.-Arch. 24 (1956) S. 81. Bei Verweis auf 
Gleichungen dieser Arbeiten wird den betreffenden Gleichungsnummern der Hinweis I baw. IL vorangestellt. 

 Siehe z. B. J. Kozeny, Hydraulik, S. 389. Wien 1953. 

2 Siehe z. B. G. Hamel, Z. angew. Math. Mech. 14 (1934) S. 129; G. Hamel und E. Giinther, Z. angew. 
Math. Mech. 15 (1935) S. 255; V. V. Wedernikow, Z. angew. Math. Mech. 17 (1937) S. 155; G. Hamel, Z. angew. 
Math. Mech. 18 (1938) S. 39; H. F. RoBbach, Ing.-Arch. 12 (1941) S. 221; J. Kozeny, a. a. O. S. 329ff. 

* Vgl. auch E. P. Nemecek, Der Strémungsdruck in Béschungen. Diss. T. H. Wien 1953, in der eine 
andere Behandlungsweise der instationaren Grundwasserstrémung am Beispiel der Durchstroémung eines 
Dammes vorgeschlagen wird. 


° J. Dupuit, Etudes théoriques et pratiques, 2. éd. (1863) S. 254; Ph. Forchheimer, Hydraulik, 3. Aufl. 
S. 70, Leipzig-Berlin 1930. 
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kurz hergeleitet werden. In Abb. 1 sei z(x,y) eine gegebene undurchlassige Flache und z,(x,y,t) 
der freie Grundwasserspiegel. Bezeichnen q, baw. qy die in der Zeiteinheit langs der «- bzw. y-Rich- 
tung durch die Seitenflachen des dargestellten Elementarquaders eintretenden Fliissigkeitsmengen, 
so liefert eine Volumsbilanz 


Cqx 
ox dx ne 


Ody Ok, 
By ie tia dx dy. (1) 
Dabei sei n das hier als ortsunabhangig vorauszusetzende relative Porenvolumen. 


Fiir den Fall eines homogenen, anisotropen Grundwassertragers werde nun vorausgesetzt, dab 
der Durchlassigkeitsaffinor (I, 8) symmetrisch ist 


und die vertikale z-Achse mit einer seiner drei 
Hauptachsenrichtungen zusammenfalle. Als x- und 
y-Richtung seien seine beiden anderen Hauptachsen- 
richtungen gewahlt. Dann ist gemaB (I, 18) 


OPS pale Bas vy = hy Jy, 
wobei die Standrohrspiegelgefalle J, bzw. Jy nach 
der Dupuitschen Naherung von z unabhangig sind 


und durch — @z,/0x bzw. — 0z,/0y ersetzt werden. 
Daraus ergeben sich fiir q, und q, die Ausdriicke 


Oz 
qx ae (z5 oat %,) k,., z= dy > 

OZ (2) 
Gy = — (% — %) hy By Ge: 


Abb. 1. Zur Ableitung der Bewegungsgleichung fiir den Grund- 
wasserspiegel in einem einachsig widerstandsfreien Grundwasser-~ 


Durch Einsetzen von (2) in (1) erhalt man die fol- 


gende Gleichung fiir die Spiegelbewegung: traiger. 
Oz 0 Oz id Oz 
n = == ke ne C Psy Z| L yy By |e z,) “1 : (3) 


Fir den Sonderfall des isotropen Festkérpers (k 


= ky =k) und horizontale undurchlassige 
Schicht (z, = 0) erhalt man daraus die Gleichung 


xx 


(3a) 


die bereits in der hydraulischen Literatur bekannt ist?. 

Wir zeigen nun, daB aus den exakten Gleichungen (I, 16b) und (I, 30b), wenn hierin die ¢-Achse 
mit der vertikalen z-Achse, die ¢- bzw. 7-Achse mit der x- bzw. y-Achse des Koordinatensystems 
der Abb. 1 zusammengelegt werden, beim Grenziibergang k,,—> co die soeben hergeleitete Glei- 
chung (3) folgt. 

Aus (I, 16b) 


Ch Oh Oh 
XX Oy2 hy Oy = k,, 622 =0, (4) 
wo h die Standrohrspiegelhohe ist, folgt fiir k,,—> co die Beziehung 
we : 
E g =A) = SO ae = konst. (5) 
Oz OZ 


Die Komponenten eines Normalenvektors 11 auf die undurchlassige Flache 2, = 2,(x, y) sind 
— 6z,/@x, — @z,/dy, 1. Da die Flache z,(x, y) Stromlinienflache sein muB, gilt die Beziehung 


Oz, py Oz e. 
Dems Oyen Uy Gaal (eae (6) 
Aus dem verallgemeinerten Darcy-Gesetz (I, 12) ergibt sich fiir die Hauptachsenrichtungen 
ch ch ea, ck 
Ui Kenn ine vy > pen kyy éy ? — — k., nee: (7) 


1 Ph. Forchheimer, a. a. O. S. 104. 
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Aus (6) und (7) erhalt man 


oh 02, oh Oz, Chmas 0 8 
Keg dx Ox ' YY Oy ay ee Oz 4 (oF 


giiltig langs der undurchlassigen Schicht z,(x, y). Aus (8) und (5) sieht man, daB die Grikbe dh/ cz 
fiir den Grenzfall k,,-—> co gegen Null geht. Es ist dann 


h = h(x, y) (9) 


von z unabhangig, wie es der Dupuitschen Annahme entspricht. Aus (4) folgt 


Fiihrt man hierin den Grenziibergang k,,—> oo durch, dann hangt wegen (9) der Integrand auf 
der rechten Seite nur yon x und y ab, die Integration ist leicht durchfiihrbar. Man erhalt 


oh * ee 2h 
fascia as es 2) (hes age ths Tok (4a) 


Ox" YY oy* 


Der Wert an der unteren Grenze fiir die linke Seite kann aus (8) entnommen werden. Damit ergibt 


sich 


GLI Gy: ®X Oy Ox JY Oy oy 


3 a CDS lp Ob Oe oh Oz, a y a se Ch 10 
elite, (Fa Bis oe (= Ox Ox ! hyy oy dy (2, 2) (k XX Ox moi Ay oy? . ( ) 
Fiihrt man den Grenzwert (10) in die strenge Gleichung fiir die Spiegelbewegung (I, 30b) 
Oz? un 1 7: ch k Oh 02. gp lp a (11) 
langs der Spiegelflache z = z, 


ein, so ergibt sich nach einfacher Umformung fiir den Spiegel (h = 2) die Gleichung (3). Damit 
ist gezeigt, da die Dupuitsche Naherung als strenge Lésung fiir ein Festkérpermodell, das in 
vertikaler Richtung unendlich groBe Durchlassigkeit besitzt, anzusehen ist. Eine Schwierigkeit 
bei der Behandlung dieses einachsig widerstandsfreien Festkérpermodells besteht darin, daB bei 
der Ableitung der Gleichungen fiir die Grundwasserbewegung vorausgesetzt wurde, der Widerstand, 
den der Festkérper der strémenden Fliissigkeit entgegensetzt, sei so groBb, dafS bereits bei sehr 
kleinen Strémungsgeschwindigkeiten ein quasistationdarer Zustand erreicht wird. Diese Voraus- 
setzung ist aber hier fiir die vertikale Richtung nicht mehr erfiillt. Daraus ergeben sich bei Vor- 
handensein von nichtvertikalen Trennflachen zwischen Luft und Grundwassertrager gewisse 
Schwierigkeiten, auf die jedoch hier nicht weiter eingegangen werden soll. Sind diese Flachen ver- 
tikal, so treten diese Schwierigkeiten im stationaren Fall nicht auf. Da jedoch langs solcher ver- 
tikalen Flachen ¢h/¢cz = 1 sein muB, so kann bei dem einachsig widerstandsfreien Festkérpermodell 
wegen (9) auf diesen Flachen keine Hangquelle auftreten; die Einmiindung des Grundwasserspiegels 
erfolgt im stationaren und instationdren Fall stets in den Spiegel der ruhenden Fliissigkeit bzw. 
bei vollkommener Absenkung an der undurchlassigen Schicht. lm instationaren Fall sind demnach 
Anfangsbedingungen, die eine Hangquelle beinhalten, nicht zulassig, da diese mit einer quasi- 
stationaren Behandlung nicht vertraglich ist. Aus dem Vorstehenden ersieht man, daB die Aus- 
bildung einer Hangquelle in einer vertikalen Trennflache Luft—Festkérper nur durch den Wider- 
stand des Festkérpers in vertikaler Richtung bedingt ist, und da daher die Anwendung des 
Dupuitschen Ansatzes auf den isotropen Festkérper bedenklich erscheinen muB. 


3. Anwendung der Relaxationsmethode auf Probleme der stationaren und instationiren Grund- 
wasserstrémung. a) Dimensionslose Form der Gleichungen. Um die Relaxationsmethode 
anwenden zu kénnen ist es zunachst zweckmabig, die Ausgangsgleichungen in dimensionslose Form 
zu bringen. Wir beschranken uns hier auf einen homogenen Grundwassertrager, der isotrop oder 
anisotrop sein kann. In letzterem Falle setzen wir einen symmetrischen Durchlassigkeitsaffinor 
voraus (orthotroper Festkérper). Wir gehen aus von den Gleichungen (I, 16b) 


Ch Ch Ch 
Kee Gen + ln Ga + Kee Gea = 0 (12) 
und (I, 30b) 
al, tea ch oh oC ch o¢ 
TORS Fs hess 13 
oo a7 B 
a WY OC a o§ 0§ oh CH) Oty lings der Spiegelflache € = ¢, ( ) 
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die beide fiir = Hauptachsensystem (§, 7, ¢) des Durchlassigkeitsaffinors gelten. Wegen der vor- 
elie 02 peace ts des Grundwassertragers ist n in (13) ortsunabhangig. Im anisotropen 
fee ; : : 
ae os ) durch eine volumtreue Affinverzerrung des Bereiches vermittels der Transformations- 
Bon 


= te g = V 7 ey) ee 
SRO pee tii c= rpmecemees - (14 
VKqyy VK V Ron op 


mit | = key kn k-- in die Laplace-Gleichung 


Bb PRS iF 
ee oe 
iiber. Um die hierin auftretenden Gréfen dimensionslos zu machen, schreiben wir 
Pa ete OnSHyty baal, (16) 
worin H eine beliebige charakteristische Lange ist}. 
Damit geht (15) iiber in 
Ch* Ch* eh* 
OE *2 én 4 act a 0 z (15a) 


Die gleichen Substitutionen fithren wir in (13) ein und erhalten 


att oh*  dh* OC¥ —ah* AC# 


oT ee GE* OE* — 7)* ON* lhanes der Spiegelfliche | ase) 
Die dimensionslose Zeit t* ergibt sich hierbei zu 
Seem: 
eS (13b) 


cE ie 
In einem beliebig abgegrenzten Volumen des Strémungsbereiches befindet sich die Flissigkeits- 
menge V, = ff{fnd&dyd¢. Fiihrt man darin die Substitution (14) ein, so erhalt man Vp = 


{Jn dé dy dé. Dies bedeutet, daB das relative Porenvolumen n bei dieser Transformation invariant 
eibt. 

Die auf die Zeiteinheit bezogene Fliissigkeitsmenge Q, welche ein beliebiges Flachenstiick F 
im Innern des Grundwassertragers durchstrémt, ist gegeben durch? 


Q=fvo-d& =fvuz,dyn dl + fv, dl dé + fu,d& dy. (17) 
F F F, F, 
Setzt man die Beziehungen 
oh oh eh 
0g = — hee ae Yq = — Kan Gy > Daag rae cts) 


die aus (I, 10) und (I, 12) fir das Hauptachsensystem folgen, unter Verwendung der Transfor- 


mationen (14) und (16) in (17) ein, so erhalt man 
Drea oh* oh ni Vesa 
i itn VKur fl a&* Hynde sgh | én* Basie: lags dg* Z =H VK Q*. (17a) 


* * ak 
Ft x se 


oh* 


Bei Vorhandensein eines freien Grundwasserspiegels 2,(x, y,t) gilt langs diesem die Rand- 
bedingung h = z,, wenn die z-Achse vertikal nach oben gerichtet ist. Fallt iberdies eine Haupt- 
achse des Durchlassigkeitsaffinors im Falle eines homogenen orthotropen Festkérpers in die z-Rich- 
tung, ist es zweckmafig, statt (14) eine Affinverzerrung zu wahlen, die die z-Koordinate unver- 
andert laBt. Wir benutzen dann die Transformationen 


pe eee a yir7, ne (19a) 


und 
Ths x= Hx*, AN a Leas PB ma (HLS (19b) 


1 Alle dimensionslosen GréSen sollen durch einen Stern gekennzeichnet werden. 
2 Fe, F,, Fr bedeuten die Projektionen des Flachenstiickes auf die drei Koordinatenebenen. Bei mehr- 


facher Uberdeckung der projizierten Flachen ist die Flache in Teilstiicke mit einfachen Uberdeckungen der 


Projektionen zu zerlegen. 
3 
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Ersetzt man in (12) &7,¢ durch x, y, z und benutzt die Transformationen (19a) und (19b) so er- 
halt man 


PE 5 Within CP amas (20) 
ax%2 | ay*? I Bz*2 
Die Gleichung (13) geht mittels der gleichen Transformationen iiber in 


ack  [ah* Gh* Oz — ah* Oak 


(21) 
ae . —— 7 5 oe + Ae a..* ° 
at* oz” Ox* Ox* dy* oy lings der Spiegelflache 
jedoch ergibt sich hierbei fiir die dimensionslose Zeit t* 
kz 
al Sree (22) 


Da die Transformation (19a) nicht volumtreu ist, bleibt das relative Porenvolumen n hier nicht 
invariant. Man erhalt in diesem Fall 


ee. Php wo Sy hey 
Ves | | n dx dy dz = | [> hed dx dy dz = F[ fe dx dy dz ; 


daher ist 


(23) 


Die Gleichung (17) fiir die auf die Zeiteinheit bezogene Fliissigkeitsmenge Q geht durch Einfiihren 
von (7), (19a) und (19b) tber in 


= — ‘1 = A) lng (e A * a * “ 
Q =— H*Vky, ‘lJ dy* dz* a ck dz* dx* -- [ Gedo dy* 


oy * iz 
Bibs ae FX 
=? Vey, Oe bas (24) 


Fiir den isotropen Festkorper (k,,, = ky, = k,, = k) werden die Transformationen (14) und (19a) 


zu identischen Transformationen, die Gleichungen (13b) und (22) gehen beide iiber in 
i k 
pier m= 
hae: Jani? ee 


und die Gleichungen (17a) und (24) gehen beide tiber in 
os al on dy* de* II pa ds* dae [{ on dx ae H2kQ*. (26) 


* 
ia zy F} 


b) Die Relaxationsmethode fiir den ebenen Fall. Um die Anwendung der Relaxations- 
methode auf ebene Probleme der Grundwasserbewegung zu erlautern, wahlen wir ein spezielles 
Beispiel. Es soll die instationdare Durchstrémung 
eines Dammes (Abb. 2) untersucht werden. Urspriing- 
lich habe der Fliissigkeitsspiegel zu beiden Seiten des 
Dammes die gleiche Héhe z =z, = H. Auf der linken 

TITTTTOTOTI TOTTI Seite werde nun plétzlich die Fliissigkeit unter das 
ME EEN . Niveau z=0 abgesenkt und soll auch spater nie 
i pe Real neat manne ste mehr das Niveau z = 0 erreichen. Die Spiegelhéhe 

auf der rechten Seite werde konstant gehalten. Das 


Damm-Material werde als homogen und isotrop (Durchlassigkeit k) vorausgesetzt. Demnach 
setzen wir 


Eat he ular oir ott) oA Glcp ae a Wy oA piel ge. (27) 


Zur Anwendung der Relaxationsmethode! ersetzt man bekanntlich die zu lésende Differential- 
gleichung durch Differenzengleichungen und iiberzieht die (x, z)-Ebene mit einem Netz von Ge- 
raden? x* = wa*, z* =ya*, wo uw und y ganze Zahlen sind und a = H a* die passend zu wah- 
lende Maschenweite des Netzes bedeutet. Die zweiten Ableitungen yon h* an einem Knoten 


ee a 


* R. V. Southwell, Relaxation Methods in Theoretical Physics, Oxford 1946; D. N. Allen, Relaxation 
Methods, McGraw-Hill 1954. 


* Fur das betrachtete Beispiel wird ein quadratisches Netz verwendet. 


“XVI. Band 1958 Heinrich u. Desoyer: Probleme der stationaren u. instationaren Grundwasserstromungen 35 
,v) des Netzes werden naherungsweise ersetzt durch 


Be gee , 
sat), ge te — 2M Bae) 


eh* 1 (28) 
(a) mre ( ieee —2 hi. » + hs a 1) : 
1,0 


Die ce leichung fiir das ebene Problem geht dann bis auf einen Fehler von der Gro®enord- 
nung a** iiber in 
* 


“iu — 1,» ie he 1,» te er = The ot att 4 Lees = Tuy = ) A (29) 


Wenn in (29) Naherungswerte statt der Werte der gesuchten Lésung eingesetzt werden, so wird 
der Rest Ri» i. a. ungleich Null. Geht man von irgendwelchen Naherungswerten in den Knoten 
pas. so wird (29) zunachst i. a. nicht in allen Knoten erfiillt sein, und es ergeben sich gewisse Reste 
Ry.». Durch systematisches Abandern der h*-Werte an 
den Knoten kénnen diese Reste allmahlich bis zur ge- 
wiinschten Genauigkeit abgebaut werden. Es ist zweck- 
maBig, die jeweils dem Betrag nach gréBten Reste ab- 
zubauen. Andert man den Wert h* in einem Knoten um 
Ah*, so andert sich der Rest R* an diesem Knoten um 
—4Ah*, die Reste in den vier umliegenden Knoten 
andern sich um Ah*, wie aus (29) leicht ersichtlich ist. 
Wenn in gréBeren Gebieten des Netzes die iiberwiegende 
Zahl der Reste gleiches Vorzeichen aufweisen, wendet 
man zur Erzielung rascherer Konvergenz die Verfahren 
des Uberrelaxierens oder der Gruppen- und Blockrela- 
xation an!. td 

Bei Problemen der Grundwasserbewegung sind die Abb. 3. Randstern mit verkiirzten Abstiinden. 
Werte von h* langs des Grundwasserspiegels und langs 
der Trennflachen des Festkérpers gegen Luft oder ruhendes Wasser vorgegeben. Am Rande 
kénnen Knoten mit verkiirzten Abstinden gemaB Abb. 3 auftreten. In diesem Fall ist (29) zu 
ersetzen durch + 


a eee | - * Wa hisse * 20 fb 2\ igs 2 Bee OL Se 
Prt + ay 7 1+ py Bate + 1 + By ioe bn A+ By) © (5; | a a eee nae 


Der Fehler liegt in der GréBenordnung a**. Abinderung des Wertes hj, um Ahj,, andert den 
Rest Ri, nach (29a) um 
2) 9 
AR‘, =—(4- + 4—-\Ahj,. 29b 
Lb? By if Bu LO? ( 9 ) 


Langs einer undurchlassigen Flache muB grad h* senkrecht auf die Flachennormale stehen. 
Bei Vorhandensein von geradlinig begrenzten undurchlassigen Schichten wahlt man zweckmabig 
diese Begrenzung als Netzlinie und kann sich daran das ganze Netz mit allen Knotenwerten und 
Resten gespiegelt denken, wodurch diese Begrenzung wie gefordert automatisch zur Stromlinie 
wird ?, 

Es soll nunmehr die Anwendung der Relaxationsmethode auf das zu Beginn dieses Abschnittes 
erlauterte Beispiel der instationaren Durchstrémung eines Dammes gezeigt werden. Wir denken 
uns den Fliissigkeitsspiegel links vom Damm zur Zeit t = 0 plotzlich unter das Niveau z = 0 
abgesenkt. Fiir diesen Zeitpunkt ermittelt man die Verteilung der Standrohrspiegelhéhe?. Am 
aweckmaBigsten beginnt man dabei mit einem weitmaschigen Netz, in das man, von den gegebenen 
Randwerten ausgehend, geschatzte Werte eintragt und diese durch das Relaxationsverfahren aus- 
gleicht (Abb. 4). Die Werte h* der bezogenen Standrohrspicgelhéhe wurden mit 10* multipliziert, 
um Dezimalzahlen zu vermeiden. Abb. 5 zeigt ein aus Abb. 4 entwickeltes feinmaschigeres Netz‘, 
in das auch Linien konstanter Standrohrspiegelhéhe und konstanten Druckes eingezeichnet sind. 


PA ler. ataeaes One Ss. Olve 


2 Allen, a.a.O. S.72. Bei krummliniger Begrenzung der undurchlassigen Schicht siehe Allen, a. a. O. 
S. 74—85. 

3 Es wird dabei angenommen, daf der Anlaufvorgang, wie er in I, S. 81 erwahnt wurde, bereits abgeklungen 
ist und sich ein quasistationarer Strémungszustand ausgebildet hat. Die Anlaufzeit wird als vernachlassigbar 
klein angesehen. 

eWAllensasa. Ols. Ol, 

oie 
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In Abb. 5 sowie in den folgenden abgebildeten Netzen wurden alle Reste bis auf den Betrag < 2 
herabgedriickt. Aus der Verteilung der h*-Werte laBt sich nun die dimensionslose Sinkgeschwin- 
digkeit — 0z$/@t* in den Knotenpunkten des Spiegels aus Gleichung (13a), die fiir den ebenen, 
isotropen Fall die Form 

dh* —Oh* Get 


O2¥ | 
cap kina * * 
Oz ox™ Ox fir 2* = ae 


Tat 


annimmt, berechnen, wobei t* mit der tatsiachlichen Zeit t durch die Beziehung (25) zusammen- 
hangt. Da die Bestimmung der Differentialquotienten auf der rechten Seite von (30) auf graphische 
Art meist zu ungenau ist, empfiehlt es sich, diese rechnerisch aus Parabeln zweiter Ordnung zu 
ermitteln, die durch je drei Punkte gelegt werden. Diese Sinkgeschwindigkeiten kann man nun 
benutzen, um die Spiegellage zu einem etwas spateren Zeitpunkt zu ermitteln. Dabei wird ange- 


(30) 


nommen, daf sich die Sinkgeschwindigkeit wahrend kleiner Zeitintervalle nicht wesentlich andert. | 


Le a ee 
— 


70000 _-< 10000 1000 70000 


Abb. 4. Ausgangsnetz fiir die Verteilung der Standrohrspiegelhéhen (104+ h*) fiir den Anfangszustand der Sickerstrémung. 


Eine Abschatzung der Anderung der Sinkgeschwindigkeit bei Veranderung der Spiegelflache ist 
von vornherein schwer méglich. Man wird zunachst ein Zeitintervall At* so wahlen, da die 
Spiegelanderung Az3, die nach der Gleichung 


s One e 
Ass = oat At (31) | 


berechnet werden kann, im ganzen Bereich klein bleibt gegen die Maschenweite des Netzes. Nun — 


wird der neue Spiegel eingezeichnet und das Netz wieder abgeglichen. Kin klein gehaltener Absenk- 
schritt wirkt sich hierbei auch insofern vorteilhaft aus, als das neuerliche Abgleichen des Netzes 
wegen der geringen Veranderungen der Knotenwerte am Spiegel entsprechend rasch bewerkstelligt 
werden kann. Nach Abgleichen des Netzes kénnen die neuen Sinkgeschwindigkeiten nach (30) 
berechnet werden. Nun hat man auch die Méglichkeit, diese Werte mit den vorhergegangenen zu 
vergleichen und gegebenenfalls, wenn unzulissig grof{e Anderungen aufgetreten sein sollten, das 
Verfahren mit einem kleineren Wert von At* zu wiederholen. 

Auf diese Weise kann man nun Schritt fiir Schritt den Verlauf der Spiegelabsenkung verfolgen. 
In Abb. 6 ist eine Reihe von aufeinanderfolgenden Spiegellagen eingezeichnet, die auf diese Weise 
bestimmt wurden. Die zugehérigen Zeitintervalle sind zwischen den Spiegellagen eingetragen!. 

Da der Ubergang zur stationaren Spiegellage asymptotisch erfolgt, kann dieses Verfahren nicht 


bis zur Erreichung des stationdéren Zustandes fortgesetzt werden. Wenn die Sinkgeschwindigkeit | 
nach (30) im ganzen Bereich sehr kleine Werte erreicht hat, empfiehlt es sich, von hier aus die 
Lage des stationaren Spiegels nach der Bedingung aufzusuchen, daB dieser iiberall auf den Linien 
h* = konst. senkrecht stehen muf, weil ein stationarer Spiegel Stromlinienflache ist. Dabei ergibt | 


sich auch die Lage der Hangquelle. Abb. 7 zeigt die auf diese Art gefundene stationare Spiegellage 
sowie die Verteilung der Standrohrspiegelhéhen und Linien h* = konst. 
Man ist nun auch in der Lage, die auf die Raumeinheit bezogene Resultante 7 - © der Schnitt- 


krafte im Festkérperskelett zu ermitteln. Diese Grofe ist fiir die Beurteilung der Rutschungsgefahr | 


maBgebend. Nach Gleichung (IT, 9) gilt: 
ee Vee ean MI) AUxapAIs (32) 


Hierin ist © der Spannungstensor fiir das Festkérperskelett, 2 das Verhaltnis des von Fliissigkeit 
bedeckten Anteiles einer Schnittflache zur gesamten Schnittflache (6 und J hangen von der Schnitt- 


fiihrung ab); y;, baw. y,, sind die spezifischen Gewichte von Festkérpermaterial bzw. F liissigkeit ; 


n ist das relative Porenvolumen und U das Potential der Massenkriifte je Masseneinheit. 


' BE. P. Nemecek, a. a. O., hat die instationare Durchstromung eines Dammes im Modellversuch verfolgt. 
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Bei kérniger Struktur des Festkérpers sind zur Beurteilung der Rutschungsgefahr Schnitte zu 


betrachten, die die Kérner des Festkérperskelettes nicht schneiden. Fiir solche ist bei woraus- 


gesetzter Punktberithrung der Kérner 2 — 1. Man kann d fiir n = konst =gz die 
ease dcx ee ann fiir n onst. und U = gz die 
7-6 =VO (33) 

mit 
C= eat hee (33a) 


schreiben, Um ein Bild der Verteilung von /Q zu gewinnen, kann man die Kurvenschar Q — konst. 
mit konstanter Differenz des Kurvenparameters zeichnen. Mit h = ply, + 2 und (27) 1aBt sich 
(33a) in dimensionslose Form bringen 


Q 
ae a eae (33b) 


mit 


x = (1—n) ie 1] . (33c) 
Vw 

In Abb. 8 und 9 wurden die Linien Q* = konst. und die zugehérigen orthogonalen Trajektorien 

fiir den Ausgangszustand bzw. stationaren Endzustand der Strémung eingezeichnet. Dabei wurde 

x = 1,1 gewahlt, ein Wert, der sandigem Material entspricht. Da die Linien Q* = konst. in gleichen 


4Z 
me ap ae Nes bye 
ungestorrer Grundwassersplege/ 
T 
| Vg 
he 
} | 
as 
*s 
AS 
PATS r ee Tt a s 
GITTITTATTATIPTT TT ITIP TOV I RT eee r yi 
undurchlassige Schicht ak oe ie Ta - 
5 e Tt = 
Abb. 10. Skizze zum gewahlten Beispiel des stationiren Zuflusses zu einem Abb. 11. Zur Ableitung der Knotengleichungen fiir 
Brunnen. rotationssymmetrische Netze. 


Intervallen gestaffelt sind, ist ihr Normalabstand ein reziprokes Maf fiir die auf die Raumeinheit 
bezogene Resultante der Schnittkrafte im Festkérperskelett. Wie Abb. 8 zeigt, ist diese Resultante 
zu Beginn der Absenkung in der Umgebung des Punktes 4 am gréfSten, wahrend sie nach Er- 
reichung des stationaren Zustandes (Abb. 9) ziemlich gleichmaBbig tiber den Str6mungsbereich ver- 
teilt ist und nicht mehr so hohe Werte erreicht. 


4, Die Relaxationsmethode fiir den zylindersymmetrischen Fall. Als Beispiel fiir eine zylinder- 
symmetrische Stroémung wahlen wir den stationaren Zuflu8 zu einem vollkommenen Brunnen mit 
begrenzter Reichweite R (Abb. 10). Wir setzen voraus, da die Durchlassigkeiten ortsunabhangig 
sind und fiir alle horizontalen Richtungen denselben Wert besitzen. In vertikaler Richtung kann 
die Durchlassigkeit einen anderen Wert haben. Demgemaf setzen wir in (12) 


oh eae ays =| (34) 
Keg = Reg = Bey 5 Ieee = Iie J 
Durch Transformation auf Zylinderkoordinaten wird dann aus (12) 
oh 1 ch o*h ; 
y = 35 
k,, (a r | kg ez? 0. ( >) 
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und fiihren dimensionslose GréBen h*, r*, z* nach den Gleichungen 


hR=HAh*, r=Hr >. 2S (37) 


ein. Mit (36) und (37) wird aus (35): 
a tae ta 0 - (38) 


or*? if 


Zur Anwendung der Relaxationsmethode iiberziehen wir den Meridianschnitt mit einem Netz von 
aquidistanten Geraden rj, =“ a*, z =v a*, wobei a = H a* die Maschenweite bedeutet. 

Um in einheitlicher Weise die Knotenpunktsgleichungen sowohl fiir Sterne mit gleichen als 
auch ungleichen Abstanden abzuleiten’, entwickeln wir die Funktion n*(r*, z*) im Knotenpunkt 
(u,v) in eine Taylor-Reihe (Abb. 11). Wenn wir nach dem dritten Glied abbrechen, erhalten wir 
fiir die Werte h* in den Nachbarknoten die folgenden Darstellungen: 


: ‘ oe Te ale eae ; 
fe aay == les =p = Fs x = 9 Or *2 b Pas SI 2 ( 9) 
oh* i Gee 
Beet ee ae ie 2 eer rek 40 
éh* 1/@h*\ as 
Iie $1 = Iii» sas wie ala), 3 pede (41) 
oh* LY GA 
Saal tN eae 2 
tee 7 ve aaah: & )., , , ) 
Aus (39) und (40) erhalt man bis auf Glieder vierter, bzw. dritter Ordnung in a* 
oPh* 1 : 
Ga = aa ( ety + bps,» — 2 hyy) (43) 
By 
und 
oh* 1 
a a 2 a* (his. Lupe nee 7) 3 (44) 
Thy? 


aus (41) und (42) ergibt sich bis auf Glieder vierter Ordnung in a* 
O*h* 1 
a] = gag (hyn 41 + Itgy—1 — 2 hin») - (45) 
~ - / by 


Setzt man (43), (44), (45) in die Differentialgleichung (38) ein, so erhalt man nach Multiplikation 
mit a*? bis auf GréBen dritter Ordnung in a* 


( = | ye +(1 ee a] hisaje + hay—1 + hgyti1—4 hi» = Ry, == Oy (46) 
17 


Gleichung (46) muB fiir alle reguléren Knoten erfiillt werden. 
Wird im Zuge des Relaxationsverfahrens der Wert hi, um Ahj,, abgeandert, so andert sich 
der Rest Ri, um 


ARty == 4 ARS, (47) 


Durch Anschreiben der der Gleichung (46) entsprechenden Knotengleichungen fiir die vier 
umliegenden Knotenpunkte erkennt man, daB die obige Anderung von hi, um Ahy, folgende 
Abanderung der Reste in diesen Nachbarknoten erfordert: 


* 
ARS =(1 Ae: Ati (48) 
pie 
° a* 
A Rh + Lon = ( ae yee) Ati ? (49) 
“Tut 
AR ea = AR eae (50) 


Am Rand des Bereiches kénnen Sterne mit héchstens zwei verkiirzten oder verlangerten Abstanden 


auftreten (Abb. 12). In diesem Fall bleiben.(39) und (42) unverandert, wahrend (40) baw. (41) zu 


+ Fur Sterne mit ungleichen Abstanden findet sich beim zylindersymmetrischen Fall im Schrifttum keine 
Knotengleichung, die an Genauigkeit der Formel fiir normale Sterne entspricht. 
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ersetzen sind durch 
oh* lL / @h* 
hehe Slee Br a* + i (a7) (B, a*)?—... (40a) 
»? u,v 


bzw. 


a ah* 1 /@h* ‘ 
hn = hye (| By a" + > (a (oF hea (41a) 
tv Ly? 


A2h * 
Berechnet man aus (39) und (40a) die Werte (| und (| sowie aus (41a) und (42) den 
My? My? 
*h* ; : : : : 
Wert (i) und setzt diese in (38) ein, so ergibt sich nach Multiplikation mit a*2 bis auf GréBen 
A 


dritter Ordnung in a* 


2 l a* he Z re ue h* 2 : 
OI (| arg) rept tare A) Me pp Me 
2 Oe Sera ies, 
ee ; h* | 4 fac = ‘I eee eR a 
a (+ 6y)bu ™ 2/ treba alee Ps )he Sa epee vo. 
Abanderung des Wertes he, um Ah; andert den Rest Ri,» um 


x i 1 G* th==/py _ 

ARG A 2 pag get # ) Abi (47a) 
Die Reste in den anliegenden Knoten im Inneren des Bereiches andern sich dabei gemaB (49) 
und (50), wenn diese Knoten Mittelpunkte von regularen 
Sternen sind. 

An Hand des eingangs erwahnten Beispiels (Abb. 10) 
soll nun kurz gezeigt werden, wie die Lage des statio- 
naren Grundwasserspiegels beim ZufluB zu einem Brun- 
nen ermittelt werden kann. Der stationdre Zustand kann 
immer als jener Grenzzustand aufgefaBt werden, dem 
eine beliebige instationaére Ausgangsstrémung asymp- 
totisch zustrebt. Man wird daher am einfachsten von 
einer auf irgendeinem Wege gewonnenen Naherungs- 
lésung fiir die stationdre Spiegellage ausgehen. Fiir diese 
wird das Netz ausrelaxiert und daraus analog wie im 
vorigen Beispiel die noch vorhandene Geschwindigkeit 
des Spiegels aus der Beziehung 


* 
Zz 


| 
| 
| 


a 6 s Abb. 12. Randstern mit verktirzten Abstinden 
Oz, oh* oh* Oz (5 1) im rotationssymmetrischen Netz. 


fae az* | @r* Or*|_« 


ermittelt. Gleichung (51) folgt aus (21) durch Transformation auf Zylinderkoordinaten. t* ist 
dabei durch (22) definiert. Analog wie im vorigen Beispiel kann die Annaherung des Spiegels an 
die stationare Endlage in einigen Schritten verfolgt werden. Die endgiiltige Festlegung des statio- 
naren Spiegels kann schlieBlich wieder mit Hilfe der Bedingung erfolgen, dafs dieser auf den Linien 
h* = konst. senkrecht stehen muB. Abb. 13 zeigt einen auf diese Art ermittelten stationdren 
Spiegel’. Es sind die Linien h* = konst. eingezeichnet. Die in der Zeiteinheit dem Brunnen zu- 
flieBende Wassermenge kann nun leicht ermittelt werden. Die Fliissigkeitsmenge Q, die in der 
Zeiteinheit durch eine Zylinderflache mit dem Halbmesser r hindurchtritt, ist 


Oee2 nr ke. Qed. (52) 
2=0 
Mit (36) und (37) ergibt sich daraus 
= 
Q@ _n*« | h* 4x 
kg =O = 2m | Geet (53) 
2*=0 


1 Fiir fiinf verschiedene Absenkungen des Wasserspiegels im Brunnenschacht hat H. Schmidt nach der 
hier entwickelten Methode die stationare Lage des Grundwasserspiegels und die zugehérigen Entnahmemengen 
ermittelt. (Uber eine Anwendung der Relaxationsmethode zur Behandlung von Grundwasserstromungen, 


Diss. T. H. Wien, 1956.) 
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Die GroBe Q* muB, im stationaren Fall, von r* unabhangig sein. Dies kann als Kontrolle beniitzt 
werden. Das obige Integral kann, etwa nach der Simpsonschen Regel, leicht ausgewertet werden. 
Fiir das durchgefiihrte Beispiel ergab sich Q* = 0,864. 

Die numerischen Berechnungen zu den angegebenen Beispielen wurden mit einer vollauto- 
matischen Rechenmaschine durchgefiihrt. Es ware jedoch empfehlenswert, fiir Aufgaben dieser 
Art elektronische Rechengerate einzusetzen. Man koénnte in diesem Fall entweder die Differenzen- 
gleichungen fiir die einzelnen Knoten als lineares Gleichungssystem einer Lésung zufiihren oder 
man benutzt zur Lésung der Potentialgleichung die sog. Monte-Carlo-Methodet. Diese hat den 
Vorteil, da& man mit ihrer Hilfe den Wert der Potentialfunktion in einem beliebig vorgegebenen 
Knoten des Bereiches ermitteln kann, ohne daB es nétig ist, die Potentialwerte in den tibrigen 
Knoten zu berechnen. Da fiir die Ermittlung der Absinkgeschwindigkeit des Spiegels nur die Ab- 
leitung der Potentialfunktion am Spiegel selbst benétigt wird, kénnte man auf diese Weise die 
Berechnung der Werte im Inneren des Bereiches vermeiden, sofern sie nicht aus anderen Griinden 
von Interesse sind. 


5. Zusammenfassung. Im AnschluB an zwei vorangegangene Arbeiten, in denen Grundlagen 
fiir die Behandlung von stationaren und instationaren Grundwasserstrémungen entwickelt wurden, 
wird hier die Anwendung auf praktische Probleme der Grundwasserstrémungen gezeigt. Nach 
einer kritischen Betrachtung der Dupuitschen Naherungslésung werden die Grundgleichungen 
der instationéren Grundwasserbewegung zur Anwendung der Relaxationsmethode auf dimensions- 
lose Form gebracht. Am Beispiel der instationaren Durchstrémung eines Dammes wird die An- 
wendung der Relaxationsmethode auf ebene Probleme erlautert und zugleich gezeigt, wie man sich 
ein Bild von der Beanspruchung des Festkérperskelettes verschaffen kann. Als Beispiel fiir die 
Lésung eines zylindersymmetrischen Problems mittels der Relaxationsmethode wird der stationare 
ZufluB zu einem Brunnen als Grenzfall, dem eine instationére Strémung asymptotisch zustrebt, 
behandelt. 


(Eingegangen am 25. Marz 1957.) 


Anschrift der Verfasser: Prof. Dr. G. Heinrich und Dr. K. Desoyer, Wien IV, Karlsplatz 13, Technische Hoch- 
schule 


1 Siehe z. B. Ed. F. Beckenbach, Modern Mathematics for the Engineer, S. 293 ff., McGraw-Hill 1956 (dort 


weitere Literaturangaben). 


XXVI. Band 1958 1 rs . tome : 
Weidenhammer: Auswanderungserscheinungen in SchwingungsmeBgeraten 43 


Auswanderungserscheinungen in Schwingungsmefgeriten 


Von F. Weidenhammer 


iv Fragestellung. Bei SchwingungsmeBgeraten und Anzeigegeraten, deren mechanische 
Systeme aus einem Drehschwinger bestehen (Abb. 1), ist eine unerwiinschte Auswanderungs- 
erscheinung seit etwa 20 Jahren bekannt!. Unter der Wirkung von Erschiitterungen kénnen diese 
Drehschwinger erzwungene Schwingungen um Mittellagen ausfiihren, die nicht die vorgesehenen 
urspriinglichen statischen Ruhelagen sind. Das MeBsystem wandert aus dem Sollwert seiner Mittel- 
er ee ee eh eucian ee ee cde 

Bee Kraftespiel ein- 
greifen. Bei standig wachsender Erschiitterungsamplitude und -frequenz kann schlieBlich ein 
Moment der Tragheitskrafte mit einem zeitlichen Mittelwert von so groBem Betrage erzeugt 
werden, da8 der Drehschwinger vollstandig in die Erschiitterungsrichtung hineingezogen wird, 
falls nicht Anschlage den GréBtausschlag zuvor begrenzen. Fiir das erschiitterte Pendel tcehware: 
feld ist diese Erscheinung aus zahlreichen theoretischen und einigen experimentellen Untersuchun- 
gen seit langem bekannt®. Die Ergebnisse dieser Arbeiten kénnen allerdings nicht unmittelbar fiir 
SchwingungsmeBgerate und Anzeigegerite verwertet werden, da bei diesen Geraten das Riicksteli- 
moment durch eine Feder und nicht durch das Schwerefeld allein geliefert wird. 


Erschlitferungs- 
riehtung 
Abb. 1. Gedampfter Schwingweggeber mit Spiralfeder 
und Ausschlagbegrenzung. Abb, 2. Der Drehschwinger im Koordinatensystem, 


Diese geschilderte Art des kinetischen Versagens ist bei Geratezeigern (wie z. B. den Milli- 
amperemetern) nur durch die Erschiitterung als Stérung bedingt und k6énnte an sich durch eine 
geeignete Schwingungsabschirmung stets behoben werden. Praktisch findet diese Méglichkeit der 
Entstérung allerdings auf Grund der meist beengten Platzverhaltnisse schnell eine Grenze, so dal 
auch in diesem Falle eine genaue Kenntnis der Auswanderungsvorgange wichtig ist. Bei Schwin- 
gungsmeBgeraten mit einem Drehschwinger ist grundsatzlich keine Abhilfe méglich, da das Ver- 
sagen im mechanischen Prinzip des MeSsystems begriindet liegt. Nur der Ubergang auf ein trans- 
latorisch arbeitendes MeBsystem kénnte hier das Auswandern verhindern. Da jedoch Schwingweg- 
meBgerate nur dann eine verzerrungsfreie Anzeige liefern, wenn sie tief abgestimmt sind, mu man 
Gerate mit méglichst kleinen Gerateeigenfrequenzen baucu. Aber gerade diese unumgangliche 
Forderung an das MeBgerat verlangt konstruktiv fast immer den Ubergang auf drehschwingende 
Horizontalpendel mit groBer Drehmasse und weichen Spiralfedern. Im weiteren wird daher nur 
von Schwingungsmefgeraten dieses Types die Rede sein, obgleich alle folgenden Ergebnisse auch 
fiir Anzeigegerate mit einem drehschwingungsfahigen Systcin der gleichen Art gelten. 

Das Ziel der vorliegenden Untersuchung ist die Berechnung der Mittellagen der eintretenden 
erzwungenen Schwingungen abhangig von Frequenz und Amplitude der Erschiitterung und den 
gegebenen mechanischen Geratedaten. Diese Mittellagen werden sich im allgemeinen in der un- 
mittelbaren Nahe der statischen Ruhelage finden. Bei wachsenden Erschiitterungsintensitaten 


1 F. vy, Grundherr und K. Klotter, Jb. dtsch. Luftf.-Forschung 1938, 5. 631. % 
2 Zusammenfassende Darstellung bei K. Klotter, Technische Schwingungslehre, 2. Aufl., S. 376, Berlin/ 


Géttingen/Heidelberg: Springer 1951. 
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kénnen jedoch die erzwungenen Schwingungen um diese Lagen instabil werden, so daB ersatzweise 
Schwingungen um ganz anders gelegene Mittellagen eintreten. Die Berechnung der méglichen 
erzwungenen Schwingungen ist daher durch eine Untersuchung ihrer Stabilitat zu erganzen. 


2. Die Bewegungsgleichung. Als Bewegungsgleichung eines Drehschwingers mit Federriick- 
stellung (Abb. 1) hat man streng eine partielle Differentialgleichung fir die Schwingungen des 
Federkontinuums aufzustellen und den angehangten starren Pendelkérper und seine Bewegungs- 
méglichkeit in einer nichtlinearen zeitabhangigen Randbedingung zu beriicksichtigen. Da jedoch 
eine derartige Formulierung des Problems kaum Aussicht auf eine Lésung bietet, werden ent- 
sprechend den Gepflogenheiten der Schwingungslehre die Tragheitseigenschaften der Feder ver- 
nachlassigt. Dann bleibt ein Drehschwingungsproblem mit einem Freiheitsgrad bestehen, und die 
entstehende Bewegungsgleichung wird eine gewohnliche nichtlineare Differentialgleichung, deren 
Lésung naherungsweise méglich ist. 

Um die Bewegungsgleichung fiir einen Schwinger gemaB Abb. 2 aufzustellen, werden die 
Lagrangeschen Gleichungen herangezogen. Als Koordinate wird der Winkelausschlag g(t) gegen 
die Horizontale benutzt. Die Erschiitterung des Gerates wird durch q(t) in ihrer Weg—Zeit-Ab- 
hangigkeit beschrieben und erfolge geradlinig unter dem Winkel 6 gegen die Horizontale. Um dann 
zunichst die kinetische Energie T angeben zu kénnen, hat man das Integral 


T=5 | b2 dm . (1) 
(™) 
iiber alle Massenelemente dm des Pendelkérpers auszuwerten. Zu dem Zweck werden voriber- 
gehend eine Koordinate € auf dem Pendelkérper und die beiden Achseneinheitsvektoren t, j ein- 
gefihrt, so daB sich der Geschwindigkeitsvektor ) aus Relativ- und Fihrungsgeschwindigkeit zu 
b =i(—é sing + qeos 6) +1 (Epcos + q sin 0) 


ergibt. Damit wird aus (1) 


T= [dm + sino) pa | Edm +% | am, 


wofir sich mit der Drehmasse J, der Gesamtmasse M und dem Schwerpunktsabstand s vom Dreh- 
punkt schreiben laBt: 


ee bes es 
T= 7¢+sMsn(GO—9)-Gg7>@. (2) 


Die potentielle Energie besteht aus einem von der aufgezogenen Spiralfeder herriihrenden Anteil 
und der Arbeit, die bei Bewegung des Pendelkérpers gegen das Schwerefeld zu leisten ist. Es gilt 
also 


V=V,+-,(~—or) + Messing, (3) 


wenn V, einen belanglosen Energiebetrag bezeichnet, der weder von g noch von @ abhangt. Fiir 
die Spiralfeder wird in (3) das Hookesche Gesetz als geltend angenommen, so daf c die Drehfeder- 
konstante bezeichnet und gy, den konstanten Winkel angibt, fiir den die Drehfeder spannungslos 
wird. Dieser Federaufzugswinkel bestimmt sich aus der Forderung, daB das Pendel die horizontale 
Lage als statische Ruhelage haben soll. Aus der Momentengleichgewichtsbedingung 

0 

= =e (p—¢r) + Mgscosp =0 


ergibt sich fir diese Ruhelage @ ¢,, = 0 der erforderliche Federaufzugswinkel zu: 
_Mgs 


iM (1) 


In den Formeln wird jedoch weiterhin gy, nicht durch den Sollwert (4) ersetzt, da sich etwaige 
durch Einbaufehler bedingte Abweichungen von diesem Sollwert als sehr wesentlich fiir die Sta- 
hilitatsaussagen erweisen werden. 


Die bei Schwingungsmefgeraten angestrebte geschwindigkeitsproportionale Dampfung wird 
durch das Moment 


Yr 


0—— be (5) 


mit der Dampfungsmomentenkonstanten b in die Rechnung aufgenommen. 
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Die Lagrangesche Gleichung 


mit dem aus (2) und (3) zu bildendem kinetischen Potential L — T — V liefert nun mit dem Dimp- 
fungsmoment (5) die Bewegungsgleichung des Drehschwingers in der Form 


JG +bp+c(p—o,z) + Mgscosp +s MGsin(’—g) =0. (6) 
Die Erschiitterungsfunktion q(t) wird darin als periodisch mit der Erregerfrequenz j und der 
Amplitude q) vorausgesetzt. Dann ist 
q(t) = qo L(t) , (7) 
worin fiir rein harmonische Erschiitterung einfach J7(t) = cost 2u setzen ist. Im allgemeinen 
Falle tritt jedoch in (7) fiir J7(t) eine Fourierreihe auf. 
Nach Division durch die Drehmasse J erkennt man in der Schwingungsgleichung (6) in 


op = VE (8) 


die Eigenkreisfrequenz fiir ungedampfte kleine Schwingungen des Pendels um die horizontale 
Gleichgewichtslage. Diese Frequenz ist als Gerateeigenfrequenz besonders wichtig. Mit der redu- 
zierten Pendellange / tritt eine weitere Kreisfrequenz 


eMs ah 
oY oF ° 


auf. Dies ware die Kigenkreisfrequenz fiir den Drehschwinger allein als Schwerependel in vertikaler 
Lage ohne Federriickstellung und ohne Dampfung. Das Erregerglied mit der Erschiitterungs- 
beschleunigung g erhalt mit der Amplitude g, aus (7) und der reduzierten Pendellange / die relative 
Erregeramplitude 
q 

Baars (10) 
Diese relative Amplitude ist fiir die vorliegende Problemstellung stets sehr klein, so da mit 
é<_ ein kleiner Parameter in der nichtlinearen Schwingungsgleichung vorkommt. Zur Kenn- 
zeichnung des Dampfungsmomentes wird das dimensionslose Lehrsche DampfungsmaB verwendet: 


re beet (11) 


Da man bei SchwingungsmeBgeraten die fiir eine treue Anzeige giinstigste Dampfung D ®& 1/2 
anstrebt, hat (11) einen Wert von der GréBenordnung eins. Da Auswanderungserscheinungen erst 
bei Erschiitterungsfrequenzen w auftreten, die wesentlich tiber den Eigenfrequenzen w, und W, 
liegen, sind die mit (8) und (9) gebildeten Frequenzverhaltnisse 
On @ 

Xp = a > Xo —= = (12) 
besonders wichtige GréBen, die ebenfalls klein gegen eins sind. Mit den Abkirzungen (7) bis (12) 
und mit einer dimensionslosen Zeitzahlung 


tT=ot (13) 
ergibt sich aus (6) die grundlegende nichtlineare Schwingungsgleichung 
gp’ +2 Dapy' +ah(p — gx) + a2 cos p -+ ell (zr) sin (6 — gy) = 0, (14) 


in welcher Akzente nun Differentiation nach der dimensionslosen Zeit t bezeichnen. Hieraus geht 
die aus der Theorie der Schwingungsmefgerate gelaufige lineare Differentialgleichung der erzwun- 
genen Schwingungen hervor, wenn man sich auf kleine Schwingungen um die horizontale Gleich- 
gewichtslage und auf rein harmonische Erschiitterung in der Vertikalen beschrankt. Dann ist 
e II" sin (2/2 —y) = —é cost cosy ~ —ecost, und fiir die horizontale statische Ruhelage 
gp = 0 folgt aus (14) — ae gp + «2 = 0, was genau der Bedingung (4) entspricht. Also bleibt dann 
mit der dimensionslosen Zeitzahlung (13) nur die bekannte lineare Schwingungsgleichung 


gl +2 Dapy tary =e cost 


bestehen, woraus jedoch der beobachtete mehrdeutige Schwingungszustand nicht mehr berechnet 
werden kann. Man ist daher zum Riickgriff auf die nichtlineare Gleichung (14) gezwungen. 
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Die Durchrechnung des Problems mit der Gleichung (14) allein setzt voraus, daB die Beschran- 
kung auf nur einen Freiheitsgrad durch Vernachlassigung der Tragheitswirkungen der Riickstell- 
feder naherungsweise zulassig ist. Diese Vernachlassigung etwaiger Resonanzeigenschaften der 
Feder ist durchaus problematisch, und K. Klotier! hat durch ein Beispiel eines elastischen Schwin- 
gers mit erschiitterter Drehachse auf diese Tragheitswirkungen hingewiesen. Doch soll in der vor- 
liegenden Arbeit von den Tragheitswirkungen abgesehen werden, wie dies in der Schwingungslehre 
meist iiblich ist. Dies hat insofern eine gewisse Berechtigung, als in einer etwas andersartigen 
Untersuchung? eines Gerates mit Schraubenfedern die Schwingungseigenschaften der Feder durch 
Erweiterung der Rechnung auf zwei Freiheitsgrade beriicksichtigt wurden und dabei fiir die In- 
stabilitatsgrenze bei héherer Erschiitterungsfrequenz genau die gleichen Werte berechnet werden, 
wie sie sich aus der Gleichung (14) ergeben werden. Daher erscheint die vorliegende Formulierung 
durchaus geeignet, wesentliche Eigenschaften der Auswanderungserscheinung zu untersuchen. 


3. Die erzwungenen Schwingungen. Zur Berechnung der durch die Gleichung (14) beschrie- 
benen erzwungenen Schwingungen mit hoher Erregerfrequenz eignet sich besonders ein direkter 
Fourier-Ansatz in der von K. L. Stellmacher® angegebenen Weise, bei welcher sich die Koeffizienten 
aus den nichtlinearen Gleichungen rekursibel berechnen lassen. Da die vorliegende Gleichung (14) 
besonders durch die Einbeziehung der Dampfung komplizierter gebaut ist als die bereits behan- 
delten Beispiele vom Typ ungedampfter Pendelschwingungen, wird auf Konvergenzuntersuchungen 
ganz verzichtet und daher auch in formal etwas einfacherer Weise vorgegangen. Zur Berechnung 
der erzwungenen Schwingungen aus Gleichung (14) wird fiir den Schwingungsausschlag y(t) eine 
Fourier-Reihe in der Form 

Cc 
Y(t) =% + P(t) mit Gg = konst. und Pi Dante Oe (15) 
n = — co 

angesetzt. In (15) bezeichnet g den Hauptbestandteil des zeitunabhangigen Gliedes der komplex 
angesetzten Fourier-Reihe. Der in @(t) nochmals vorkommende nullte Fourier-Koeffizient cy 
soll nach einem Vorschlag von K. L. Stellmacher nur als Korrektur zu gy berechnet werden. Dem 
Koeffizienten @. kommt besondere Bedeutung zu, da er bis auf kleine Korrekturen die Mittellage 
der eintretenden Schwingung beschreibt. Ist namlich der Winkel g von dem Winkelausschlag 
der statischen Gleichgewichtslage verschieden, so ist der Schwinger ,,ausgewandert‘“*. Da man 
diese Mittellage bei jedem Experiment in erster Linie beobachtet, wird gy zur Kennzeichnung der 
eintretenden Schwingung verwendet. 

Zur Durchfiihrung der Rechnung soll zunachst der wichtige Fall der rein harmonischen Er- 
schiitterung und dann ein spezieller Fall einer allgemeineren Erschiitterungsfunktion behandelt 
werden. 

a) Harmonische Erschiitterung. In (14) ist nun //'’ = — cost einzusetzen, so daB zu- 
nachst fir Schwingungen um die noch unbekannte Mittellage g) die durch Entwicklung nach @ 
entstehende nichtlineare Schwingungsgleichung gilt 


IT ¢ / 2 | 2 : oD 2 ; P 
O" +2 Doa«,D' + oF (® + % — Hp) + 02 |cos Y — siD Po 71 — ©°S Po x + sin Py 37+: | 
: ce pe 3 
+ € cos T |sin (% — 0) + cos (% — 0) in sin (%— 0) 31 7 69s (Y— 9) oe : | = .0..0( 1460) 


Tragt man hier den Fourieransatz fiir @(t) ein und vergleicht die Faktoren von e‘"* fiir festes 
N £~ 0, so findet man die Koeffizientengleichungen 


I 


2 * 6 : T fee 2 a Be 
— N® ey +12 Doay N cy + OF Cy + Of |— sin Op Cy — COS My 


F It 
; SIN Mo % DS, >, CN—m—n Cm CO) Onc 


m n 


: : 1 
+ COS (Po a 0) (en—1 Bla CN t 1) aa (%— 0) Dy o> (ene oe —] ifs CN—m | 1) Cm 
ven m 
— es (M—?) = 2d Dene ae EON m—n + 1) Cm Cn ar ei aE) 
(DN sea isnt a NS SI 


‘ Vel. S. 356 des in FuBnote 2 auf S. 43 genannten Buches. 
* FE’. Weidenhammer, Z. angew. Math. Mech., 37 (1957) S. 286. 


_ 3 KL, Stellmacher, Z. angew. Math. Mech. 34 (1954) S. 105. (Die Lektiire dieser Arbeit ist durch zahl- 
reiche Druckfehler sehr erschwert.) 
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in denen die Summationen iiber alle positiven und negativen ganzzahligen Werte der Indizes zu 
erstrecken sind und e, , das Kroneckersymbol bezeichnet: 


Conn tee LUE Tt = ON Coe ee Oe MUL Ut SN. 


Fiir den nullten Fourierkoeffizienten q+ cy findet man aus (15, N) fiir N = 0 die besonders 
wichtige Gleichung 


pas, hoe 
m mon 


2 vee 2 i \ i ; 4 1 : 1 : 
Xb (Po + %— Fr) +93 [eos Po — 8in Y  — 008 Gy J) Cm em + sin Gy SS e—m—n lm ey fo 


“ F i 
= ae Joos (% — 6) (c_, + ¢,) — sin (po —- 0) ra re eee Pte Gn oe a) ey 


1 a 
egos (Po — 0) = 2 Pa (Con —n— 1 Al Cm —— nk 1) Cm Cy 


m nm 


0, (15,0) 


aus der sich die Mittellage durch g) und cy) berechnen libt. 

Um das nichtlineare Gleichungssystem fiir die Fourierkoeffizienten gy, und cy(N =0, + 1, + 2,...) 
rekursibel auflésen zu kénnen, werden die Koeffizienten nach dem kleinen Parameter ¢ von (10) 
entwickelt. Zu dem Zweck sind Voraussetzungen iiber die GréBenordnungen der itbrigen Gleichungs- 
parameter von (14) unumganglich. Es gilt also eine Beziehung iiber die GroBenordnungen der 
Frequenzverhaltnisse «1, und a, zu ¢ festzulegen. Hierfiir greift man zweckmaSig auf die Unter- 
suchungen des erschiitterten Pendels im Schwerefeld zuriick. Bei diesen hat sich stets ergeben!, 
daB Auswanderungen aus der statischen Ruhelage dann eintreten, wenn die Erschiitterungs- 
amplitude gy und die zugehérige Frequenz ~ so gro® sind, daB der Intensitatsparameter 


von der GréBenordnung eins ist. Also mu das Frequenzverhaltnis x) von der GréBenordnung ¢ 
sein. Diese Voraussetzung iiber die GréBenordnung wird sich auch fiir das vorliegende, abgeanderte 
Problem als richtig erweisen. Da im allgemeinen bei MeBgeraten «,. und m, die gleiche GréBen- 
ordnung haben, sind auch die Frequenzverhaltnisse von der gleichen GréBenordnung ¢, was durch 
die Schreibweise 

oS yee op = Op € (16) 
zum Ausdruck gebracht werden soll. Fiir den Dampfungsbeiwert gilt dann entsprechend 

2K 12 DR ee (17) 


da 2 D fir die in MeBgeraten angestrebte giinstigste* Dampfung die GréBenordnung eins hat. 
Die Fourierkoeffizienten c, sollen nun ebenfalls nach Potenzen von ¢ entwickelt werden: 


e =ceM 41 Me4 24... (Oe Le 2) cn (18) 


Dabei soll in den Entwicklungen (18) fiir n = 0 das von ¢ unabhangige Glied c fehlen, da hierfiir 
seiner besonderen Bedeutung wegen das Mittellagenglied q mit der Gré®enordnung eins schon 
in den Ansatz (15) eingetragen wurde. Dann bietet (18) die Méglichkeit, cy als Korrektur zur 
Mittellage g) durch 
c= cM) ee c®) i 

zu berechnen. 

Geht man mit den Entwicklungen (18) in die Koeffizientengleichungen (15) ein, so kann man 
diese Gleichungen mit den Ansatzen (16) und (17) rekursibel lésen. Zunachst erhalt man aus den 
Gleichungen (15, N) fiir die Potenzen ¢° 


— N?2c% —0 mit den Lésungen cQ= 0, ING aie Siena eae as (19) 
Der Vergleich der Vorfaktoren der Potenzen ¢! liefert mit (19) die Gleichungen 
= N2 of) a sin (~ — 6) (e1, w + e—1, n) = 9, 


so daB sich die Lésungen 
le ( é < 
pies ols > sin (= 0). t= 0 SUN 2s eos) (20) 


1 Vel. S. 379 des in FuBnote 2 auf Seite 43 genannten Buches von K. Klotter. Der dort verwendete Er- 
schiitterungsparameter ist bis auf den Faktor 2 gleich dem Kehrwert des hier verwendeten Intensitatspara- 
meters. 
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ergeben. Durch Koeffizientenvergleich findet man fiir die Potenzen ¢* mit c® — 0 die Gleichungen 
—N* e8) +42 Diy Nc) +4 cos (yp —8) (a+ Hex) =0 (N=41,42-.), 


deren Lésungen sich mit (20) wiederum unmittelbar angeben lassen: 


‘aah hailed } 

c®) =i D Xp sin (py — 6) + ace (Yo — 9) ef? | 
a 1 

c®), = —i D&gsin (py — 0) + % 008 (Po — 0) cf), | (21) 
lh . 

e8), = ce = a5 el 2 (% — 9) ; 

Ce ==) far IN tale oy eae J 


Grundsatzlich kann man die Berechnung der Fourierkoeffizienten cy (N ~ 0) auf diesem Wege 
bis zu beliebigen Potenzen von ¢ weiter fortsetzen. Doch geniigt es fiir die vorliegende Aufgaben- 
stellung, die Rechnung hier abzubrechen. im 

Die noch verbleibende Gleichung (15, 0) fiir den nullten Fourierkoeffizienten liefert auf dem 
gleichen Wege der Entwicklung fiir die Potenzen ¢° und ¢! nur leere Aussagen. Fiir ¢? hingegen 
findet man die besonders wichtige Beziehung 


= it 
OF (Po — Pr) + BG cos py + J C08 (Y — 8) (c, + ef?) = 0, 


wobei beachtet wurde, daB fiir den Sollwert my, nach (4) auch 


Mgs 8 a2 


Yr 


c Op ae 
geschrieben werden kann, so daB y, als von der GréBenordnung eins erkannt werden kann. Kleine 


Abweichungen von diesem Sollwert werden sich als wesentlich erweisen, andern aber diese Gréfen- 
ordnung nicht. Mit den Werten c®, aus (20) folgt die wichtige Mittellagen-Gleichung 


= = tae 9< 
M(o) = “F (Yo — Pr) + 0G cos Po + ge 2 (% — 6) = 0. (22) 


Aus dieser transzendenten Gleichung lat sich ~) und damit die Mittellage der eintretenden er- 
zwungenen Schwingung bis auf kleine Korrekturen c( berechnen. 

Die Gleichung (22) soll daher noch in eine fiir die Zahlenrechnung geeignete Form gebracht 
werden. Zu dem Zweck wird der Federaufzugswinkel y, mit seinem Sollwert (4) und einem Fehler A 
in der Form 


a2 
2 ar (23) 
in (22) eingesetzt und der abgespaltene Faktor ¢? wieder zugefiigt, so daB nach Division durch «2 
at. lee ‘ 
~~ (Po — A) + cos yy —1 +6 sin 2 (py) — 6) =0 (22, 1) 
0 
entsteht. Durch 
8 gw gw” 
6 a2, Can Samed Ce) 


ist hier der schon erwahnte Intensitatsparameter zur Kennzeichnung der Erschiitterung eingefiihrt 
worden, der sich als das geeignete MaB zur Beurteilung der Auswanderungen erweisen wird. Der 
Parameter ¢ entspricht dem von K. Kloiter! bei der Untersuchung des Pendels mit oszillierendem 
Aufhangepunkt verwendeten Erschiitterungsparameter, der sich also auch bei dem vorliegenden 
andersartigen Problem als wesentlich erweist. Die iibrigen Gleichungsparameter of/o2 = w/w? 
und A sind Geratekonstanten, wihrend der Winkel 6 die Lage der Erschiitterungsrichtung gegen 
die Horizontale angibt. Im allgemeinen hat man fiir jedes Gerat mit seinen festen Daten fiir jeden 
Wert ¢ einen oder mehrere Werte von @ durch Auflésung der transzendenten Gleichung (22,1) 
zu berechnen. Lediglich fiir den Sonderfall eines nullfehlerfreien Gerates (A = 0) und fiir den in 
der Praxis haufigsten Fall der lotrechten Erschiitterung (6 = 90°) kann eine wichtige Lésung von 


* Vgl. S. 379 des in FuBnote 2 auf Seite 43 genannten Buches. 
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(22,1) sofort angegeben werden. Aus 

Or = 

aq Po t cos pH —1— = sin2y, = 0 (22,0) 
folgt namilich in diesem Falle sofort die Lésung Po = 9, welche allerdings nur besagt, daB das 
Gerat in der gewiinschten Weise Schwingungen um die statische Ruhelage Y) = 0 ausfiihren kann. 

Die weitere Auswertung von (15,0) beziiglich der Potenzen ¢® fihrt auf 

a es 1 , , Ls 
age) a OG SIN Po cM) a8 2 Coe (Po aes 0) (c®), ig c®)) sage sin (Po i= 0) pak (lay i ee.) CDiE=i0 ’ 


m 
m 


woraus sich mit den Werten c{), c® aus (20), (21) die erste Korrektur cl) zur Mittellage gy bestimmen 
laBt. Man findet zunachst 


= ee 1 er : : 

ape) — a2 sin gp, c + z Sin (Y — 0) cos (gy — 0) cf!) _ ea sin (~ — 0) sin (@) — 6) ce) = 0 
und daraus endlich 

= 25 8 1 ; 
cM) lag — 32 sin gy -+ 5 cos 2 (M — 9)| = c) M'(y,) = 0. (25) 

Wenn also y als Lésung yon M(y) = 0 nach (22) bestimmt wird und nicht auch die erste Ab- 
leitung M'(y,) fiir diese Wurzelwerte verschwindet, dann ist c{) =0 und eine Korrektur zu q 
ist mindestens von der GréSenordnung ¢*. Allerdings ist ihre Berechnung langwierig und die 


Rechnung soll daher nur so weit getrieben werden, da man ersehen kann, daB der Wert c?) nicht 
verschwinden muB. Aus (15,0) folgt fiir die Potenzen ¢! 


= Sse = 1 1 
Xp cP) — o2 sin gy c®) — a2 cos Mo = yO OM + = cos (% — 6) (c®), + c®)) 


hor. : 5 ’ , 
Gag A sim (Po Fe 0) DS (Cees SF Be ea) cl) =n Cos hes Ge co | ) c®)] 
1 
+ 12 cos (Yo ey 0) iy wa (c,, n 1 ae e)_, n-+ a cf) cf) = 0 ° 


m n 


worin wiederum c®) = 0 fiir alle N unmittelbar beriicksichtigt wurde. Mit den Werten (20) und (21) 
gewinnt man daraus nach einiger Rechnung 


sail= ae Je Ck : : 1 ‘ : ; 
c®) |x — a? sin gy — a. sin? (~) — 0)| — ~ COS ‘Py sin? (y, — 0) + =z 608 (p — 6) (c®, + c®)) 


— 5 sin® (% — 6) sin 2 (yp) — 0) = 0, 
was nur mit Kenntnis der c®) und c®), allgemein ausgewertet werden kann. Da diese Koeffizienten- 
anteile jedoch nur durch langere Rechnung aus (15, N) fiir ¢? gewonnen werden kénnen, soll c?) 
nicht allgemein, sondern nur fiir den Spezialfall eines nullfehlerfreien Gerates mit lotrechter Er- 
schiitterung und auch nur fiir yg, = 0 angegeben werden. Dann kann man die Korrektur c?) ohne 
Kenntnis von c® und c®), berechnen, da der Vorfaktor dieser Koeffizienten cos (py — 6) = cos 2/2 
verschwindet. Man findet dann c® 4 0 aus 


2) (ai = a) oe) A; 


sofern nur der schon in (25) aufgetretene Faktor M'(y, = 0) =a~— 1/2 4 0 ist. Die Bedeutung 
dieser Ungleichung als Stabilitatsaussage wird sich spater zu erkennen geben. 

Die Auswertung der Fourierkoeffizientengleichungen (15,N) und (15,0) kann auf dem an- 
gegebenen Wege bis zu beliebiger Genauigkeit vorgetrieben werden, wobei allerdings die Kon- 
vergenzfrage offen bleibt!. Beschrankt man sich auf quadratische Potenzen von ¢, so sind die 
erzwungenen Schwingungen durch die zur Fourierreihe (15) gehorige reelle Darstellung 


@(t) = Py + cc? + € sin (Py — 0) cost —€ 2 Day sin (Y — 0) sin T 
ah = sin 2 (~) — 0) cos 2 t + O(e*) . (26) 


1 Fiir Konvergenzuntersuchungen vgl. man die in Fufnote 3 auf Seite 46 genannte Arbeit yon K, L. Stell- 


macher. 
4 
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gegeben. Dab y(z) die Ausgangsgleichung in der angegebenen Naherung tatsachlich erfiillt, kann 
durch Einsetzen von (26) in (14,1) leicht erprobt werden. Da sich die Mittellage @ (und abhangig 
davon auch c??)) mehrdeutig aus (22,1) berechnet, gibt es mehrere, voneinander verschiedene Schwin- 
gungsméglichkeiten. Welcher Schwingungszustand tatsachlich eintritt, kann daher nur eine Sta- 
bilitatsuntersuchung entscheiden. 

b) Periodische Erschiitterung durch Schubkurbelantrieb. Um den Einflub héherer 
Harmonischer kleiner Amplitude in der Erschiitterungsfunktion kennen zu lernen, soll eine Zwangs- 
fihrung durch einen Schubkurbelantrieb betrachtet werden (Abb. 3). Diese spezielle Anordnung 
wurde bereits von K. L. Stellmacher! beim Schwerependel betrachtet und verdient besonderes 
Interesse im Hinblick auf die meist iibliche Priifung der Gerate auf einem zwanglaufig angetrie- 
benen Riitteltisch. Aus der Kinematik des Kurbeltriebes ist die an Stelle von (7) tretende Er- 
schiitterungsfunktion und die zugehérige Beschleunigung bekannt?. Mit dem Stangenverhaltnis 
2 =r/L ergibt sich in diesem Falle die Beschleunigungsfunktion zu 

3 


j3 y) 
q’ =—ril'"(t) =—r [cost Sic (; 47) eos 27— | cos Arey ens (27) 


worin wiederum die dimensionslose Zeitzahlung (13) verwendet wird und auBerdem hohere als 
vierte Potenzen von 7 < 1 unberiicksichtigt bleiben sollen. An die 
Stelle der Amplitude gq, ist hier der Kurbelhalbmesser r getreten, 
so daB& mit der reduzierten Pendellange / fiir die relative Erreger- 
amplitude (10) nun 


eintritt. Nimmt man an, daB die reduzierte Pendellange / des Ge 
rates und die Pleuelstangenlinge L von gleicher GréBenordnung 
sind, so ist der Quotient 

l 


geen y py 


von der GréBenordnung eins, so da sich A als von der Gréfen- 
ordnung € ergibt: 


4 ig 
A= ib =4E 
Abb. 3. Periodische Erschiitterang = _TJamit kénnen die Koeffizienten der Beschleunigungsfunktion (27) 
durch Schubkurbelantrieb, 5 


durch den kleinen Parameter ¢ ausgedriickt werden. Tragt man so- 
dann (27) in (14) ein, so erhalt man in diesem Falle einer allgemeineren Erschiitterungsfunktion 
die Bewegungsgleichung 


g +2 Dapg! + oF (Y — pp) + 04 Cos Y 


a 
+ (eeos 4 + acos2 x +5 at (cos 21 — cos 42) ...)sin (p— 6) ==aiOu 


Unter den gleichen Voraussetzungen und auf die gleiche Weise wie im vorstehenden Abschn. 3a) 
lassen sich auch in diesem Falle die erzwungenen Schwingungen berechnen. Einzelheiten der 
Rechnung bieten kein Interesse, so daB nur das Ergebnis mitgeteilt werden soll. An Stelle von 
(26) ergibt sich bis auf einen Zusatzterm in der zweiten Harmonischen die gleiche Fourierentwick- 
lung: 


P = Po + cP ce? + € sin (py — 5) cos t —e 2 Dag sin (yy — 0) sin t 
cL & sin (@) — 6) + - sin 2 (%@ — ®) cos 2T + O(c). (26,1) 


Hierin bestimmt sich bemerkenswerterweise die Mittellage gy aus der unveranderten Gleichung 
(22,1), so da’ auch bei der Diskussion dieser Mittellagengleichung keine neuen Erkenntnisse ge- 
wonnen werden. Die Korrektur c) zu qm) verschwindet unter den gleichen Voraussetzungen wie 


in (25). Auf die Berechnung von c?) wird verzichtet, da trotz langerer Rechnung keine neuen Er- 
kenntnisse zu erwarten sind. 


* Vel. S. 118 der in FuBnote 3 auf Seite 46 genannten Arbeit. 


* Vel. C. B. Biezeno und R.Grammel, Technische Dynamik, Bd. 2, S. 227, 2. Aufl., Berlin/Géttine 
Heidelberg: Springer 1953. J a Histgstoook Ses toa 
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eee eee roan die Setanta Schwingungen (26) erzwungene Schwingungen 

gen sind, gibt es fiir ein Gerat bei gleicher Erschiitterungsfrequenz 
und gleicher Amplitude mehrere Schwingungsméglichkeiten. Um die wichtige Frage zu ent- 
scheiden, welche der méglichen Schwingungen tatsachlich eintreten, hat man ihre Stabilitat zu 
ers: Aus Griinden der rechnerisch einfachen Durchfithrbarkeit wird die infinitesimale 
ean Da Pe ere zur Entscheidung der Stabilitatsfragen 
en Schwingungen aus (26) und (26,1) wird also aus 
linearen Differentialgleichungen fiir eine kleine Stérung u(r) der abgeanderten Schwingung 9 =@ 4-u 
berechnet. Falls eine kleine, jedoch sonst allgemeine Stérung u(r) sich als mit der Zeit exponentiell 
abklingend ergibt, ist die auf Stabilitat zu priifende Schwingung q stabil in dem Sinne, daB bei 
Kleiner Abanderung der Anfangshedingungen die erzwungene Schwingung @ mit iiberlagerten 
kleinen Stérungen endlicher Amplitude bestehen bleibt und also nicht etwa eine andere der be- 
rechneten méglichen erzwungenen Schwingungen @ eintritt!. 


ae Stérungshewegung u geniigt einer Differentialgleichung, die sich aus der Bewegungs- 
gleichung fiir p durch EKinfiihrung von ¢ = ~ + wu mit anschlieBender Linearisierung in bezu 
u ergibt. Aus (14) erhalt man so die homogene lineare Differentialgleichung 


wu’ +2 Dapw' + [ab —o? sin —e [7 (t) cos (py — 4)] u=0, (28) 


g auf 


in welcher die zu untersuchende Schwingung g(t) mit ihrer bekannten Fourierentwicklung ein- 
zutragen ist und sich alle nichthomogenen Gleichungsterme herausgehoben haben, da g(t) Lésung 
von (14) ist. Um die Stabilitatsgleichung (28) in die auswertbare Form einer Differentialgleichung 
mit periodischen Koeffizienten iiberzufiihren, wird gemaB (15) y(t) nach Potenzen von @(t) um 
die Mittellage qm entwickelt, was auf 


2 


cD 
tT 6 . 5 
u’ +2 Da,pw + si — Xf SIN Py — AF COS Po t + 04 SiN Py 57 + 


—ell''(t) (cos (Go — 0) — sin (% — 6) = — cos (Y — 5)e +-:. | u. ==0 (231) 


fiihrt. Hier hinein ist nun //’'(r) und ¢(t) einzusetzen, was zunachst fiir den Fall rein harmonischer 
Erschiitterung geschehen soll. 


a) Harmonische Erschitterung. In diesem Fall sind //'’ = — cost und g(t) mit ihren 
Fourierentwicklungen (26) als bekannte Zeitfunktionen in (28,1) einzusetzen. Wird fiir die 
Naherung (26) mit einem Fehler der GréBRenordnung ¢? verwendet, so wird die Gleichung fiir u 
eine Naherungsdifferentialgleichung mit einem Fehler der GréBenordnung ¢*. Sie ist eine Hillsche 
Differentialgleichung 


u? -2Doa,u' + V(r) u =0 (28,2) 


mit der periodischen Koeffizienten-Funktion ¥/(t). Zur Vereinfachung dieser linearen Differential- 
gleichung kann man die neue abhangige Veranderliche w(t) durch die Transformation 


Ga ie e FT w(r) (29) 
einfiihren, so daB® sich fiir w(t) die einfachere Differentialgleichung 
wl! + (W(t) — D2 ak) w = 0 (30) 


ergibt. Die in (30) neu aufgetretene Koeffizientenfunktion soll in der fiir die weitere Rechnung 


zweckmaBbigen komplexen Form 


VG Dae — >, dae" (31) 


n=— oo 


geschrieben werden und ihre Koeffizienten d, werden wiederum nach Potenzen von € geordnet. 


Vernachlassigt man vierte und héhere Potenzen von é, so bricht die Fourierreihe (31) nach der 
dritten Harmonischen ab und man erhalt nach einiger Rechnung die Koeffizienten mit den 


1 Vgl. z. B. A. Liapounoff, Ann. Faculté des Sciences de Toulouse, II, 9 (1907) S. 203—469, Neudruck 


Princeton 1949. 
4 * 
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Werten 


eae te 
dy = d®) e® = jak (1 — D*) — of sin q Bias sin? (~)— 0)| €? , 


bo| &! 


d,=d_,=dMe+ do) & cos (Py — | ey aie a COS Yo Sin (Py — 0) — —- sin (% — 9) 


= “s sin 2 (~ — 6) sin (~ — 0) — % cos (Py — 0) sin? (Py — | a 


aR 


Leen aD : 
d, = d®@ 2 + de = ie aq PR (Po — | ey E Paro) sin? (gy — 6)| €°, 


ef f OFF sin? (5 — | &, 


d = d®), roe aL d®), se = = sin? (Po = 0) 


: : i 5 
dda ed = fe A, sin 2 (— d) sin (Yq — 0) — rg €08 (Yo — 9) sin? (Yo — | &, 


dm = 0 fir |m| >3. 

Nach dem Floquetschen Theorem! setzt sich die allgemeine Lésung der linearen Differential- 
gleichungen (31) im vorliegenden allgemeinen Falle aus zwei linear unabhangigen Teillésungen — 
der Form 

Ce el oh. Go eee lr) (32) | 
zusammen. In (32) bezeichnet g@ den charakteristischen Exponenten, wahrend V,(t) und V,(t) 
zu der Koeffizientenfunktion Y — D? xj gleichperiodische Funktionen sind. Mit zwei Integrations- 
konstanten C, und C, laBt sich aus den Teillésungen (32) eine allgemeine Lésung zunachst fiir w — 
und damit auch fiir u zusammensetzen. Uber das Wachstum dieser Stérbewegung u entscheidet 
zufolge der Transformation (29) der charakteristische Exponent g, so das die Berechnung von 0 
die Hauptaufgabe der Stabilitatsuntersuchung ist. Zu dem Zweck geniigt es, eine Teillésung in 
der Form 


cO 


i a ve ei (a + n)t (32,1) 
n=—o 


mit den Fourierkoeffizienten v, zu betrachten. Da die Koeffizientenfunktion durch eine Fourier- — 
reihe mit nach Potenzen von ¢ geordneten Koeffizienten der Differentialgleichung (30) gegeben ist, — 
liegt es nahe, auch die Fourierkoeffizienten v, der Teillésung (32,1) durch solche Entwicklungen — 
schrittweise zu berechnen. Eine derartige rekursible Lésung wird sich als méglich erweisen, wenn 
man auch den charakteristischen Exponenten in gleicher Weise entwickelt. 

Tragt man den Ansatz (32,1) in (30) ein, so wird man fiir die Koeffizienten vy auf das homogene 
lineare Gleichungssystem 


—(N®=b 2 9 Nebo?) oe De ne =0 (N=0,+1,+2,...) (33) 


— — 
gefihrt, worin der gesuchte Exponent ¢ als Parameter enthalten ist. Tragt man die ¢-Entwicklungen — 


y= 09 4+ Me +o@eht--- (n=0,+1,42,...) 4 


n 


= 0 a 0) er 0) Pyle pes a 
in die Gleichungen (33) ein und beachtet den Aufbau der Koeffizienten (31), so lassen sich diese 


Gleichungen rekursibel lésen. Zunachst findet man fiir N = 0 wegen d® — 0 fiir alle n | 
— or = 0, (33,00) 


was durch 9 = 0 und entsprechend dem homogenen Charakter der Gleichungen (33) mit beliebigem 
v gelost werden kann. Diese Lisung entspricht durchaus bekannten Ergebnissen iiber die Mathieu- 
sche Differentialgleichung, mit welcher (30) fiir kleine Parameterwerte eng verwandt ist. Es wird — 
v) — 1 gewahlt. Fiir N A 0 ergibt sich weiter aus (33) mit 62 30 | 


— NN? — 0. a= 0) =U (N ae ee), (33, N, 0) 


Weiter liefert fiir Potenzen ¢’ und N = 0 (33) zunachst nur eine leere Aussage. Dagegen folgt 
fiir N ~ 0 aus (33) 


— NP of— 2% of + ¥ (dQ_, of + dQ_, o) = 0. (33, N, 1) 


‘ Vel. z. B. J. Horn, Gewohnliche Differentialgleichungen S. 96, 5. Aufl. Berlin 1948. 
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(0) ie ne 2 A 
: lee 0 fiir N a i kommt das unbekannte 9 in (33, N, 1) nicht vor. Da es auBerdem nach (31) 
eine Bestandteile d‘Y und nur die Anteile dQ, d®), gibt, ist die Ausrechnung einfach und fihrt 


mit v® — 1 auf 
. v(t) = d{), ve), = d®, , w= 0 fir N=+2,4+3,.. 
Der Vergleich der Potenzen ¢? in (33) fiir N — 0 fiihrt auf 


SO By SIU, re) dO ru) ee 0; 


: : 
was mit den berechneten Werten von vo und v® = 1 zur Bestimmung von o dient. Es ergibt 
sich die wichtige Beziehung ; 


gO — dl, of) + al al, 4d) = a + 2 ah al, ey 
woraus sich spater die Stabilitatsfragen entscheiden lassen. 
Zunachst soll noch das in der Entwicklung von 9 folgende Glied 9® berechnet werden. Zu dem 
Zweck muB vorbereitend zunachst (33) auch fiir N+ 0 und & ausgewertet werden. Aus 
— N?rQ@— 2 No r— py + y (aW_, v® + dQ@_,, ) =0 


erhalt man fiir die sogleich benétigten Koeffizientenanteile v® und v®, 


v?) = —2 eM rf), po ae Oy on 


Damit kann @® aus (33) fiir N = 0 durch Betrachtung der Potenzen ¢* berechnet werden. Die 


Gleichung 
ee Oo ay (dave) dh, 0) 4 d®) se) 0 


reduziert sich bei Beachtung der Entwicklungen (31) und mit den Lésungen v‘!) und v® auf 

2 of gf) — Sal, a2) = af a2 + a, of 
Hieraus folgt 

2 0!) o® = 2 of (d®, v®) — dM Pe Nea Ue 
weil vo) = v™, und d® = d®, ist. Da 0 40 berechnet wurde, ist also 0° = 0. Mithin gilt 
fiir den charakteristischen Exponenten nach (34) 
@ =e + (6%) = yd?) + 2 dO, de + 0), 

was spater ausfiihrlich diskutiert werden wird. 

Zunachst soll jedoch gezeigt werden, da® auch fiir den behandelten Fall der allgemeineren 
Erschiitterung diese Beziehung gilt. 

b) Periodische Erschiitterung durch Schubkurbelantrieb. Auf Grund der Erschiit- 
terungsfunktion (27) stellt sich die Zwangsschwingung (26,1) ein, fiir die man wiederum die Stabili- 
tatsgleichung (28,2), jedoch mit einer etwas anderen Koeffizientenfunktion Y erhalt. _Beschrankt 
man sich auf quadratische Potenzen von ¢, so ist lediglich bei d, gegeniiber (31) die Anderung in 


dA de doe > cos (y — 0) — s sin? (~) — 0)| € 


zu beachten. Da die Gleichungen (33,00), (33,N,0) und (33,N,1) davon unbeeinfluBt bleiben, gilt 
der Ausdruck (34) fiir den charakteristischen Exponenten auch jetzt noch. Eine weitergehende 
Rechnung bestatigt ferner auch hier 0°) = 0. 

Die im Vorstehenden geleistete Berechnung des charakteristischen Exponenten gestattet die 
Beantwortung der Stabilitatsfrage in dem hier bendtigten Umfange sowohl fiir rein harmonische 
Erschiitterung als auch fiir den behandelten Fall der allgemeineren Erschiitterung durch einen 
Schubkurbelantrieb. 

Zufolge von (29) und (32) kommt es zur Beurteilung des Zeitverhaltens einer Stérung wu auf 
den Stabilitatsindex 

I = Re (— Dag + 70) (35) 
an. Nur wenn I < 0 wird, ist die zugehbrige Schwinguug p(t) stabil im dem hier verwendeten 
Sinne. In (35) hat 9 den durch (34) gegebenen Wert. Um die Abhingigksit von den eigentlichen 
mechanischen Gri8en des Problems zu erkennen, wird 0? mit den urspriinglichen Bezeichnungen 
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in die stets reelle Form 

02 re oll e® = ab (1 — D®) —a2sin yy + 7 008 2 (% — 6) = M'(y) — D? ak (36) 
gebracht. Falls M' — D? a} > 0 gilt, ist fir das Gergt und die gewahlten Werte von Erreger- 
frequenz und Amplitude [<< 0 und die Stabilitat der Schwingung um die Mittellage qm) daher 
gesichert. Falls jedoch M'— D? a < 0 ist, wird der zweite Bestandteil in (35) reell und kann 
dann Instabilitat anzeigen. Dieser Instabilitatsfall tritt ein, wenn aus (35) mit (36) 


I = Re(—Doy + i /M'(gy) — D¥ak) 2 0 
folgt und also die beiden Ungleichungen 


M'(y) — D?ak <0, — Dar + JD? axb— Mg) = 0 


gleichzeitig erfiillt sind. Dies ist gegeben, wenn nur M'(p)) < 0 ist. 
Mithin herrscht keine Stabilitat mehr, wenn 


a2 ie 
ag 88 Po + = 00s 2 (Y — 9) < 0 (37) 
ist, worin ¢ wieder den Intensitatsparameter von (24) bezeichnet. 

Uberraschenderweise ist fiir die Entscheidung iiber die Stabilitat in dieser Naherung die Kennt- 
nis der Dampfung nicht erforderlich. Die Gleichung (37) gilt, wie auseinandergesetzt, auch fir 
den betrachteten Fall allgemeinerer Erschiitterung. Eine etwaige Abweichung des Federaufzugs- 
winkels und eine dadurch bedingte statische Ruhelage des Pendels auBerhalb der Horizontalen 
geht iiber die Mittellagengleichung (22,1) in qp ein und ist infolgedessen auch von starkstem Ein- 
fluB fiir die Stabilitatsaussage (37). Zahlenbeispiele werden die Bedeutung dieses Nullfehlers fiir 
die Stabilitat der Schwingungen deutlich machen. 

Ohne auf Zahlenrechnungen zuriickzugreifen, laBt sich die Stabilitat eines MeBgerates mit 
vertikaler Erschiitterung (0 = 90°) und mit fehlerfreiem Federaufzugswinkel y, nach (4) unter- 
suchen. Gemaf (22,1) mit 1 = 0 besteht hier die Méglichkeit der erzwungenen Schwingungen um 


g(t) den Sollwert der statischen Ruhelage gy =0. Dann besagt (37) 


oS —_—“" wegen 2 (~,— 0) = a einfach, daf Instabilitat eintritt, wenh 
d=350° nz, 7 
Zee 
a : (38) 


ist. Ausgedriickt in den urspriinglichen mechanischen GréBen ergibt 
diese Ungleichung zufolge der Abkiirzungen (8), (9), (10) und (24), 
daB fir Instabilitat 


Pp uy = V2 lor (39) 


sein muB. Demnach mu die Geschwindigkeitsamplitude der har- 


monischen Erschitterung gréBer sein als j2 mal dem Produkt aus 
reduzierter Pendellange und der Eigenfrequenz des Gerates. Diese 


ig Bedingung fiir Instabilitat der Schwingungen um die statische 
Abb. 4, Hecuersal erpebslients Sehwere- Gleichgewichtslage ist auffalligerweise formal genau die gleiche, 
pene el zum erg ereh, 


wie man sie von einem verwandten Problem des erschiitterten 
Schwerependels her kennt. Wird namlich ein Schwerependel rechtwinklig zu seiner statischen 
Ruhelage erschiittert (Abb. 4), so kann die Instabilitatsgrenze auch bei diesem andersartigen 
Problem aus (39) berechnet werden, wenn man @, durch die beim Schwerependel allein vorhandene 


Kigenfrequenz @,) = yell ersetzt. Man kann dann (39) auch in die Formel von K. Klotter! 
UE Ore eg 


iiberfiihren. Die Verwandtschaft der Probleme des erschiitterten Schwerependels und des er- 
schiitterten MeSgerates ist also noch zu erkennen. Allerdings sind die Auswanderungslagen des 
MeBgerates wesentlich durch die Eigenschaften der Riickstellfeder bedingt. Da das Schwerependel 
keine Federriickstellung besitzt, sondern allein durch die Schwere zuriickgefiihrt wird, hat die 
Verwendbarkeit der Untersuchungen des Pendels als Modell fiir Auswanderungserscheinungen von 
Mef- und Anzeigegeraten beziiglich der zahlenmaBigen Ergebnisse enge Grenzen. 


* Vel. S. 379 der in FuBnote 2 auf Seite 43 genannten zusammenfassenden Darstellung. 
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Die vorstehenden Rechnungen beschranken sich auf die Berechnune 
finitesimale Stérungen der zu priifenden Schwingungen. Tatsachlich ae 
Stabilitat gegen endliche Stérungen den praktischen Gegebenheiten besser angepaBbt, da man 
beispielsweise sogar beim stehenden MeBgerat schon durch eine hinreichend groBe Anfangsaus- 
lenkung das MeBpendel iiber eine instabile Federgleichgewichtslage hinausfiihren kann, so daB eine 
andere und stabile Gleichgewichtslage angenommen wird (vgl. Abb. 5, statische Auslenkung 
groBer Po sau, = 41,4° fihrt auf stabile Gleichgewichtslage Postar. = 242,4°). Doch wiirde ahs 
derartige Stabilitatsuntersuchung im GroBen mathematisch so weitlaufig werden, da die Be- 
schrankung auf die Untersuchung der infinitesimalen Stabilitat geboten erschien. Bei einer experi- 
mentellen Priifung der Stabilitatstheorie hat man entsprechend vorsichtig zu arbeiten, um gréBere 
Stéramplituden zu vermeiden. ‘ 


der Stabilitat gegen in- 
e eine Untersuchung der 


3. Numerische Ergebnisse’. Die allgemein abgeleiteten Ergebnisse sollen nun durch spezielle 
Zahlenbeispiele illustriert werden. Zu dem Zweck werden die Auswanderungsmittellagen qy ab- 
hangig von der jeweiligen Erschiitterungsintensitat ( berechnet. In dieser Rechnung wird stets 
streng harmonische Erschiitterung (//(t) = cost) in der Lotrichtung (6 = 90°) vorausgesetzt, 
weil diese Erschiitterungsbelastung die in der Praxis haufigste sein diirfte und auBerdem dieser 
Fall zahlenmaBig am einfachsten durchzufiihren ist. Die jeweiligen fiir das MeBgerat kennzeich- 
nenden Geratedaten sind die Kigenfrequenzen «, und «9 (aus (8) und (9), die durch den Quotienten 
@/@; mit verschiedenen Werten in die Rechnung aufgenommen werden. Als fiir die Zahlenergeb- 
nisse sehr wesentlich wird sich, wie schon erwabnt, ein etwaiger Nullfehler durch fehlerhaften 
Federaufzug erweisen. Die Rechnung darf also nicht auf das fehlerfreie Gerat (A = 0 nach (23)) 
beschrankt bleiben, sondern muf auch Werte A = 0 fir nicht streng horizontale statische Ruhe- 
lage beriicksichtigen. Da der Mittellagenwinkel gy die erzwungenen Schwingungen kennzeichnet, 
sind nach der Berechnung von g als mehrdeutige Funktion des Intensitatsparameters ¢ alle Schwin- 
gungsmoglichkeiten des Gerates bekannt. Eine abschljeBende Rechnung hat dann zu entscheiden, 
ob die jeweiligen Schwingungen stabil sind oder als instabil keine physikalische Realitat besitzen. 

Zur Bestimmung der Mittellagen q, sind die Wurzeln der transzendenten Gleichung (22,1) fiir 
von Null an wachsenden Intensitatsparameter ¢ aufzusuchen. Zu dem Zweck geniigt im allgemeinen 
fiir die angestrebte begrenzte Genauigkeit der Zahlenrechnung eine grafische Auswertung, wobei 
Qo als Schnittabszisse der beiden Kurven F, = €/4 sin 2 y und Fy = (yy — A) alo? + cos yy — 1 
bestimmt werden. Die Variation des Parameters ¢ verlangt also lediglich die Zeichnung von sin- 
Kurven verschiedener Amplituden, wahrend ein Nullfehler nur eine Parallelverschiebung der 
Kurven F, erforderlich macht. Fir jedes Gerat mit anderen Eigenfrequenzen hat allerdings 
of/x2 = w/w? einen anderen Wert, so daB die Kurve F, neu zu zeichnen ist. Da die Stabilitats- 
untersuchung fiir ein Wertepaar g, ¢ auf die Frage nach dem Vorzeichenwechsel der linken Seite 
von Gleichung (37) hinauslauft, kann man nach den speziellen Wertepaaren , ¢ fragen, fiir welche 


gerade die Stabilitatsgrenze erreicht wird. Diese Werte qo, ¢ berechnen sich aus 


y 2 = 
s — sin Y — ey cos 2g, = 0 (40) 
0 

und kénnen als Kurven zur Beschreibung der Stabilitatsgrenzen in die Diagramme der Mittel- 
lagen ¢,() eingetragen werden. Durch (40) ist also die Stabilitatsentscheidung wieder auf die 
Lésung einer transzendenten Gleichung zuriickgefiihrt. Die Auswertung geschah wiederum gra- 
fisch durch Aufsuchen der Schnittpunkte der Kurve G, = a%/aj;— sing) und G, = C/2 cos 2 gy. 
Fiir Kurventeile, bei welchen die grafische Auswertung keine zufriedenstellende Genauigkeit ergab, 
wurden einzelne Werte von ¢ fiir vorgegebene Werte von gy) unmittelbar aus den Gleichungen 
(22,1) und (40) berechnet. 

In Abb. 5 sind die Mittellagenwinkel gy, abhangig von der Erschiitterungsintensitat ¢ fiir einen 
bewahrten Schwingweggeber mit w}/w? = 0,346 aufgetragen. AuBer den Kurven fiir das febler- 
freie Gerat (A = 0) sind auch die Werte ¢(C) fiir das gleiche Gerat bei etwa vorhandenen Null- 
fehlern (A = —15°, —10°, —5°, +5°, +10°, 4 15°) eingetragen. Die eingezeichneten Stabilitats- 
grenzkurven teilen die (qo, ¢)-Ebene in Gebiete auf, in denen Wertepaare liegen, die entweder zu 
stabilen oder zu instabilen Schwingungen gehéren. 

Man entnimmt dem Diagramm zunachst fiir ein fehlerfreies Gerat (A = 0) die wichtige Tat- 
sache, daB zu einem bestimmten Intensitatswert der Erschiitterung ¢ = qj w?/g 1 stets mehrere 


1 Bei den numerischen Auswertungen haben mich die Herren F’. Leiss und R. Buttmi wirksam unter- 
stiitzt, was durch eine Zuwendung der Karlsruher Hochschulvereinigung e. V. in dankenswerter Weise ermog- 


licht wurde, 
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Werte der Mittellage gy gehéren. Es besteht also fiir jede beliebige Amplitude und Frequenz der 
Erschiitterung ein -mehrdeutiger Schwingungszustand. Allerdings sind in wirklichen Geradten im 
allgemeinen nur solche Schwingungen mdglich, deren Mittellagen q in unmittelbarer Nachbar- 
schaft der Sollage gy = 0 liegen, da meist Anschlage die Bewegungsméglichkeit des Schwingers 
begrenzen, so da® dieser bei groBeren Betragen von @p anliegt, ohne noch schwingen zu kénnen. 
Betrachtet man das 

ea fehlerfreie Gerat bei aus 
270° der Ruhe anwachsender 
Erschiitterung, so findet 
man den Schwinger zu- 
nachst in der statischen 
Ruhelage gy =9. Er be- 
ginnt sodann bei Er- 
°F Se) schiitterung mit kleiner 
s Za ee Intensitat € mit der glei- 
fi chen Lage q als Mittel- 
Va lage zuschwingen. In der- 
4 folgenden Betrachtung 
wird vorausgesetzt, daB 
sich diese Schwingungen 
ungehindert ausbilden 
kénnen, wozu also et- 
waige Anschlagbegren- 
zungen entfernt werden 
miissen. Steigert man 


120° jee ad 
a W/, 
; Y die Intensitat € weiter 


90°} Y stabil und iiber den aus (40) 


r % folgenden Wert 
L es 6 = 2 af/a? = 0,692 
60° ae ‘y, 


780° 


hinaus, so ist zwar @) = 0 
immer noch Lésung von 
(22,1), doch hat man die 
Stabilitatsgrenze in ein 
instabiles Gebiet hinein 
tiberschritten. Die zu- 
gehérigen Schwingungen 
um die gewiinschte Mit- 
tellage ~ = 0 sind also 
instabil und der Schwin- 
ger kann auswandern. 
Eine kleine Stérung wird 
ihn im allgemeinen zu 
Schwingungen um nega- 
tive, d.h. unter der Hori- 
zontalen liegenden Mit- 


instabil 


oO 
Abb. 5. Auswanderungs-Mittellagen ~, abhangig von der Erschiitterungsintensitat € mit Feder- tellagen Po veranlassen. 
aufzugsfehlern 4 = 0, +5°, +10°, +15° als Kurvenparametern fiir harmonische Erschiitterung Bei wachsender Erschiit- 


in Lotrichtung. . Liars : 
. terungsintensitat wird 


die zugehérige Mittellage gy) mehr und mehr aus dem Wert q) = 0 auswandern und schlieflich 
bei tiber alle Grenzen anwachsender Intensitat den Wert g) = — 90° annehmen. Bei der voraus- 
gesetzten vollen Bewegungsfreiheit des Schwingers kénnte im Instabilitatsfall die Auswanderung 
aus der Mittellage ~) =0 auch durch die Uberkopflage g = + 90° in Schwingungen mit Mittellagen 
in der Nahe der stabilen Gleichgewichtslage @p gia, = 242.4° hinein erfolgen. Bei groBen Intensi- 
taten nahert sich die Mittellage dieser Schwingung dem Wert +270°, so daf sich der Schwinger 
schlieBlich wiederum der gleichen Vertikalen (— 90°) nahert. Dieses Verhalten ist durchaus ver- 
standlich, denn bei unbegrenzt wachsender Intensitat wird der Schwinger auf jeden Fall vollends 
in die vertikale Erschiitterungsrichtung gezogen. Bei grofen Intensitaten (etwa ¢ >3) weist 
die Abb. 5 tatsachlich auch stabile Schwingungen in der Uberkopflage aus. Bei unbegrenzt 


\ 
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wachsender Erschiitterungsintensitat nimmt die Mittellage @  schlieBlich den Wert y, = 90° 
an, so daB das Hineinwandern in die Erschiitterungsrichtung bei hinreichend groBen Tuieeiaien 
also auf allen Wegen méglich ist. 

Fiir die Praxis sind nur kleine Erschiitterungsintensitaten ¢ und wegen der Ausschlagbegrenzung 
auch nur kleine Mittellagen y von wesentlichem Interesse. Richtet man in der Abb.5 die ee 
merksamkeit auf den dann verbleibenden kleinen Ausschnitt in der Nahe des Ursprungs, so stellt 
man hier einen besonders grofen Einflu® etwaiger Fehler beim Aufzug der Feder fest. Fin Auf- 
zugsfehler gemaf (23) hat fiir das nicht erschiitterte Gerat eine statische Ruhelage zur Folge 
in welcher der Schwinger nicht genau in der Horizontalen liegt, sondern davon um einen Winkel 
Post, abweicht. Dieser Winkel berechnet sich aus (22,0) fiir ¢ = 0 wiederum als Lésung einer 
transzendenten Gleichung. Die Ergebnisse dieser Rechnung fiir verschiedene aired ale A 
sind in der Tabelle 1 angegeben. Ferner sind in die Tabelle auch die Werte Po stat, fur diejenigen 
statischen Gleichgewichtslagen aufgenommen, die nicht in der Nahe der Sollage nee 0 liegen. Fiir 
das betrachtete Gerat gibt es im allgemeinen drei statische Gleichgewichtslagen, von denen jedoch 
die im ersten Quadranten oberhalb @p 5,4, = 20,3° liegenden nicht stabil sind. Diese statische In- 
stabilitat erkennt man mit Hilfe des Potentials (3) daran, daB fiir die in Betracht kommenden 
Werte P = Post. 1M ersten Quadranten CV/0y?< 0 wird. Bemerkenswerterweise gibt es fir 
eine um ein Weniges zu weit aufgezogenen Riickstellfeder (4 > + 10°) bereits keine statische 
Gleichgewichtslage in der Nahe der Sollage yy = 0 mehr, wahrend fiir 1 = 10° gerade beide Gleich- 
gewichtslagen im ersten Quadranten zusammenfallen. Als einzige statische Ruhelagen verbleiben 
dann nur praktisch uninteressante Lagen im dritten Quadranten. 


Tabelle 1 
A —15° =i se eee 102", Ise 
| | 
| | | | 
Sie Y= 82 ase ‘oe 5908 | 20.2" a 
> Stas. 54,5° 50,7° 46,4° 41,4° 35,0° 20,2° a 
238,1° 239,6° 241,0° 242,4° 243,8° 245,2° 246,6° 


Selbst der EinfluB eines nur kleinen Nullfehlers auf die Auswanderungserscheinung ist von 
iiberraschender GréBe. Betrachtet man z. B. ein Gerait mit einer nur geringfiigig zu weit auf- 
gezogenen Riickstellfeder, so wird dadurch der stabile Arbeitsbereich des Geriates so stark ver- 
kleinert, daB seine praktische Brauchbarkeit in Frage gestellt ist. Z. B. entnimmt man der Abb. 5 
fiir einen Federaufzugsfehler A = + 5°, fiir welchen also der Schwinger in der statischen Ruhelage 
um nur 5,9° aus der Horizontalen herausgehoben ist, bereits ein starkes Absinken der Stabilitats- 
grenze von etwa £ = 0,7 auf etwa € = 0,2. Fiir gréBere Aufzugsfehler im gleichen Sinne ist das 
Gerat itiberhaupt nicht mehr zu stabilen Schwingungen um die statische Sollage fahig. Ganz andere 
Konsequenzen hat ein zu schwacher Aufzug der Riickstellfeder. Wie man aus der Abb. 5 ersehen 
kann, hebt auch die kleinste negative Winkelabweichung A die Instabilitat in dem Sinne vdllig auf, 
daB nun ein stetiges Anwachsen des Betrages der Mittellage zu negativen Werten von gp hin ein- 
tritt. Das bedeutet fiir den Schwinger bei wachsender Erschiitterungsintensitat eine allmahliche 
Auswanderung der Mittellage gy, auf den Wert gy = —-90° zu. Ein solches Verhalten ist zwar fiir 
Anzeigegerate durchaus unerwiinscht, doch liefert ein Schwingweggeber in diesem Falle immer 
noch auswertbare Anzeigen. Falls man also Toleranzen fiir den Bau von Schwingweggebern hin- 
sichtlich der Horizontallage des Gebers in der Ruhelage vorgeben mu8, sind bei vertikaler Er- 
schiitterung Lagen des Schwingers unter der Horizontalen mit zugehérigem zu schwachem Feder- 
aufzug unbedingt zu bevorzugen. Fiir die iibrigen stabilen Mittellagen gy haben die Aufzugsfehler 
nur den aus der Abb. 5 ersichtlichen geringen EinfluB im Sinne einer kleinen Verschiebung der 
Zahlenwerte gp, die fiir das fehlerfreie Gerat berechnet werden. 

Die vorstehenden Ergebnisse gelten nur fiir ein spezielles, durch das Frequenzverhaltnis 
wr|o? = c/M g s = 0,346 gekennzeichnetes Gerat. 

In einer weiteren Rechnung wurden daher die Gerate mit verschiedenen Daten in den Grenzen 
0 < wy/w? < 2 untersucht. Dabei wurden nullfehlerfreie Gerate vorausgesetzt und auch die Rech- 
nung auf den vorzugsweise interessierenden Bereich in der Nahe der Sollage yy = 0 beschrankt. 
Die Ergebnisse der Berechnung von gy und der zugehorigen Stabilitatsuntersuchung sind in Abb. 6 
aufgetragen. In dieser Abbildung sind diejenigen Kurvenaste (6). die zu einer instabilen Schwin- 
gung gehdren und daher keine physikalische Bedeutung haben, gestrichelt eingezeichnet. Man ent- 
nimmt der Darstellung zunachst die einleuchtende Tatsache, daB mit klemer werdendem w/w? 
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auch die Stabilitatsgrenze fiir die Schwingungen um die Sollage absinkt. Aus Gleichung (40) ist 
der Wert der kritischen Erschiitterungsintensitat zahlenmaBig bekannt: ¢ = 2 wp/@>. Demnach 
sind gerade Schwingweggeber mit der bei diesen Geraten erwiinschten tiefen Eigenfrequenz 


Or = yelJ besonders anfallig fiir Auswanderungserscheinungen und man muB sich daher bemihen, 


die Frequenz @,) = Veil durch eine groBe reduzierte Pendellange / von gleicher GroBenordnung 
wie wp klein zu halten. Aus konstruktiven Griinden kommt man hier allerdings bald an eine Grenze. 
Gerate mit hohen Eigenfrequenzen kénnen in der Praxis wohl kaum auswandern, da die bei ihnen 
erforderlichen groBen Erschiitterungsintensitaten selten erreicht werden. 

Der Verlauf der y(¢)-Kurven fiir wachsende Werte wy/w; ist durchaus plausibel. Betrachtet 
man zunachst Gerate mit sehr weicher Riickstellfeder, so ist w/w? klein und man erhalt auBer der 
stets vorhandenen Kurve gy = 0 eine Kurve vom Typ der fiir p/w? = 0,346 gezeichneten. Sie 
beginnt in der instabilen Ruhelage @ s,,, = 41.4° und trennt dann auf der ¢-Achse den stabilen 


> 


% 


Abb. 6. Auswanderungs-Mittellagen gy) in der Nahe der Sollage g~, = O fiir verschiedene Gerate. Kurvenparameter ist wher. 


gestrichelte Kurventeile gehéren instabilen Schwingungen zu. 


vom instabilen Bereich. Dabei gehért die @-Kurve von diesem Durchgang an selbst zu stabilen 
Schwingungen und mit wachsenden Intensitaten nahert sich gw seinem Grenzwert ~ = — 90° 
fiir unendlich groBe Intensitaten. 

Die strichpunktiert eingetragene Kurve fiir w}/w?= 0 hat keine wirkliche Bedeutung, sondern ist 
als gedachter Grenzfall fiir eine Feder mit verschwindender Federsteife aufzufassen. Denn fiir 
iiberhaupt fehlende Feder verschwindet auBer m, auch y, und man hat g, nicht aus der Gleichung 
(22, 1), sondern aus der noch nicht umgeformten Gleichung (22) zu berechnen. Man erhalt dann 
richtig @y = + 90° als mégliche Ruhelagen und wird auf die Theorie des erschiitterten Schwere- 
pendels zuriickgefiihrt. 

Verfolgt man in Abb. 6 die Kurven fiir wachsende Parameter w}/@?, so kommt man schlieBlich 
fiir wF/w? = 2/7 = 0,036... auf eine Kurve, die als instabile statische Ruhelage gerade die Uber- 
kopflage gy = 90° besitzt. Diese Kurve ist insofern die letzte Kurve des beschriebenen Types, 
als fiir gréBere Parameterwerte nur noch zwei aus der -Achse verzweigte Kurvenaste der g-Kurven 
entstehen. Der zu kleineren Werten von g, gehérende Zweig nahert sich mit wachsender Intensitat 
wieder dem Wert y, = — 90°. Er gehdrt stets zu einer stabilen Schwingung. Der obere Ast der 
Po" Kurve lauft auf den Wert gy) = 90° zu und ist bei groBen Intensitaten ebenfalls einer stabilen 
Schwingung zugeordnet. Dieses Ergebnis zeigt die Stabilisierbarkeit des Schwingers in der Er- 
schiitterungsrichtung bei groBen Intensitaten an und ist auch aus der Pendeltheorie gelaufig. 

Die in Abb. 5 und 6 zusammengestellten Ergebnisse der Zahlenrechnung diirften die wesent- 


lichen Aussagen der vorliegenden Theorie illustrieren, so da anschlieBend nur noch iiber einige 
experimentelle Ergebnisse berichtet werden soll. 
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6. Experimentelle Ergebnisse. Um einige wesentliche Ergebnisse der vorstehenden Theorie 
zu bestatigen, wurden von G. Benz und H. Heidenhain einige Versuchsreihen mit einem Modell 
durchgefiihrt. Ziel dieser Versuche war vor allem, die bemerkenswerte Abhangigkeit der Aus- 
wanderungserscheinungen von dem Federaufzugswinkel auch in einem Experiment sichtbar zu 
machen. Dabei sollten nur die Ergebnisse der Theorie und der Zahlenrechnuns nachgepriift werden 
nicht aber die Voraussetzung der Vernachlassigbarkeit der Hedenichwingwnped fiir die Nnteatante 
derungserscheinungen. ‘ 

Zur Durchfiihrung der Versuche wurde ein kleines Horizontalpendel sorgfaltig in Kugellagern 
gelagert und mit einer verstellbaren Torsionsdrahtriickstelluug in die horizontale Daceroetanet. 
Der Torsionsdraht wurde in Richtung quer zur Erschiitterung so steif gewahlt, daB er nicht zu 
nennenswerten Querschwingungen angeregt 
wurde. Die ganze Anordnung wurde auf 
einem Riitteltisch fest verschraubt und 
konnte in vertikaler Richtung erschiittert 
werden (Abb.7). Dabei war die Amplitude 
durch die Exzenterstellung und die ge- 
wahlte Untersetzung des Antriebes durch 
einen Kurbeltrieb fest gegeben. Hingegen 
konnte durch den zwangslaufigen Antrieb 
des Tisches die Erschiitterungsfrequenz 
unabhangig von der Amplitude in weiten 
Grenzen tiber die Drehzahl des Antriebs- 
motors geregelt werden. Die fiir die Messun- 
gen erforderliche hohe Frequenzkonstanz 
und der Ubergang auf jeweils nahe benach- 
barte Frequenzen machte allerdings die Ver- 
wendung eines gittergesteuerten Gleich- 
richters zur Speisung des Gleichstrom- 
Antriebsmotors erforderlich. Da sich kleine 
Abweichungen vom rein harmonischen Ver- 
lauf der Erschiitterungsbeschleunigung in 
der Theorie als nicht wesentlich erwiesen  / 
haben, war die Verwendung des Rittel- (“7777777 UL 
tisches in dieser Hinsicht unbedenklich. Abb. 7. Versuchsanordnung (schematisch). 

Auf den Einbau eines besonderen Dampfers 
in das Modell wurde verzichtet, weil sich die Dampfungskrafte in der Theorie als nicht wesentlich 


Bibs Horizontale 


ee 


elek tron. Lihlwerk 


gezeigt hatten. 

Die Messungen wurden mit fester Erschiitterungsamplitude (qg = 5 mm) und mit von Null an 
wachsenden Frequenzen ausgefiihrt. Dabei wurden die Frequenzen durch Zahlung mit einem 
elektronischen Zahlwerk wahrend eines Zeitintervalles von 30sec. Dauer bestimmt. Die beob- 
achteten Mittellagenwinkel (d. h. gy) + ¢)) wurden jeweils aus einer photographischen Aufnahme 
mit 10 sec. Belichtungsdauer entnommen. Zu dem Zweck wurde die Horizontale durch eine Wasser- 
waage kenntlich gemacht und mit photographiert. Nach Entwicklung des Films lieBen sich dann 
unter dem VergréBerungsgerat die Abstande des UmriBbildes des Schwingers von der Horizontalen 
messen und daraus der Mittellagenwinkel geniigend genau berechnen. 

Um das Frequenzverhaltnis ,/m, fiir das Schwingermodell zu ermitteln, wurden die Eigen- 
frequenzen @, und @, des Schwingers einzeln mit Hilfe von Ausschwingversuchen bestimmt. 
Zunachst wurde das Modell als Horizontalpendel durch eine Anfangsauslenkung aus der horizon - 
talen Gleichgewichtslage zu kleinen Schwingungen veranlaBt und daraus die Kigenfrequenz w », 
gefunden. Zur Bestimmung der Eigenfrequenz jw, dagegen wurde die Federriickstellung vollstandig 
abgeschaltet, so da® der Schwinger als reines Schwerependel schwingen konnte. Um bei diesen 
Ausschwingversuchen die Zeitabhangigkeit des Winkelausschlages registrieren zu kénnen, wurde 
der Pendelkérper in den Strahlengang einer beleuchteten Photozelle eingebracht, so dafs die Ab- 
deckung der Zelle als elektrische Spannungsschwankung auf dem Bildschirm eines Kathodenstrahl- 
oszillographen sichtbar gemacht werden konnte. Der Ausschwingvorgang lief sich so auf dem 
Bildschirm zugleich mit einer Zeitmarke photographieren und im Hinblick auf die gesuchten Figen- 
frequenzen @, und @, und die reduzierte Pendellange auswerten. Zur Herabsetzung der Lager- 
reibung wurde bei diesen Ausschwingversuchen zwar das Modell vorsichtig erschittert, aber dennoch 
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konnten die Frequenzen w, und @, nicht sehr genau bestimmt werden. Diese Tatsache hat jedoch 
keine Einschrankung der Versuchsergebnisse zur Folge, da der gefundene Wert w/w; = 0,353 
zweifelsfrei sicherstellt, daB zu den Versuchen Kurven vom Typ der berechneten und in Abb. 5 - 
wiedergegebenen gehéren. Die reduzierte Pendellange ergab sich aus @, zu 1 = 10,6 cm. 

In Abb. 8 sind die gemessenen Werte des Mittellagenwinkels (d. h. gp + Cy) abhangig von der 
Erschiitterungsintensitat ¢ = q? ?/[? w% fiir fiinf verschiedene kleine statische Vorauslenkungen 
aufgetragen. Darin kennzeichnen die Pfeile bei positiven Werten von qo diejenigen Intensitaten, 
bei denen mit Sicherheit keine stationaren Schwingungen zu finden waren. Jedoch waren bereits 
bei etwas kleineren Intensitaten keine reproduzierbaren MeBwerte fiir gy mehr zu erhalten, da bei 

diesen Werten von € die Stabilitat 
der Schwingungen schon verschwin- 
dend klein war. Hingegen erreicht 
man auf den zu negativen MeBwerten 
von @ gehérenden Kurven keine 
Stabilitatsgrenze. Jedoch konnten 
créBere Erschiitterungsintensitaten 
als etwal = 0,5 mit dem verwendeten 
Schiitteltisch nicht erreicht werden. — 
Die MeBwerte der Abb. 8 lassen den 
gleichen Verlauf und die gleiche Ab- 
hangigkeit von einem fehlerhaften 
Federaufzug erkennen, wie die in 
Abb.5 (ohne die Korrekturen cy) auf- 
getragenen, berechneten Mittellagen 
Qo. Die Betrage der Korrekturen cy, 
waren nach der Rechnung in der 
GréBenordnung ¢? = 0,0022 ~ 0,07° 
zu erwarten und konnten daher als 
praktisch bedeutungslos und auber- 
halb der MeBgenauigkeit liegend un- 
beachtet bleiben. Die Ubereinstim- 
Abb. 8. Gemessene Auswanderungs-Mittellagen (+ ¢) in Abhiingigkeit von der mung von Rechnung und Messung ist 

Erschiitterungsintensitat ¢ fiir verschiedene Federaufzugsfehler 4 bei vertikaler sls , 
Eevehanvenuig: damit in dem hier angestrebten Um- 

fange gezeigt. 

7. Zusammenfassung. Drehschwingungsfahige Systeme wie Schwingweggeber oder Drehspul- 
systeme in elektrischen MeBgeraten kénnen unter der Wirkung von Erschiitterungen sehr hoher 
Frequenz erzwungene Schwingungen um andere Mittellagen als ihre statischen Sollagen aus- 
fiihren. Dieses Phanomen ist seit langem als ,,Auswanderungserscheinung“* bekannt, und man hat 
das Verhalten des erschiitterten Schwerependels haufig zum Verstandnis dieser Erscheinung heran- 
gezogen. Da aber das Schwerependel als Riickstellkraft allein die Schwerkraft besitzt, fiir den 
Schwingweggeber die Riickstellung jedoch durch eine Feder erfolgt und die Schwerkraft in anderem 
Sinne in das Kraftespiel eingreift, wurde zur Berechnung der Auswanderungserscheinungen eine 
erweiterte Theorie entwickelt. Die notwendig nichtlineare Rechnung erfaBt eine Riickstellung 
durch eine Spiralfeder, wobei jedoch in der Rechnung von den Schwingungseigenschaften dieses 
Federkontinuums abgesehen wird. Die Untersuchung des dann verbleibenden nichtlinearen 
Schwingungsproblemes mit einem Freiheitsgrad und die anschliefende Stabilitatsuntersuchung 
verwendet ahnlich einem Vorschlag von K. L. Stellmacher direkte Fourieransatze und bestimmt die 
Koeffizienten durch eine nachfolgende Stérungsentwicklung. Die Mittellagen Po der ausgewanderten 
Schwingungen und die Stabilitat dieser Schwingungen werden abhangig von der Intensitat (€ nach 
(24)) der Erschiitterung angegeben. Wahrend sich der Einflu8 hinreichend kleiner héherer Har- 
monischer in der Erschiitterungsfunktion als nicht wesentlich erwiesen hat, ist ein etwaiger Feder- 
aufzugsfehler von bestimmendem Einflu8 auf den stabilen Einsatzbereich des Gerates und den 
Verlauf der Auswanderung. Einige Versuchsreihen bestatigen diese wesentliche Abhangigkeit der 


Auswanderungserscheinung von den etwa vorhandenen Einbauungenauigkeiten der Feder auch 
experimentell. 
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Aus dem Institut fiir Mechanische Schwingungstechnik der Technischen Hochschule Karlsruhe 
(Eingegangen am 28. Marz 1957.) 
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Die Berechnung offener Rahmentragwerke nach dem Reduktionsverfahren 


Von S. Falk 


1. Einleitung. Die vorliegende Arbeit stellt die Fortsetzung und Verallgemeinerung meines 
ersten Aufsatzes! dar und behandelt offene, aus lauter geraden Stiicken zusammengesetzte Rahmen- 
tragwerke. Solche offenen Rahmen lassen sich durch einen einfachen Kunstgriff stets auf einen oder 
mehrere Durchlauftrager zuriickfiihren, und es zeigt sich, daB bei ebenen (raumlichen) Tragwerken 
héchstens drei (sechs) lineare Gleichungen aufzulésen sind, unabhaingig von der statischen oder 
geometrischen Unbestimmtheit und der sonstigen Beschaffenheit des Tragwerkes. An drei Bei- 


spielen wird das Verfahren vorgefiihrt. Formeln aus meiner ersten Arbeit, werden mit dem Buch- 
staben I zitiert. 


2. Der Begriff der Reduktion. Gegeben sei ein im Endpunkt E eingespannter starrer Trager 
nach Abb. 1. In jedem beliebigen Schnittpunkt S 1aBt sich die im Bereich A S angreifende gegebene 
Belastung reduzieren auf eine Kraftesumme St, und eine Momentensumme Siz. Beide Vektoren 
zerlegen wir in Tangentenrichtung (Langskraft und Drillmoment) und in zwei zueinander senk- 


M, 


ts os M7 
®; | g; 
~ ~ Ihe 
ie RM; Mu; 
A 
~ i; Smeal 
Abb. 1. Starrer durchlaufender Rahmen mit KraftgréBen. Abb. 2, Gerader Balken mit KraftgréBen. 


rechte Normalrichtungen (Querkrafte und Biegemomente). Diese sechs Komponenten werden vom 
Balkenteil A S auf den Teil S E tibertragen; wir nennen sie die KraftgréBen im Schnitt S. Die 
Reduktion der Kraftgréfen gelingt, weil die Vektoren 1, und Jt, im Anfangspunkt A gegeben 
sind; man braucht ja, in 4 beginnend, nur nach den Gleichgewichtsbedingungen 


®, = K_1+K 
Me Nt, + Mm Li, x @_, 


(1) 
fortzuschreiten (Abb. 2). 


Nun sei der Trager elastisch verformbar. Dann treten zu den sechs KraftgréBen noch sechs 
DeformationsgréBen hinzu, die diesen paarweise zugeordnet, ,,konjugiert“ sind, und zwar handelt 
es sich um die folgenden sechs Paare: 


Langsverschiebung ened) 
Riel Langskraft ; 
iE rn i Langsverdrehung (Drillung) oP, 
Langsmoment (Drillmoment) fe 


Querverschiebung (Durchbiegung) vv , 


ie ~ | Querkraft Q,> 9 
pc ichians Drehwinkel (Neigung) P > co) 
Biegemoment M, m 


Querverschiebung (Durchbiegung) w, 


: , | Querkraft Q.. 
2-Richtung Drehwinkel (Neigung) he 
Biegemoment Me 


Die Vorzeichen sind nach Abb. 3 gewahlt. Die KraftgréBen bezeichnen wir allgemein mit k,, die 
Deformationsgré$en mit d;. Wenn eine Unterscheidung nicht nottut, schreiben wir einfach g;. 


1 §. Falk, Ing.-Arch. 24 (1956) S. 216. 
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Den Begriff der Reduktion erweitern wir jetzt sinngemaB auf alle 12 GréBen (2): zu den sechs 
statischen Gleichungen (1) gesellen sich ja noch sechs geometrische, die mit jenen durch das Hooke- 
sche Gesetz verkniipft sind, und diese insgesamt zwélf Gleichungen geben eindeutig an, wie die 
zwolf GréBen g, von Punkt zu Punkt zu reduzieren sind. Allerdings — und das ist eine einschnei- 
dende Voraussetzung — mii®ten simtliche Kraft- und DeformationsgréBen im Anfangspunkt A 
gegeben sein, was in Wirklichkeit nie der Fall ist. Die Kernfrage des Reduktionsverfahrens lautet 
demnach: wie findet man die Kraft- und DeformationsgréBen im Anfangspunkt A eines Balkens ? 
Ihre Beantwortung stellen wir fiirs erste zuriick. 


3. Die grundlegenden Differentialgleichungen. Die GréBen (2) sind bekanntlich durch 12 lineare 
Differentialgleichungen erster Ordnung verkniipft, die sich zu den folgenden vier Differential- 
gleichungen zusammenfassen lassen: 


I 


—[E F(x) u(x)!" = p(x), 
[EJ ,(x) w(x)" = ¢.(x) » | 
—[G I(x) O(x)!" = tx) , | 

[E (x) v(x)" = qyla) » 

wobei Striche Ableitungen nach x bedeuten. 

Hier stehen auf den rechten Seiten die gegebenen Belastungsfunktionen p(x), q,(x), t(«) und 
q,(x). Wir wollen vorerst annehmen, y und z seien die ither den ganzen Integrationsbereich kon- 
stanten Haupttragheitsrichtungen, auch mégen 


die Steifigkeiten FF, EJ,, GIund EJ, konstant 


(3) 


YM, 


Abb. 3. Vorzeichen der Kraft- und DeformationsgréBen. Abb. 4, Gerader Balken mit Druck- und Drehfedern, 


sein (Voraussetzungen, von denen wir uns spater wieder befreien werden); dann namlich lassen 
sich die Lésungen von (3) unmittelbar hinschreiben. Als Integrationskonstanten wahlen wir die 
Kraft- und DeformationsgréBen an der Stelle x = 0 [vgl. (2) bis 1(9)]. 

Weiter sei im Punkt | des Balkens (x = 1) eine Anzahl von Druck- und Drehfedern angebracht 
(Abb. 4); die von diesen auf den Balken tibertragenen KraftgréBen sind nach dem Hookeschen 
Gesetz 

Ke 02. (tah, 2.21 Olas (4) 


wo die k; zu den d; konjugiert sind. Diese Federreaktionen miissen wir noch zu der im Bereich 
0 <x </ angreifenden Belastung hinzuschlagen und haben damit den fiir das Reduktionsver- 
fahren grundlegenden Zusammenhang zwischen den Kraft- und DeformationsgréRen. an zwei 
Schnitten 0 und 1 (genauer gesagt unmittelbar rechts davon) in der Entfernung / voneinander 


I [ p® dé = L(x), 
u(l) = u, = uy — EF Ly + u(l) , wich? 
A mit a Z (5) 
Liar ee hh Sean a ee de =u) 


te 
Hess 
— 


pl) = Wy 2 
M,(l) = M,, = Miyake) = Cue 
0.() = 0.1 = co + OD — 6, 0 (6) 
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mit eae) de == O42) [QAé) dé = M(x) 
0 
My) be @E 
| ey, dé = — (x), ae, ) dé = w(x) 
0 v oat 
f(&) dé = T(x) , 
WH) =0,=%—2, +H), ae | 
a mit eas (7) 
nid) = 1 = Ty) + Fi) —C,9, -— oe dé = §(x) | 
‘oa 
I? B & 
ol) =, =% +ly— 5 EJs M,9 + CEJ Qyo + v(l) , } 
l i? 2s 
y)) =y= 5 A ey a a My + CEJ. Oyo + pl) . 
M,(!) = M,, mae Mi LQyo sa M,(1) Cy VY ? 
Q,(2) == Qo a Ovo are Tie = Cc, Vy c (8) 
mit f q() dé = Q(x). J O,(é) dé = — M(x), 
QO 0 
Mey cs * le Ps 
- da aa) | one) 


0 
4. Die Leitmatrizen. Nun schreiben wir die Gleichungen (5) bis (8) noch einmal in Matrizen 
hin und bekommen, wenn wir, um allgemeiner zu bleiben, die Zeiger 0 und | durch i — 1 und i 
ersetzen, ferner aus formalen Griinden die triviale Gleichung 1 =1 zu jedem der vier Gleichungs- 
systeme hinzufiigen: 
te RR a 


L “Ud 


9, = Lyi 0a» 
: (9 
Ie Vita ) 


2 oJ 
bc = eidia- 


Die vier sogenannten ,,Leitmatrizen“ mit ihren zugehérigen Vektoren sind dabei der Reihe nach 


Dee ze TAREE i tee 0 u 
Y= 0 | L | eae Ue [Eee *€00) 
0 0 pe ae y 
haw shel IDOE Sede TGS pyr 40: ann Ot 310 © 
0 Ame REVUES MS POINT gt sepa + Oto ia 9 
xv, =| 0 0) il l M, Im alanis See ee (id) 
0 0 “0 ] Q. _—— Cw 0 ke 
0 0 0 0 ett cael) ale | 
1 sah 0 ) oe 5) 
Sp : OS ale ee eye (12) 
0 rath ee “i 
i l —/ ae Sd Be 3/6 E A vo sit 0 0 \ ie \ 
iy ge ap 2/2 E ae 0 0 y 
aes | 0 0 I 5 M.1 0 —C, pen NE (13) 
Ouaped 0 1 0, | —% 0 Q, 
lo 0 0 0 Peet 0 
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Die von der gegebenen Belastung herriihrenden, mit dem Zeichen ~ versehenen ,,LastgréBen“ 
sind fiir die praktisch wichtigsten Falle in den Tabellen 1 und I (1) zusammengestellt. 

Aus den Matrizen (10) bis (13) lassen sich nun leicht Leitmatrizen héherer Ordnung zusammen- 
stellen, zum Beispiel fiir den Sonderfall der ebenen Statik (Biegen in der (z, x)-Ebene mit Langs- 
dehnung in x-Richtung) 


¢ 0 0 0 Oe Veal ae 0 0 0 u 
gees bea Ss) PES, BES, 0 Ww 0 0 0 w 
0 0 L =e eg. 20 p 0 0 @ 20 Y 
0 0 0 1 l 0 M,; 0° 0 = Cree | Mele 
0 0 0 0 1 0 0. 0 C0 Q. 
Ou s0 0 0 0 1 iD eu 0 L 
0 0 0 0 0 0 1 0 0 0 1 


oder ganz allgemein: sind 9 Paare an der Deformation des Balkens beteiligt, so kann man diese 
gemeinsam mit der Ziffer 1 zu einem Vektor zusammenfassen und aus den Matrizen (10) bis (13) 
die zugehérige Leitmatrix mit 2 9 + 1 Zeilen und 3 9 + 1 Spalten ohne Schwierigkeiten zusammen- 
bauen. Alle diese Leitmatrizen sind offensichtlich von der folgenden typischen Form: 


i Ie \ 1, | 
5 oe . 5) D . 
| 0) a) e oe 0 a 
(Qu = | - --- == ‘ = (15) 
~~ k | ? 5 t> 
oe oh) en f ky 
| 7 hy 
eA Ome O a Oy malas O 0) 1 1 


oder noch einmal ausfiihrlich in Untermatrizen der Ordnung 0 geschrieben (9 ist die Nullmatrix): 
Gide te te. =(b—)),;, | 
Gt.+C (—d),=¢—-f,. J 


Jede dieser Untermatrizen enthalt eine gewisse Gruppe von Kenngréfen, und zwar gilt folgende 
Zuordnung : 


(16) 


(Ss > Geometrie , SM —> Kraftegeometrie , | a7) 


vy —> Steifigkeit des Balkens , & — Federung von auBen . J 


In St finden wir 0 der Gleichungen (1) wieder; %} verschwindet, wenn die Steifigkeiten unendlich 
groB werden (starrer Balken), ©, wenn die auBeren Federn fehlen. Den punktiert (gestrichelt) 
eingerahmten Teil von <' wollen wir den Kern (den Hauptteil) der Leitmatrix nennen. Die letzte 
Zeile hei®t (aus spater verstandlichen Griinden) die Kontrollzeile, die die gegebenen LastgréBen 
enthaltende Spalte die Lastspalte. Die Determinanten der Kerne und Hauptteile aller Leitmatrizen 
haben offenbar den Wert 1. Die Reziproken der Kerne entstehen, wie man leicht zeigt, wenn man x 
durch — x baw. | durch —1/ ersetzt. 

Die Matrizen (11) und (13) lassen sich iibrigens noch um die Ordnung zwei erniedrigen, wenn 
man irgendeine GréBe g; mit Hilfe ihrer Konjugierten, die dann in der Lastspalte erscheint, eli- 
miniert. So entsteht zum Beispiel, wenn man in (11) die Querkrafte entfernt, die folgende Leit- 
matrix, die im wesentlichen den gesamten Formelapparat der Traversenmethode von Stewart! 
enthalt : 


—2 Y2E J, a + 3 (w;, — w;4)/l : 0 v2 
Qh ees : ~ wee 
<i | 6D fl 2). \M,-.6 FS) (w= en ne eee (i 
0 0 1 , 0 1 
mit 
G=9+30/l, M,=M,—6EJ, ar. (19) 


1 R.W. Stewart und A. Kleinlogel, Die Traversenmethode, Berlin 1952. 
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Solche, an der Deformation des Balkens zwar beteiligten, aber aus dem Kern der Leitmatrix ent- 
fernten Paare nennen wir ,,unterdriickt“. Falls von Interesse, Jassen sich die unterdriickten GréBen 
nachtraglich leicht wiedergewinnen; zum Beispiel findet man aus (11) fiir die Querkrafte Q, , und 
Q;, wenn man der Ubersichtlichkeit halber die Indizes z und y bei Q und M fortlaBt So 


Q; Sligo Q; Q; ea ; (M; M, M;.4) ? (20) 


ferner fiir die Durchbiegung w als Funktion von x im ersten Feld 


a i x ei: x ESS oe | 
(x)= x (1 ; (: tS Jor . EJ My | Fe) == i w(l) + wy + (w, — wo) 2 (21) 


und entsprechend in allen anderen Feldern. 


SchlieBlich fiigen wir noch eine Bemerkung zu der Matrizenmultiplikation 
= Ut (22) 


ein. Die Leitmatrix %; hat 30 + 1 Spalten, aber der Vektor tj nur 2@ + 1 Komponenten; 
man faBt daher die ersten g@ Komponenten von 7; als die letzten 9 des Vektors y,_, auf und alte 
pliziert sodann wie gewéhnlich. Kine solche Matrizenmultiplikation nennen wir ,,verschrankt“. 
Zwar kénnte man die Gleichung (22) auch auflésen in die beiden Gleichungen . 


t= O;t-4> t=U,if (23) 


mit dem quadratischen Hauptteil ; und der ebenfalls quadratischen ,, Ubergangsmatrix“ 11, [I (13) 
in aad Erweiterung]. was aber nur unnétige Schreib- und Rechenarbeit verursachen 
wiirde. 


Tabelle 1. Lastgréfen fiir Langsdehnung, Drillen und Biegen mit unterdriicktem Q, w 
| 


p p 7 t 
® Son alee ® eee eae) 
a hh a ae 
u(1) ie Polk | 2p ht) BF nO) | =TaGl | =th/261 
cy a P ph 7K) an 7 | th 


®@ eee Pee _ cio 


PEG Bi a = 


te —Pa a\ Mb 36 

bm ees ed = 3/6 
El) = ( r| pa ( 31] qh8/24 BJ 
= a® b? 2 
M(l) Pa(t— fr] m(1—3 7) qh2/4 


5. Der durchlaufende Trager. Der geradlinige Abschnitt zwischen zwei benachbarten Feld- 
grenzen heiBt das Feld der Nummer i. Ein Punkt F ist Feldgrenze, wenn fiir ihn mindestens eine 
der folgenden Eigenschaften zutrifft: 

1) In F tiitt eine Federreaktion nach (4) hinzu. (24) 

2) Fiir z der Kraft- oder DeformationsgréBen sind in F feste Werte vorgeschrieben (2. B: 
die Werte Null), oder allgemeiner: zwischen den Kraft- und Deformationsgréen im Punkte 
sollen z lineare Bedingungen erfillt sein. 


3) F ist ein Verzweigungspunkt. 
4) F ist ein Knickpunkt. (27) 


Demnach stellt also auch jedes freie oder irgendwie gegen die feste Umgebung gestiitzte Ende 
eines Tragers eine Feldgrenze dar. Mitunter kann es zweckmafig sein, auch andere Punkte zu 
Feldgrenzen zu erklaren, etwa Punkte, in denen die Steifigkeitsfunktion unstetig ist (abgesetzte 


Wellen und dergleichen). 


9) 
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Fiir einen durchlaufenden Trager mit n Feldern (Abb. 5) gilt nun nach (9), (15) usw. der Reihe 


nach 


t= it; te= hee ot = SG (28) 
oder eingesetzt 
tate a, Caer Seer wT Xo “ (29) 
ee 
in = Ba Eo (30) 


Sind nun an der Deformation Paare beteiligt, so miissen in den Endpunkten A und FE auch genau 


je @ Randbedingungen der Form : 
Bo) 0 Lee Oy (31) 


Sun = O(“ =1,2,.--@) be 


(bzw. Linearkombinationen dieser Gleichungen) befriedigt werden, und zwar muB in A wie in E 
aus jedem Paar gerade eine GréBe g verschwinden. Die in A verbleibenden g wollen wir die ,,Frei- 
gréBen“ (A), Ay,...A,) = a und die in E verschwindenden die ,,NullgréBen“ (B,, B,...B,) =b 
nennen. Aus (30) und (32) entsteht dann das lineare inhomogene Gleichungssystem 


Wack =p 0s (33) 
Damit haben wir die Kernfrage des Reduktionsverfahrens (Ziff. 2) fiir den Durchlauftrager — und 
wie wir sehen werden, im Prinzip bereits fiir jedes beliebige Tragwerk — beantwortet: Von den 


2 0 Kraft- und Deformationsgréfen im Anfangspunkt A eines durchlaufenden Tragers sind g gleich 
Null; die itibrigen 9 gewinnt man aus dem linearen Gleichungssystem (33). 


Abb. 5. Durchlaufender Rahmen, schematisch. Abb. 6. Knickpunkt eines Rahmens mit begleitendem Dreibein. 


Den praktischen Hergang der Rechnung zeigt das Schema (34): Die mit den Freigré8en A; 
und der Ziffer 1 multiplizierten 9 + 1 Spalten des Anfangsvektors x, werden unabhangig von- 
einander an sdmtlichen Leitmatrizen vorbeimultipliziert bis unten rechts neben Ue die o Null- 
gréBen B; erscheinen. Dann wird des Gleichungssystem (33) der Ordnung g9 gelést und mit dem 
nun bekannten Anfangsvektor y) die Rechnung wiederholt, wobei die Kraft- und Deformations- 
gréBen unmittelbar rechts von jeder Feldgrenze zahlenmaBig erscheinen. Nach (5) bis (8) sind 
dann auch Kraft- und DeformationsgréBen als Funktionen der Balkenkoordinate x bekannt, 
womit die gestellte Aufgabe erschépfend gelést ist. 


A, fe peered lai 
1 0 0 0 
0 1 0 0 
Yo Ste RRP UR RIE Metis His oc Be’ cc 
0 0 as aeeey 1 (34) 


t 
(2) (a 


Sind alle Feldgrenzen vom Typ (24), so sind die Elemente samtlicher Leitmatrizen von vornherein 
gegeben und die Rechnung bereitet nach (28) und (34) keine Schwierigkeiten. Zwischenbedingungen 
nach (25) kann man auf drei verschiedene Arten erfiillen: durch Ablésen der FreigréBen, das ab- 
gekiirzte und das erweiterte Verfahren [I (8)]. Bei abgewinkelten Tragern schlieBlich (Abb. 6) 
transformieren sich die GréBen (u,v, w) und ebenso (L, Q,>9,). (9, yp) (7; M,, M,) mit Hilfe 


1 


ie 
y 
t 
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der Matrix 


Cp0y Cy 0g ~ 6s 057 (35) 
Cf Ga. Kes iy ey Gy) 

Je nach Art und Zahl der beteiligten Paare werden zweimal @ Elemente dieser Matrix zu einer 

Gesamttransformationsmatrix S zusam- 


mengefaft; es gilt dann beim Ubergang 
an einer Knickstelle i NV 


G—=Uto> (36) = : as 
Sy Cr alis TTT: 


to = Gn tin = 
und zwar empfiehlt es sich, von vorn- 
herein die Kerne der Leitmatrizen von 
rechts mit I zu multiplizieren. 


Abb. 7. Ersatz eines durchlaufenden Rahmens durch eine Koppelfeder. 


6. Federmatrizen. Wir fragen jetzt, wie die 9 KraftgréBen mit den 0 DeformationsgréBen in 
einem beliebigen Schnitt N eines Durchlauftragers zusammenhingen (Abb. 7). Rechnen wir den 
Trager bis zum Punkt N durch, so erscheint unten rechts in (34) folgendes Zahlenschema: 


Tis eds Fen. Hoo 1 

* * * \ d, ( d, \ 

* * * d, d, 

bee eee kh ; 

* * * d, d, 

* * * i a k =Utyn- (37) 
ers 3 

sil bese ieteneh td S + 

* * jaa * iy hy 

x * th * ky ky 

0 0 Sis 0 1 il N 


Mit den o-reihigen Matrizen Si und © und den Vektoren 
= (4i435). 2 :'A,)5 
Dred, 0s wi: 0), (38) 
1 (Kip lep eae Fi) 


kénnen wir dafiir kiirzer schreiben 


Ra=d—dD,. Ga=f—t (39) 

und weiter, da $i stets regular ist, ~ 
a= (b—d), (40) 
f—f=GR1(d—d) =C(0—?). (41) 


Die iiberstrichenen GréBen stammen aus der gegebenen Belastung; besteht der Trager aus einem 
einzigen Feld, so stimmen sie natiirlich mit den LastgréBen der Tabellen 1 und I (1) tiberein. Die 


Matrix 


C15 cA sat 
A C2 cae C22 Cio 9 
teas eat Pant) ond sigue ©) ae) 
Coo 20 Go | 


heiBt die Federmatrix des Rahmenschnittes N, ihre Reziproke C+, falls vorhanden, die EinfluB- 
matrix. Nach dem Satz von Maxwell- Betti ist die Federmatrix, wie in (42) bereits ausgedriickt, 
symmetrisch zur Nebendiagonale; die Elemente c;; sind negativ oder gleich Null. Diese Feder- 
matrix stellt die natiirliche Verallgemeinerung der Matrix © in (15) dar. Dort waren noch alle 
o Gleichungen nach (4) entkoppelt. Wollten wir hier aber die elastische Wirkung des abgeschnit- 
tenen Rahmenteiles im Punkte N durch eine einzige Feder ersetzen, so miiBte diese die Matrix 
(42) haben, lieBe sich also nicht einfach aus Druck- und Drehfedern wie in Abb. 4 zusammensetzen, 
5%* 
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sondern erst verwirklichen durch ein komplizierteres (gedachtes) Gebilde, das wir als Koppelfeder 
bezeichnen und nach Abb. 7 symbolisieren wollen. Die auf den Punkt N von der Feder tibertra- 
genen Anteile der gegebenen Belastung deutet ein Pfeil an; wir lassen ihn fort, wenn samtliche 
iiberstrichenen GréBen in (37) verschwinden. 

Da nun gewisse Bauelemente der Praxis immer wieder vorkommen, lohnt es sich, Kataloge 
von Koppelfedern zusammenzustellen. So zeigt die Tabelle 2 Federmatrizen nebst Kehrmatrizrn 
jt am einfeldrigen Balken mit konstanter Biegesteifigkeit FE J, fiir die Paare (M,, p) und (Q,, w). 
Hier kann man die LastgréBen direkt aus den Tabellen 1 bzw. I.1 ubernehmen; in allen anderen 
Fallen dagegen mu man sie nach (37) berechnen, da das Katalogisieren von LastgroBen kaum 
praktischen Wert besitzt. 


Abb. 8. Zuriickfihrung eines mehrfach verzweigten Rahmens auf einen Durchlauftrager. 


7. Das Anfedern von Rahmenteilen. Jetzt sind wir in der Lage, einen beliebigen offenen Rahmen 
auf eine Anzahl von Durchlauftragern zuriickzufiihren: wir erklaren irgendeinen durchlaufenden 
Linienzug des Tragwerkes zum Haupttrager, schneiden simtliche Seitenzweige ab, ersetzen sie in 
~ den Schnittpunkten durch ihre Koppelfedern, deren Matrizen in (15) unten rechts einzutragen 
sind, und kénnen nun den Haupttrager wie iiblich berechnen. Dann aber liegen nach (40) auch die 
FreigréBen der Seitenzweige fest, womit die Aufgabe gelést ist. Auch bei mehrfachen Verzweigun- 
gen nach Abb. 8 kommt man stets zum Ziel. Bei schiefwinkligen Anschliissen sind Kraft- und 
DeformationsgréBen in die Richtungen des Haupttragers zu zerlegen. 


8. Anfedern mit Nebenbedingungen. Bei biegesteifen Anschliissen sind die DeformationsgréBen 
des abgeschnittenen Zweiges gleich denen des Haupttragers, so daB die in diesen eingeleiteten 
KraftgréBen als linke Seiten von (41) unmittelbar gegeben sind. Sollen aber z (1 <z <Q) der 
KraftgréBen fest vorgeschriebene Werte annehmen (etwa die Werte Null wie in den Anschliissen 
der Abb. 9), so werden fiir den Haupttrager auch gerade die z konjugierten DeformationsgréBen 
des Seitenzweiges bedeutungslos; diese eliminieren wir daher in den Gleichungen (41) mit Hilfe 
geeigneter Zeilenkombinationen, wodurch eine um die Ordnung z erniedrigte, ebenfalls neben- 
symmetrische ,,reduzierte’ Federmatrix entsteht. Ahnlich geht man vor, wenn 2 lineare Bedin- 


Tabelle 2, Federmatrizen fiir einfeldrige Balken mit konstanter Biegesteifigkeit EJ, Lange | 


Art der Stiitzung | | | 


bees G2) CO b, Ca = K-1 (db —d), Qa 
Ge, Wig iva, Oe | Federmatrix © | Matrix {1 | Freie Konstanten 
ban (iia, + OF ail +6/2 + 2/1 | M 
=: £ aa! ai E 
ag ane & ie 12/P = 6/? (= 12/P = ds | (one 
2a 2b Ba(* 3/P —3/l Way eg Be 50 iti. 
aa e& a& 3/8 = 3/P — 3 EV/P = SE IIe) | Q/a 
nie aaa el yee tere eee a iy, 
1 3a 3b m EI / = [}) | w 
a G @ a 0 0 0 —E Jl | M/. 
ay eS esa sa | cals? Gaines Hea 
a a Gace Woes mace | 


Die unteren Vorzeichen in & und S{~! gelten fiir die Falle b) 


a 
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ungsgleichungen zwischen . 6 : : : - 
sunese g den KraftgréBen oder noch allgemeiner z Gleichungen zwischen samt- 


lichen 20 Kraft- und DeformationsgréBen vorgeschrieben sind; doch hat dieser Fall wohl nur 
theoretisches Interesse. 


; 9. Die Addition (Parallelschaltung) von Koppelfedern. Die freien Enden yon p Durchlauf- 
tragern mégen im Punkte H biegesteif miteinander verbunden sein (Abb. 10). Dann gilt fiir 


diesen Punkt nach (41), wenn wir uns saémtliche Kraft- und DeformationsgréBen nach gemein- 
samen Richtungen zerlegt denken, 


(43) 
oder 
4A. 
mit ) 
ie (45) 
=0 H 
Sona aan es 
CALLS, 


Abb. 9. Rahmenanschliisse mit Nebenbedingungen. Abb. 10. Ersatz eines offenen Rahmens durch eine Koppelfeder. 


Man kann daher den ganzen Rahmen der Abb. 10 durch eine einzige (belastete) Koppelfeder im 
Punkte H ersetzen, und zwar auch dann, wenn infolge Nebenbedingungen einige oder alle Feder- 
matrizen reduziert sein sollten. Sind insbesondere alle p Zweige unbelastet, so wird einfach 


fy = (2'G;) Dy = ©, dy (46) 


Das heiBt: die von der gewéhnlichen Druck- und Drehfeder her bekannte Tatsache, daB sich bei 
Parallelschaltungen die Federkonstanten addieren, iibertragt sich auch auf die Federmatrizen. 


LaBt man jetzt in H gegebene KraftgréBen f angreifen, so folgt aus der Gleichgewichtsbedingung 


f +- &, On —() (47) 
sofort 


dy =—Cpf, (48) 


und damit ist nach (41) auch schon der ganze Rahmen berechnet: ein Vorgehen, das in diesem 
einfachen Falle stark an die Deformationsmethode erinnert. Auch dort werden ja die 9 aus dem 
Zusammenhang eliminierten DeformationsgréBen d;,, aus g linearen inhomogenen Gleichungen 
ermittelt. 


10. Rahmen mit unverschieblichen Knoten. Sind die Knoten (Feldgrenzen) eines Tragwerkes 
unverschieblich, so kann man alle Kraft-Verschicbungspaare unterdriicken und allein mit den 
Momenten-Drehwinkel-Paaren rechnen. Bei ebenen Tragern hat man dann lediglich das Paar 
(M,,y) mit der Leitmatrix (18), in der die w; verschwinden. Aus der Momentenlinie folgt der 
Querkraftverlauf durch Ableiten (20), die Biegelinie aus (21). 

11. Praktische Durchfiihrung und Beispiele. Das Einflechten der Federkonstanten nach I. (22) 
ist nun durch die verschrankte Matrizenmultiplikation iiberfliissig geworden. Um auch die Rechen- 
proben mit in das Schema (34) einzubeziehen, addieren wir zur letzten Zeile der Leitmatix die 
negativen Summen der ersten 2 Zeilen, so daB die Zeilensummen dieser neuen Matrix gleich 
Null bzw. Eins sind. Addiert man nun zur letzten Komponente des Vektors y; die Summe seiner 
ersten 2 9-Komponenten, so kommt bei richtiger Rechnung unten in jeder der homogenen, mit 
einer Freigré8e multiplizierten Spalte eine Null, in der inhogenen (Last-) Spalte dagegen eine Eins 
heraus, gegebenenfalls bis auf vertretbare Abrundungsfehler, die zu streichen sind, denn weiter- 
rechnen darf man natiirlich nur mit den Zablen 0 oder 1. Weitere willkommene Kontrollen stellen 
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die Bedingungen c,, <0 und die Nebensymmetrie der Federmatrizen dar (42). Auch sollte man 
immer das Gleichgewicht am Tragwerk iiberprifen. , 

Bei ebenen Tragwerken zeichnen wir positive Biegemomente nach oben auf (negative w, Q- 
Richtung); sie erscheinen dann auf der Zugseite des Balkens. Positive Langskrafte bedeuten 
Druck-, negative Zugkrafte. 

a) Erstes Beispiel. Fiir ein offenes Tragwerk mit unverschieblichen Knoten (Abb. 11) 
ist die EinfluBlinie fiir das Biegemoment M,, unmittelbar rechts vom Pfosten 2 infolge einer wan- 
dernden Einzellast P zu berechnen. 


+ |\-M 
y D. (as 
_ C, C3 
> : 
ZL, DX; L: 
>’ a (a 


Abb, 11. Rahmen mit wandernder Einzellast, dazu Ersatztrager. 


Wir bringen rechts vom Pfosten 2 ein Gelenk an und lassen an ihm ein Momentenpaar +M, 
—M angreifen; aus der dadurch hervorgerufenen Biegelinie w(x) folgt dann bekanntlich 


Me Oe (1.01) 
Rs a 
wo @ die Winkeldifferenz (SprunggréBe, Tabelle I,3) im Gelenk und w(x) die Durchbiegung an der 
Laststelle ist. Da die Knoten unverschieblich sind, rechnen wir mit der Leitmatrix (18), in der alle 
w; = 0 zusetzen sind. Mit 
ae 


werden die Drehfederkonstanten der Pfosten nach Tab. 2 (Elemente c,, rechts oben) 


dh. I=1, EJ => (1.02) 


Pfosten 1 (Fall 2a): Cc, = 3 —— a 


(1.03 
Pfosten 2 und 3 (Fall la): C, = C, =4 zed 


Als auBere Belastung haben wir im Feld 2 das Moment +M = P* 1* = 1 mit b = 0 und im Feld 3 
das Moment —M = —1 mit 6 =/=1 (Tabelle 1,3); in die Lastspalte des dritten Feldes ist 
auBerdem die SprunggréBe ~ aufzunehmen. Freigrée ist hier das Einspannmoment M). Damit 
laBt sich das Zahlenschema | hinschreiben; die Zeilensummen in den Leitmatrizen sind der Uber- 
sichtlichkeit wegen fortgelassen worden. Die Elemente rechts von {,, zum Beispiel finden wir so 


(—2)- (1) + (2): (—5) +040: (*) =— 12, 
(1,5)- (1) + (2) (5) 41-04 (—4) 2) 59,8. 
Die Zwischenbedingung M, = 0 (Gerbergelenk!) tibergehen wir zunachst, rechnen also nach dem 


erweiterten Verfahren I (8c). Aus den beiden eingerahmten Gleichungen ergibt sich dann zum 
SchluB 


1 


M,= 595M,+ 1 =0>M,=—~,=—0,0168, 
i 1.04 
My 450M = 0 eee Oy 


Damit lat sich nun die rechte (End-) Spalte des Schemas hinschreiben, womit Drehwinkel und 
Biegemomente fiir alle Feldgrenzen gegeben sind (Zahlenschema 1). 
Nach (21) sind jetzt die Biegelinien in allen drei Feldern wegen w, = 0, 1=1, 


Por = Yu + P = 0,2017 — 0,7472 = — 0,5455, EJ = > 
im Feld 1: w(x) = — x (1 — x) [(1 + x) - 0 — x(—0,0168)] +0, 
Feld 2: w(x) = x(1 Z| (! L 2) ( 0,0168) — x(0,0840)| ois (1.05) 


Feld 3: w(x) = — x (1 — x) [(1 +x) (—0,5455) — x - 0] + [(—1) 2? — «3(— 1) L]. 
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M, 1 


Zahlenschema 1 0 0 0 Po 
( tl 0 [—noss ) = (ri 
0 1 il 1 
—2 1 Orn lis 0 cata: 0 — 0,0168 V1 
oo —2 0 | —3 —5 0 0,0840 } = (i 
0 0 Ly 0 0 ae, 1 1 
—2 2 Gg). 0 ete 0 0,2017 Pa 
Fe | Toe ae 1-1 —s| ( 505 — I ) 9,000 = (¥. 
0 0 oe 0, at nalts 1 1 
—2 1—2q—1 | 0 83,5 29 0,0910 Ps 
o( 3 —2 3942 | — ‘ ( 489 19) $000 = (s.) 
| 0 0 | 0 0) at 1 1 
Das nach (1.01) durch —@ = 0,7472 dividiert und geordnet gibt 
Feldjl: M,, = (0,0225 & — 0,0225 &) Pl, 
Feld 2; M,, = (—0,0618 & + 0,1124 £ + 0,0225 £) Pl, (é = 7 (1.06) 
Feld 3: M,, = (0,6085 & — 1,3383 & + 0,7301 é) Pl. 
Den SOS a zeigt Abb, 12. zy Ly Dy 
b) Zweites Beispiel. Das offene 


verzweigte Tragwerk von Abb. 13 ist 
zu berechnen. Wie ist die Rechnung 
abzuadndern, wenn die Anschliisse in B 
und C nicht biegesteif, sondern gelenkig 
sein sollen? Das Tragwerk sei in Langs- 
richtung starr (1/E F= 0), die Biege- 
steifigkeit des Mittelpfostens ABCD sei 
2 E J, die aller iibrigen Felder E J. Ge- 
gebene Werte sind 
t= 2.000... BE J.— 80 tm? , 

q= 11 t/m ean 
mor ait PJ P© 2K Jil, “and 
wenn wir im Hinblick auf die Nenner 
in Tabelle I,1 vorerst 


Mpp [PL] 


Ce 


Po # PK 
dF 24 P |! = 480 t/m Abb. 12. Momenten-EinfluBlinie fiir den Rahmen der Abb, 11. 
wahlen, 


Pleas S125 g=24, dP =gip4, Pt =qhi24. (2.02) 
Das Schema der Anfederung zeigt die Abb. 13 rechts. 
1. Berechnung des Zweiges I nach Abb. 14. Die AnschluBkrafte und -momente der Pfosten 


12 und 13 finden wir in Tabelle 2, Fall 1a vor; da beide Pfosten unbelastet sind, verschwinden 
die ~-Glieder, und es verbleibt 


(o)= i Be aD a Aa (")= ee m) ( \= (x Sse Lh (i = 12,13). (2.03) 


Der Pfosten 11 ist gelenkig angeschlossen (Nebenbedingung!) wir eliminieren daher den Winkel Pu: 
indem wir die mit —1,5 multiplizierte erste Zeile der Federmatrix (2.03) zur zweiten addieren 


und M,, = 0 setzen. 


pa ae 7a ala (2.04) 
0-0), \-6 —3/ \prohn ~ 


100); ‘1,5(0 — My) oe ae. 1,5(w zis 7) os ee @ —15 M +1,5 w), — 1,5 wy 
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oder mit den Werten aus Tabelle I,1, Spalte 3 und aus Formel (2.02) 
Q,, = (24 — 1,5 - 24/2 + 1,5 - 24/12) —15 wy =9—1,5 uy. (2.05) 
Da die w,Q-Richtungen der Pfosten mit der (L, u)-Richtung des waagerechten Tragers I zusammen- 


fallen, setzen wir in (2.03) und (2.05) wy = Wy, = Wy3 = Uz, ferner Py = Py Pis = P2 und bekom- 
men als eingeleitete KraftgréBen (die eingeklammerten Werte iibersehen wir vorerst) 


VE ea ee) » Oy =9— 1 a; (= — 8,467) , 
My = — 3 us — 2 gy, (= — 0,721), Qe = — 6u,;— 3, (= — 1,614), (2.06) 
Mig = —3 vy — 24 (= — 0.882) Qu — 6 5 — 8 go(= — L781) 


Ly = 50; =9 — 13,5 us —3 (g, + 92) (= 5,072). (2.07) 


Abb. 13. Mehrfach verzweigter offener Rahmen mit Haupttrager. 


Somit hat der Zweig I das Aussehen der Abb. 14, rechts. Wir unterdriicken die Paare (L,u) und 
(Q, w) und rechnen mit den Leitmatrizen (18), wobei die w; = 0 zu setzen sind. Unten rechts im 
Zahlenschema 2.1 finden wir neben %;, die beiden Gleichungen 


— 649, + 15u,=93; 117 g)—27 us = M; 


vor und berechnen daraus 


My == 1,926 0, 0422 ae (2.08) 
Go = — 0,0156 ws + 0,2344 us . (2.09) 


My=-3u;  My=-3u, 
—— 


: : Us, L3 
y 
om 7 


Abb. 14, Teil I des Rahmens von Abb. 13. 


Mit den Werten py, = — 2 Gp, 2 = 11 Y) — 3 uz aus dem Schema 2.I wird nun nach (2.07) und 
(2.09) 
L, — 9 = — 13,5 ug — 3 (9 gy — 3 us) = — 4,5 ug — 27 (—0,0156 g, + 0,2344 ug) 

= — 10,828 u; + 0,422 », , (2.10) 


und das gibt zusammen mit (2.08) 
M, P ( 0,422 —1,828 us) 
Deo ee 0.698 0.422 a (2A) 
Wie nach (42) zu erwarten, ist diese Federmatrix nebensymmetrisch; die Elemente c,, sind negativ. 
2. Berechnung des Zweiges II. Wir stellen den abgeschnittenen Zweig auf den Kopf, damit 


das freie Ende rechts liegt (Abb. 15), was zwar nicht unbedingt notwendig, aber fiir die Anschau- 
ung bequem ist. Fiir die Pfosten 14 und 16 wird nach Tabelle 2, 2a 


+ _ 1 /—3/0,52 —3/0,5 SSG eer 
G,=G, = 2 iz 3/0,53 — ve) = a 12 cal A (2.12) 
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und fiir die Pfosten 15 und 17 nach Tabelle 2, 1b): 


l 6 aA eae 
aes Z. 2 
Gs=7(_y 4 Sie eae 


3 alt 6/1,52 — 4/1,5 io 4/3 — 4/3 
EEN P7947 58 ets) = (16) 4/3) 


Da Belastungsglieder nicht vorhanden sind und wegen Wi4 = Wiz = Wig = Wy, = Ug, ferner Y, 


= Pis = Pr Pis = Viz = Ys sind die in den Punkten 4 und 5 in den Trager II eingeleiteten Kraft- 
gréBen 


My, =—6u,—3q,(= 6,954); Qi, =— 12 u, — 6g, (= 13,908) , 
M,,= 3ug—2q,(=—5,499); Q,, =—6 ug +3, (= 10,420), 
M,, =—6ug—395(= 7,787); Qjg = — 12 ug — 6g, (= 15,574), (2.14) 
4 4 4 19 4 
M,, = = Ver aes (= — 2,330); Q,, =— 6 U6 3 V5 (= 2,974) . | 
L,= >'0,= 31, Ug—3 M4 45 s (= 42,876) . (2.15) 


Die Drehfederkonstanten sind demnach C,=3+4+2=5 und C, = 3 + 4/3 = 13/3 (Abb. 15 
rechts). Das Zahlenschema 2.II liefert die beiden Gleichungen 


1 ind 
57 5 Pa + 205 ug Pes 


2 ne 
L075 Pa 375 ug = M,. 


das heiBt 


(LN GaN Bele 
M, =— 1,878 9, + 0,517 u,, (2.16) —> Fi ales 
@, = 90,0174 o,— 0,355 ug. (2.17) 77 Zz 
Nun setzen wir in (2.15) ein und er- 
halten Abb. 15. Teil If des Rahmens von Abb. 13. 
L, = — 31,778 u, — 3 gy, — 4,667 (— 7 go, — 3 ug) 
= — 17,778 ug + 29,667 (0,0174 ~, — 0,355 ug) = — 28,299 ug + 0,917 Gg » (2.18) 
also zusammen mit (2.16) 
M, Be 0,517 ##—1,878 (Fr (2.19) 
L, — 28,299 0,517) \@5 
/ Po 1 
i 0) \ 0,0863 
Zahlenschema 2.1 () 0) 0) Site 
0 ey I 
ey) i () 0 f —2 Aes — 0,173 
2, 3 eee -2) ( Tle 3 ¥ eR Wok Waa te 
—1 lL Laps iy 4 2 0) 1 I 
sae) 1 0 | ) 1] —3u, — 0,117 
or ee PLS 6 Ging ) (- 42 Gael (2504965 |. 24, 
—] Re ae a 2 0) 1 1 
—2 if 0) /— 64 1S ats — 0,193 ; 
Si Elem eaee | | ivy — 27 uz 0,503 | = 3, 
: —] 1 1 0) haute 1 


3. Berechnung des Haupttragers ABCD. Nach Abb. 16 ist in den AnschluBpunkten B und C 
zu setzen . 
Ug = Wg, We=0, Pe=—Psi Ug =— My, We =—9, Ye =Po- (2.20) 
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Ya 1 
1 0 0,578 
Zahlenschema 2.11 —5 —3 U, LAS ole 
0) 1 1 : 
—2 ] 0) Le) —7 —3 Us, 0,300 
eal Age et 4,279 
Qs HO 3 Ug | — 13/3 33 3 Ug B ¥ i. : 
14 
—l1 {i 1+ 3 Ug | 13/8 0 i i 
1 it 
—2 1 0 573 =. 203 ug 3,681 
g 2 2 = ‘ 
ae 3 —2 0 — 107, — 373 us — 7,660 Is 
—1 1 1 0 1 1 


Somit wird aus (2.11) und (2.19), wenn wir beachten, da®B — L, (nicht + L,) mit der positiven 
Q, w-Richtung des Haupttragers zusammenfallt, 


MeN 0,422" Sd, 828\ 0. ae = 0,517, =A BIB es 2.21) 
D0) 10928 0,422 fs > Wah) een 300" Fe oreiT are 


Diese beiden Federmatrizen sind daher in 
LY, und &, (Schema 2.1II) einzusetzen. Aus 
den eingerahmten Bedingungsgleichungen 
Yi = 90 und M, = 0 berechnen wir die 
beiden FreigréBen zu 


M, =1,7452, Q,=—3,1044 (2.22) 


Ma e Me 
ea © 
ot 
cee 
a Qe 


Abb. 16. Mittelpfosten ABCD des Rahmens von Abb. 13. Abb. 17. Gerader Balken mit KraftgréBen,. 


und wiederholen mit dem nun bekannten Anfangsvektor {), die Rechnung; in der Tat wird w,) = 0, 
M,) = 0, wie es sein muf. 


4. Jetzt gehen wir zu den abgeschnittenen Zweigen zuriick. Fiir die beiden Zahlenschemata 2.1 
und 2.11 brauchen wir als erstes die Werte us, @y baw. ug, Yy Aus dem Schema 2.III bekommen 
wir nach (2.20) zunachst 


Wg = Us = 0,355; Go = 9; = — 0,193;, (2.23) 
—w, = ug =— 1,448; , =H, = 3,680, (2.24) 

weiter aus (2.09) und (2.17) 
( = — 0,0156 (— 0,193) + 0,2344 (0,355) = 0,0863 , (2.25) 
G4 = 0,0174 (3,680) — 0,355 - (— 1,448) = 0,5777,, (2.26) 


berechnen damit alle Werte p;, M; in 2.1 und 2.11 und ebenso die von den Pfosten in die waage- 
rechten Trager eingeleiteten [in (2.06) und (2.14) eingeklammerten] Kraftgré8en. Da diese aber 
gleich der im AnschluBpunkt reduzierten Gesamtbelastung eines jeden Pfosten sind, findet man 
deren Reaktionen Q, und M, in den Einspannungen bzw. Gelenken und daraus Querkrafte und 
Biegemomente der Pfosten nach Abb. 17 auf folgende Weise: 


0,=@,—9, M,=M,—(,—Qh—M, (2.27) 
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zum Beispiel fiir die Pfosten 11, 12 und 16 


On = 0, = 86 T S04 = 18 539° } 
Qi. = Q2,—9 =— 1,614, 
ioe = Vig, sa == 15.574 5 
WE My (0, — 24) 1 — 12° 0 8,467 1. 12 = 3,539, re) 
EN en = 07278 11614 S 0:893 
Mig, = Mig, — Qig (2 = 7,787 — 15,574/2 = 0. J 
SchlieBlich gibt die Gleichung (20), die hier infolge der fehlenden Belastung iibergeht in 
Q, = const = 3 (M;—M,.,), (2.29) 
RS 
SB 
Momente in 2¢q0% ‘nn 5514 


Krafre in ql = 


Ss 


) 1614 y 1781 ) 3104 ; 

> > 

q 3,533 | 0,893 i 0.949 | 145 
S Ss 


Abb, 18. Gleichgewichtssystem am Tragwerk von Abb. 13. 


4907 15,874 73,908 


Q in P*=24g0 


75,533 1614 4,787 5104 
Abb. 19. Querkraftverlauf fiir das Tragwerk von Abb. 13. 


den Querkraftverlauf der Trager I und II und damit auch die senkrechten Auflagerkrafte, d. h. 
die Langskrafte in den Pfosten. Anhand der Abb. 18 bestatigt man leicht, daB die Gleichgewichts- 
bedingungen erfiillt sind. Die Abb. 19 und 20 zeigen Querkraft- und Momentenverlauf, und zwar 
in den Dimensionen P* = q 1/24 und P* I* = q [7/24 fiir beliebige Werte von q und /. Die Dimen- 
sionen der Drehwinkel und Verschiebungen sind 1 und /* =/. Nach (2.01) miissen wir diese 
Werte bei der gegebenen Belastung noch mit 11/480 = 0,0229 multiplizieren. Beispielsweise 
hat die waagerechte Verschiebung des Zweiges I, namlich us = we = 0,355 1* (2.23) die wahre 
GréBe 
0,0229 - 0,355 - 200 cm = 1,63 cm. 
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. Sollen nun die Zweige I und II in B und C gelenkig angeschlossen werden, so ist in den 
E BREN: (2.11) und (2.19) M, = 0 bzw. M, = 0 zu setzen. Addiert man die mit passenden 
Vielfachen multiplizierten ersten Zeilen zu den zweiten, so daB ps baw. (pe herausfallt, so wird 


0 0,422, = "1828 (“s) aye 0,527 By cale ) (2.30) 
(r,—9) = = 10,731 0,0000) \ps/’? \Lg!  \— 28,157 0,000) \qe)’ 


also nach (2.20) 


ra 0,422 
L, = 0g = 9 = 10,731 ws, v3 = (F Gs) = 1309 (3 = 0,231 wg, 
* 0,517 (2.31) 
— Db, =0,=—28,157 W,, 93 = Go) =—Te78 9 = — 0,275 We . 


73,173 


M in P*1*= 24g? 


5533 0.893 O49 41S 


Abb. 20. Biegemomentenverlauf fiir das Tragwerk von Abb. 13. 


M, 0; 1 
0 0 0 0 
Zahlenschema 2.111 0 0 0 0 
| 0 0 17452, = te 
0 1 0 aa 
0 0 1 
Vat, oe ibn 01 20 0 0,5 0,167 0 0,355 
i Sh i 0 —] —=0,5 0 —0,193 
ee 08a le One 0:22 1428 3,039 1,984 0 —0,857 |= fg, 
OP 0 00 He 595 = 105828054229] | 830 oe 1,968 
—l 0 —1/2—5/3 —8 10,506 1,406 0 0 1 a ny 
Wel 1/2. 6s 1 20 0 2,047 1,490 1,448 
he aah SS 2) 0 —=1,121 —1,977 a 5 3,680 
oe ee. ee Ome el L 12, —0,517—1,878 |] —1,750 3,910 20,646 5,453 labs. 
0 0 © 1 24 =-28,200" 0517 |= 63-181- 42.147 a8 489 ie ,907 
—l1 0 —1/2 —5/3—40 —28,816 2,395 0 0 
1 —1/2 1/8 1/48 1/16 , [1,072 2,089 4,612 0,000 
0 1-1/2 —1/8.,1/2 Be Bee | Os: 3,566 
ig ae wee Oa ees |—33,340 —17,163 4,902 0,000 |= no: 
Q. vO one] 19% (16 hurl) onthe Aire ek te oe oe —4,907 
11/2 5/8 1.396" ©. 0 0 irr 


Man hat daher die Federmatrizen (2.21) in As baw. X'y zu ersetzen durch 


0 oer 0 0 
10/73 ey On emcal ase 0) (282) 
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Die beiden Koppelfedern in B und C (Abb. 16) sind somit in gewdhnliche horizontal wirkende 
Druckfedern entartet. Von Formel (2.21) an ist die Rechnung mit neuen Zahlenwerten zu wieder- 
holen. Fiir die Praxis besonders angenehm ist der Umstand, daB die beiden Federmatrizen (2.11) 
und (2.19) von zwei Rechnern gleichzeitig und unabhangig ermittelt und kontrolliert werden 
kénnen. Von der absoluten Uberlegenheit der hier gezeigten Methode iiberzeugt man sich am 
besten, wenn man die Aufgabe nach dem iiblichen Eliminationsverfahren zu rechnen versucht; 
das Tragwerk ist 18fach statisch, bzw. 13fach geometrisch iiberbestimmt. 


c) Drittes Beispiel. Gegeben ist das symmetrische und symmetrisch belastete offene Trag- 
werk von Abb. 21. Alle Stabe seien in Langsrichtung starr. Man ermittle Querkraft-, Liangs- 
kraft- und Biegemomentenverlauf. 


Abb, 21. Offener symmetrischer Rahmen mit symmetrischer Belastung, dazu Ersatztriger. 


Aus Symmetriegriinden kénnen wir uns auf die linke Halfte des Tragwerkes beschranken. Die 
Bedingungen im Mittelschnitt 2 lauten dann (H und V sind Horizontal- bzw. Vertikalkrafte) 


C= 0. (3.01) 
x = WW, s6in Xx -+- U,cosa = 0 > w, =— u, cotga, (3.02) 
Vz, =Q, cosa —L, sina =0 —> Q, = (L,/cos a) sina = H, sina. (3.03) 
Infolge der Langsstarrheit ist uy = u, = u, = u, Uy = W, = Uta, WwW, = —ucotg a, also 
. “y 10 
Wy— UW, = 0, w,— un, = utga + uctga—u (3 3) = zu. (3.04) 


Ebenso wie im zweiten Beispiel unterdriicken wir die Paaie (L,u) und (Q,w), ersetzen die beiden 
Pfosten I und II durch Koppelfedern und zerlegen die senkrecht wirkende Belastung in ihre beiden 
Anteile q cos?.« und qg sina cosa. Mit den VergleichsgréBen 


2EJ 
ee ae rE eran jee dre bad 
wird 
Yeosa=1, EJ=12, gq=l2; (3.05) 
ferner ist 
sin x = 1/10 = 0,3162 , cosa = 3/}/ 10 «= 0,9487. (3.06) 


Fiir die Lastspalte der Matrix (18) brauchen wir jetzt die Werte der Tabelle 1, 3 


3 10 


12 - (1/2) 1 BST 10 
1/2 3 


= fa a = 3 u=20u 
casi eee 24 - 1/2 8? , 


3 
lap 
Wy—w, = 9, 7 (1 — u 2) 
2 


(3.07) 


ay * 72 FS 2 3 9 les 
M,=1* =3, M,=VCPX 2, —6ES (mj —w)/f=— zu =—40u. 
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Mit den Langen und Biegesteifigkeiten der beiden Pfosten 


EJ; =15L=0,75, h=1,51=1,3cose-1* =1,5-0,9487 = 1,4230, 
i 11 W (No 
EJu =105=050, n= l= Zcosa+1* = 0,9487 = 1,7392 


bekommen wir als in den Haupttrager (Abb. 21 rechts) eingeleitete Biegemomente M; und Hori- 
zontalkrafte H, mit x, = u/cos« und nach Tabelle 2, la 


iy ye oan ze 2, Las %% ; — 2,3424 u — 2,1082 | ee) 
noe TRIO eet ee, 
sel AL peep es — 1,0454 u — 1,1499 g, 

ae Lge sue eee = (3.10) 
ae 1,7302 cosa 17392 7 eee ease 


| 


kKratte in gt/12 
Momente in gt te 
Abb, 22, Gleichgewichtssystem am Rahmen von Abb, 21. 


Die Drehfederkonstanten des Haupttragers sind demnach C, = 2,1082 und C, = 1,1499. Damit 
laBt sich das Zahlenschema ausfiillen. Die eingerahmte Bedingung ~, = 0 gibt 


Py = — 1,9632 u — 0,2784,, (3.11) 


und somit die Spalte (a) des Schemas. Als nachstes addieren wir die zweiten Spalten von (3.09) 


und (3.10): 


. H, + H, = H, = —4,4943 u — 2,2222 ¢, — 0,9917 @, 
= — 4,4943 u — 2,2222 (— 1,9632 u — 0,2784) — 0,9917 (5,7228 u + 0,1437), | (3.12) 
H, = — 5,8071 u + 0,4762 . 


Weiter wird nach (20) 


eae ae vs — MM 4 6, =[(28,5536 u —0,0394) — (— 17,1084 u -+ 0,8258) —1,5] 2 +6, 


Q, = 91,3240 u + 1.2696 , 


(3.13) 


und aus der Bedingung (3.03) 
Q, = 91,3240 u + 1,2696 = H, sina = (— 5,8071 + 0,4762) 0,3162 


folgt nun 


u = —0,0120. (3.14) 
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In ppee ciepalie (b) des Schemas finden wir jetzt die Biegemomente in der Dimension P* 1* 
qe / 12 =a P12 vor, berechnen daraus den Querkraftverlauf nach (20), weiter aus den nun 
zahlenmafig vorliegenden eingeleiteten KraftgréBen (3.09) und (3.10) Querkraft und Biegemoment 


5,608! 


44 
Lingle 1,3888 os 


72,3479 Druck 


Q in gh/t2 


56520 Druck 


Abb. 23. Querkraftverlauf fiir die linke Halfte des Rahmens 


Sea Wa eOT: Abb. 24, Lingskraftverlauf fiir die linke Halfte des Rahmens 


von Abb. 21, 


1048 


038i? 


M in ght 


2969 00308 


Abb. 25, Biegemomentenverlauf fiir die linke Hilfte des Rahmens von Abb. 21. 


(a) (b) 
Po t a 1 
Zahlenschema 3 1 0 —1,9632 —0,2784\ /—0,255 
( 29,1082 —-2,3424 wy 1,7964 | 565) 
0 1 0 il 1 


—l 1 —1+1,0454u, 1,1499 0 1 0 I 1 


= 4055 -0,1254-20'u | 14,1866 27,8513 u + 3,950] 0,0002  0,0005 0,000 
ig en 0715+ — 40 ‘ 28,5536 0.20 [0.20 
0 0 


Pee 1,5 0,3755-20 


—2 1 —]J ia) —4,1082 —2,3424u — 1 ( 53,7228 =: 0, 1437 0,075 
3 —2 3—1,0454 u 11199 11,9404 +6,3329 wu + 4,1499 Keer te 0,8258 1,031 


0 1 1 


_ 2 * a + " 
ce . . - 
Z =) “ , 
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der Pfosten und schlieBlich durch eine einfache Zerlegung in den AnschluBpunkten 0 und | auch 
die Langskrafte. (Abb. 22 bis 25.) Fiir den rechten Teil des Tragwerkes sind Momente und Langs- 
krafte symmetrisch, die Querkrafte dagegen antisymmetrisch fortzusetzen. 


12. Zusammenfassung und Ausblick. Wir haben gezeigt, wie sich mit Hilfe des Begriffes der 
Koppelfeder jeder offene Rahmen auf einen (oder mehrere) Durchlauftrager zuriickfiihren 1aBt. 
Die Zahlenrechnung wurde im Hinblick auf Kontrollen und leichte Programmierbarkeit fiir Rechen- 
automaten noch weiter schematisiert durch Einfiihrung von rechteckigen erweiterten Leitmatrizen 
und einer zweckmaBig definierten sogenannten ,,verschrankten“ Matrizenmultiplikation. Drei 
Beispiele zeigen den praktischen Gang des Verfahrens. In einer dritten und abschlieBenden Arbeit 
sollen Tragwerke mit geschlossenen Umlaufen (Stockwerkrahmen, Tragerroste, Holme usw.) be- 
handelt werden, wobei die ganze Reichweite und Anpassungsfahigkeit des Reduktionsverfahrens 
noch mehr als bisher hervortreten wird. 


(Eingegangen am 8. April 1957) 
Anschrift des Verfassers: Dr.-Ing. Sigurd Falk, Braunschweig, Celler Str. 93a 


Berichtigung 
zu meiner Arbeit im Band XXV, S. 164 des Ingenieur-Archivs: ,,Biegetheorie 
der Rotationsschale mit flacher, kreisférmiger Erzeugender* 


Von R. Windels 


Die erste Gleichung (14) muB (mit Vertauschung der Indices g und « in der Klammer) lauten 


E 


n= pap (Cae + Hep) aa 


Die Gleichung (56) mu [mit Streichung des Faktors = vor 20) lauten 


ne Be i a eat heacis 
heat VWiew+a ae Qo)? — a EO (56) 


(Eingegangen am 3. Dezember 1957.) 


Anschrift des Verfassers: Dr.-Ing. R. Windels, Hamburg 36, Jungfernstieg 49. 
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Theorie schallnaher Strémungen 


Von 
Dr.-Ing. habil. K. G. Guderley 


Mit 125 Abbildungen 
XV, 376 Seiten Gr.-8°. 1957 


Ganzleinen DM 42,— 


Inhaltstibersicht: I. Allgemeine Grundlagen. — II. Vereinfachung der Strémungsdifferenzialgleichungen, 
Abnlichkeitsgesetz fiir schallnahe Strémungen. — ITI. Linearisierte Behandlung schallnaher Strémungen. — 
IV. Exakte Lésungen der Potentialgleichung fiir schallnahe Stromungen. — V. Die Grundlagen der Hodo- 
graphenmethode. — VI. Diskussion schallnaher Strémungsfelder mit Hilfe des Hodographen. — VII. Par- 
tikularlésungen der Tricomischen Gleichung.— VIII. Strémungen mit der Macuschen Zahl 1. — IX, Stré- 
mungsfelder, die nur wenig von einer Strémung mit der Macu-Zahl 1 abweichen. — X. Hinzeluntersuchun- 
gen, die Partikularlésungen der Form Gl. (VII 3.3) benutzen. — XI. Achsensymmetrische Strémungen. — 


Literaturverzeichnis, — Sachverzeichnis. 


Das Verhalten von Gasstrémungen zeigt weitgehende Verschiedenheiten, je nachdem die Strémungs- 
geschwindigkeit kleiner oder gréBer als die Schallgeschwindigkeit ist. Der tinersats von einer Form zur 
anderen, das ,,schallnahe“ Gebiet, bereitet fiir das Verstaéndnis besondere Schwierigkeiten, sowohl vom 
physikalischen als vom mathematischen Standpunkt aus. Jedoch ist in den letzten zwélf Jahren auch hier 
ein gesicherter Kern der Erkenntnis geschaffen worden, Das vorliegende Buch stellt erstmalig die Theorie 
schallnaher Strémungen im Zusammenhang dar. Die Forschung auf diesem Gebiete ist durch eine be- 
sonders enge Wechselwirkung zwischen physikalischem und mathematischem Denken gekennzeichnet. 
Das bestimmt auch den Charakter dieses Buches. Es zeigt die mathematische Theorie, wie sie og physi- 
kalischen Vorstellungen herauswichst. Dabei fiihrt es auch in die mathematischen Fragen ein, die hier 
entstehen und die gelegentlich tiber die klassischen Fragenstellungen hinausgehen. Doch wird auch dann 
der Versuch gemacht, durch physikalische Vergleiche die mathematischen Ergebnisse anschaulich zu 
machen. Die Darstellung ist fiir den Aerodynamiker, der einige mathematische Erfahrungen besitzt, be- 


stimmt; aber auch dem Mathematiker kann das Buch als Einfiihrung in die Probleme des schallnahen 


Gebiets yon Wert sein. 
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Springer-Verlag + Berlin / Géttingen / Heidelberg. Printed in Germany. V/12/6 0,8 K.B./Z, 040, 
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